THE CONTROL VARIATIONAL METHOD FOR ELASTIC CONTACT PROBLEMS

Directory of Open Access Journals (Sweden)

Mircea Sofonea

2010-07-01

Full Text Available We consider a multivalued equation of the form Ay + F(y = fin a real Hilbert space, where A is a linear operator and F represents the (Clarke subdifferential of some function. We prove existence and uniqueness results of the solution by using the control variational method. The main idea in this method is to minimize the energy functional associated to the nonlinear equation by arguments of optimal control theory. Then we consider a general mathematical model describing the contact between a linearly elastic body and an obstacle which leads to a variational formulation as above, for the displacement field. We apply the abstract existence and uniqueness results to prove the unique weak solvability of the corresponding contact problem. Finally, we present examples of contact and friction laws for which our results work.

Application of an enriched FEM technique in thermo-mechanical contact problems

Science.gov (United States)

Khoei, A. R.; Bahmani, B.

2018-02-01

In this paper, an enriched FEM technique is employed for thermo-mechanical contact problem based on the extended finite element method. A fully coupled thermo-mechanical contact formulation is presented in the framework of X-FEM technique that takes into account the deformable continuum mechanics and the transient heat transfer analysis. The Coulomb frictional law is applied for the mechanical contact problem and a pressure dependent thermal contact model is employed through an explicit formulation in the weak form of X-FEM method. The equilibrium equations are discretized by the Newmark time splitting method and the final set of non-linear equations are solved based on the Newton-Raphson method using a staggered algorithm. Finally, in order to illustrate the capability of the proposed computational model several numerical examples are solved and the results are compared with those reported in literature.

Problems in nonlinear resistive MHD

International Nuclear Information System (INIS)

Turnbull, A.D.; Strait, E.J.; La Haye, R.J.; Chu, M.S.; Miller, R.L.

1998-01-01

Two experimentally relevant problems can relatively easily be tackled by nonlinear MHD codes. Both problems require plasma rotation in addition to the nonlinear mode coupling and full geometry already incorporated into the codes, but no additional physics seems to be crucial. These problems discussed here are: (1) nonlinear coupling and interaction of multiple MHD modes near the B limit and (2) nonlinear coupling of the m/n = 1/1 sawtooth mode with higher n gongs and development of seed islands outside q = 1

Primal and Dual Penalty Methods for Contact Problems with Geometrical Non-linearities

Czech Academy of Sciences Publication Activity Database

Vondrák, V.; Dostál, Z.; Dobiáš, Jiří; Pták, Svatopluk

-, č. 5 (2005), s. 449-450 ISSN 1617-7061. [GAMM Annual Meeting 2005. Luxembourg, 28.03.2005-01.04.2005] R&D Projects: GA ČR(CZ) GA101/05/0423 Institutional research plan: CEZ:AV0Z20760514 Keywords : primal penalty * dual penalty * contact problem Subject RIV: BA - General Mathematics

Nonlinear Displacement Discontinuity Model for Generalized Rayleigh Wave in Contact Interface

Energy Technology Data Exchange (ETDEWEB)

Kim, No Hyu; Yang, Seung Yong [Korea University of Technology and Education, Cheonan (Korea, Republic of)

2007-12-15

Imperfectly jointed interface serves as mechanical waveguide for elastic waves and gives rise to two distinct kinds of guided wave propagating along the interface. Contact acoustic nonlinearity (CAN) is known to plays major role in the generation of these interface waves called generalized Rayleigh waves in non-welded interface. Closed crack is modeled as non-welded interface that has nonlinear discontinuity condition in displacement across its boundary. Mathematical analysis of boundary conditions and wave equation is conducted to investigate the dispersive characteristics of the interface waves. Existence of the generalized Rayleigh wave(interface wave) in nonlinear contact interface is verified in theory where the dispersion equation for the interface wave is formulated and analyzed. It reveals that the interface waves have two distinct modes and that the phase velocity of anti-symmetric wave mode is highly dependent on contact conditions represented by linear and nonlinear dimensionless specific stiffness

Nonlinear Displacement Discontinuity Model for Generalized Rayleigh Wave in Contact Interface

International Nuclear Information System (INIS)

Kim, No Hyu; Yang, Seung Yong

2007-01-01

Imperfectly jointed interface serves as mechanical waveguide for elastic waves and gives rise to two distinct kinds of guided wave propagating along the interface. Contact acoustic nonlinearity (CAN) is known to plays major role in the generation of these interface waves called generalized Rayleigh waves in non-welded interface. Closed crack is modeled as non-welded interface that has nonlinear discontinuity condition in displacement across its boundary. Mathematical analysis of boundary conditions and wave equation is conducted to investigate the dispersive characteristics of the interface waves. Existence of the generalized Rayleigh wave(interface wave) in nonlinear contact interface is verified in theory where the dispersion equation for the interface wave is formulated and analyzed. It reveals that the interface waves have two distinct modes and that the phase velocity of anti-symmetric wave mode is highly dependent on contact conditions represented by linear and nonlinear dimensionless specific stiffness

Scalable algorithms for contact problems

CERN Document Server

Dostál, Zdeněk; Sadowská, Marie; Vondrák, Vít

2016-01-01

This book presents a comprehensive and self-contained treatment of the authors’ newly developed scalable algorithms for the solutions of multibody contact problems of linear elasticity. The brand new feature of these algorithms is theoretically supported numerical scalability and parallel scalability demonstrated on problems discretized by billions of degrees of freedom. The theory supports solving multibody frictionless contact problems, contact problems with possibly orthotropic Tresca’s friction, and transient contact problems. It covers BEM discretization, jumping coefficients, floating bodies, mortar non-penetration conditions, etc. The exposition is divided into four parts, the first of which reviews appropriate facets of linear algebra, optimization, and analysis. The most important algorithms and optimality results are presented in the third part of the volume. The presentation is complete, including continuous formulation, discretization, decomposition, optimality results, and numerical experimen...

Chaotic behaviour in the non-linear optimal control of unilaterally contacting building systems during earthquakes

International Nuclear Information System (INIS)

Liolios, A.A.; Boglou, A.K.

2003-01-01

The paper presents a new numerical approach for a non-linear optimal control problem arising in earthquake civil engineering. This problem concerns the elastoplastic softening-fracturing unilateral contact between neighbouring buildings during earthquakes when Coulomb friction is taken into account under second-order instabilizing effects. So, the earthquake response of the adjacent structures can appear instabilities and chaotic behaviour. The problem formulation presented here leads to a set of equations and inequalities, which is equivalent to a dynamic hemivariational inequality in the way introduced by Panagiotopoulos [Hemivariational Inequalities. Applications in Mechanics and Engineering, Springer-Verlag, Berlin, 1993]. The numerical procedure is based on an incremental problem formulation and on a double discretization, in space by the finite element method and in time by the Wilson-θ method. The generally non-convex constitutive contact laws are piecewise linearized, and in each time-step a non-convex linear complementarity problem is solved with a reduced number of unknowns

Chaotic behaviour in the non-linear optimal control of unilaterally contacting building systems during earthquakes

CERN Document Server

Liolios, A

2003-01-01

The paper presents a new numerical approach for a non-linear optimal control problem arising in earthquake civil engineering. This problem concerns the elastoplastic softening-fracturing unilateral contact between neighbouring buildings during earthquakes when Coulomb friction is taken into account under second-order instabilizing effects. So, the earthquake response of the adjacent structures can appear instabilities and chaotic behaviour. The problem formulation presented here leads to a set of equations and inequalities, which is equivalent to a dynamic hemivariational inequality in the way introduced by Panagiotopoulos [Hemivariational Inequalities. Applications in Mechanics and Engineering, Springer-Verlag, Berlin, 1993]. The numerical procedure is based on an incremental problem formulation and on a double discretization, in space by the finite element method and in time by the Wilson-theta method. The generally non-convex constitutive contact laws are piecewise linearized, and in each time-step a non-c...

Nonlinear dynamic modeling of rotor system supported by angular contact ball bearings

Science.gov (United States)

Wang, Hong; Han, Qinkai; Zhou, Daning

2017-02-01

In current bearing dynamic models, the displacement coordinate relations are usually utilized to approximately obtain the contact deformations between the rolling element and raceways, and then the nonlinear restoring forces of the rolling bearing could be calculated accordingly. Although the calculation efficiency is relatively higher, the accuracy is lower as the contact deformations should be solved through iterative analysis. Thus, an improved nonlinear dynamic model is presented in this paper. Considering the preload condition, surface waviness, Hertz contact and elastohydrodynamic lubrication, load distribution analysis is solved iteratively to more accurately obtain the contact deformations and angles between the rolling balls and raceways. The bearing restoring forces are then obtained through iteratively solving the load distribution equations at every time step. Dynamic tests upon a typical rotor system supported by two angular contact ball bearings are conducted to verify the model. Through comparisons, the differences between the nonlinear dynamic model and current models are also pointed out. The effects of axial preload, rotor eccentricity and inner/outer waviness amplitudes on the dynamic response are discussed in detail.

On a Mixed Nonlinear One Point Boundary Value Problem for an Integrodifferential Equation

Directory of Open Access Journals (Sweden)

Mesloub Said

2008-01-01

Full Text Available This paper is devoted to the study of a mixed problem for a nonlinear parabolic integro-differential equation which mainly arise from a one dimensional quasistatic contact problem. We prove the existence and uniqueness of solutions in a weighted Sobolev space. Proofs are based on some a priori estimates and on the Schauder fixed point theorem. we also give a result which helps to establish the regularity of a solution.

Nonlinear vibration of rectangular atomic force microscope cantilevers by considering the Hertzian contact theory

Energy Technology Data Exchange (ETDEWEB)

Sadeghi, A., E-mail: a_sadeghi@srbiau.ac.ir [Islamic Azad Univ., Dept. of Mechanical and Aerospace Engineering, Science and Research Branch, Tehran (Iran, Islamic Republic of); Zohoor, H. [Sharif Univ. of Technology, Center of Excellence in Design, Robotics and Automation, Tehran (Iran, Islamic Republic of); The Academy of Sciences if I.R. Iran (Iran, Islamic Republic of)

2010-05-15

The nonlinear flexural vibration for a rectangular atomic force microscope cantilever is investigated by using Timoshenko beam theory. In this paper, the normal and tangential tip-sample interaction forces are found from a Hertzian contact model and the effects of the contact position, normal and lateral contact stiffness, tip height, thickness of the beam, and the angle between the cantilever and the sample surface on the nonlinear frequency to linear frequency ratio are studied. The differential quadrature method is employed to solve the nonlinear differential equations of motion. The results show that softening behavior is seen for most cases and by increasing the normal contact stiffness, the frequency ratio increases for the first mode, but for the second mode, the situation is reversed. The nonlinear-frequency to linear-frequency ratio increases by increasing the Timoshenko beam parameter, but decreases by increasing the contact position for constant amplitude for the first and second modes. For the first mode, the frequency ratio decreases by increasing both of the lateral contact stiffness and the tip height, but increases by increasing the angle α between the cantilever and sample surface. (author)

Introduction to nonlinear finite element analysis

CERN Document Server

Kim, Nam-Ho

2015-01-01

This book introduces the key concepts of nonlinear finite element analysis procedures. The book explains the fundamental theories of the field and provides instructions on how to apply the concepts to solving practical engineering problems. Instead of covering many nonlinear problems, the book focuses on three representative problems: nonlinear elasticity, elastoplasticity, and contact problems. The book is written independent of any particular software, but tutorials and examples using four commercial programs are included as appendices: ANSYS, NASTRAN, ABAQUS, and MATLAB. In particular, the MATLAB program includes all source codes so that students can develop their own material models, or different algorithms. This book also: ·         Presents clear explanations of nonlinear finite element analysis for elasticity, elastoplasticity, and contact problems ·         Includes many informative examples of nonlinear analyses so that students can clearly understand the nonlinear theory ·    ...

Nonlinear dynamics of contact interaction of a size-dependent plate supported by a size-dependent beam

Science.gov (United States)

Awrejcewicz, J.; Krysko, V. A.; Yakovleva, T. V.; Pavlov, S. P.; Krysko, V. A.

2018-05-01

A mathematical model of complex vibrations exhibited by contact dynamics of size-dependent beam-plate constructions was derived by taking the account of constraints between these structural members. The governing equations were yielded by variational principles based on the moment theory of elasticity. The centre of the investigated plate was supported by a beam. The plate and the beam satisfied the Kirchhoff/Euler-Bernoulli hypotheses. The derived partial differential equations (PDEs) were reduced to the Cauchy problems by the Faedo-Galerkin method in higher approximations, whereas the Cauchy problem was solved using a few Runge-Kutta methods. Reliability of results was validated by comparing the solutions obtained by qualitatively different methods. Complex vibrations were investigated with the help of methods of nonlinear dynamics such as vibration signals, phase portraits, Fourier power spectra, wavelet analysis, and estimation of the largest Lyapunov exponents based on the Rosenstein, Kantz, and Wolf methods. The effect of size-dependent parameters of the beam and plate on their contact interaction was investigated. It was detected and illustrated that the first contact between the size-dependent structural members implies chaotic vibrations. In addition, problems of chaotic synchronization between a nanoplate and a nanobeam were addressed.

Nonlinear Elliptic Boundary Value Problems at Resonance with Nonlinear Wentzell Boundary Conditions

Directory of Open Access Journals (Sweden)

Ciprian G. Gal

2017-01-01

Full Text Available Given a bounded domain Ω⊂RN with a Lipschitz boundary ∂Ω and p,q∈(1,+∞, we consider the quasilinear elliptic equation -Δpu+α1u=f in Ω complemented with the generalized Wentzell-Robin type boundary conditions of the form bx∇up-2∂nu-ρbxΔq,Γu+α2u=g on ∂Ω. In the first part of the article, we give necessary and sufficient conditions in terms of the given functions f, g and the nonlinearities α1, α2, for the solvability of the above nonlinear elliptic boundary value problems with the nonlinear boundary conditions. In other words, we establish a sort of “nonlinear Fredholm alternative” for our problem which extends the corresponding Landesman and Lazer result for elliptic problems with linear homogeneous boundary conditions. In the second part, we give some additional results on existence and uniqueness and we study the regularity of the weak solutions for these classes of nonlinear problems. More precisely, we show some global a priori estimates for these weak solutions in an L∞-setting.

A nonlinear oscillatory problem

International Nuclear Information System (INIS)

Zhou Qingqing.

1991-10-01

We have studied the nonlinear oscillatory problem of orthotropic cylindrical shell, we have analyzed the character of the oscillatory system. The stable condition of the oscillatory system has been given. (author). 6 refs

Study on the contact ratio of base mat of reactor buildings considering nonlinear soil-structure interaction effects

International Nuclear Information System (INIS)

Aihara, S.; Atsumi, K.; Ujiie, K.; Emori, K.; Odajima, M.; Masuda, K.

1983-01-01

The objective of this paper is to evaluate the nonlinear soil-structure interaction effects resulting from base mat uplift for static lateral loads. Nonlinear soil-structure interaction effects are modeled through the use of equivalent soil-structure interaction frictional and axial springs, which properties are determined by results of experimental data. It is assumed that normal stresses in compression and corresponding shear stresses, and friction, can occur in the area of contact between the embedded structure and soil. The remaining parts of the structure and soil are based on elastic analysis. A two-dimensional finite element method with incremental loadings is applied. The substructuring technique is used to reduce computation time. The results of this method with respect to the contact ratio of the base mat are compared with the values obtained by static elastic calculation which is simply derived from an overturning moment and a vertical load of the structure. This analytical concept will be developed into dynamic problems, and then it will be possible to state whether or not this concept can represent a true alternative for the contact ratio of the base mat of a structure. (orig./HP)

Overlapping Schwarz for Nonlinear Problems. An Element Agglomeration Nonlinear Additive Schwarz Preconditioned Newton Method for Unstructured Finite Element Problems

Energy Technology Data Exchange (ETDEWEB)

Cai, X C; Marcinkowski, L; Vassilevski, P S

2005-02-10

This paper extends previous results on nonlinear Schwarz preconditioning ([4]) to unstructured finite element elliptic problems exploiting now nonlocal (but small) subspaces. The non-local finite element subspaces are associated with subdomains obtained from a non-overlapping element partitioning of the original set of elements and are coarse outside the prescribed element subdomain. The coarsening is based on a modification of the agglomeration based AMGe method proposed in [8]. Then, the algebraic construction from [9] of the corresponding non-linear finite element subproblems is applied to generate the subspace based nonlinear preconditioner. The overall nonlinearly preconditioned problem is solved by an inexact Newton method. Numerical illustration is also provided.

Mathematical Modeling of Contact Problems of Elasticity Theory with Continuous Unilateral Contact

Directory of Open Access Journals (Sweden)

I. V. Stankevich

2015-01-01

Full Text Available The work [1] presents the formulation and numerical solution of the problem concerning the unilateral discrete contact interaction of an elastic body and a rigid half-space. However, many parts and components of engineering structures have a pronounced continuous contact within a given surface [2, 3]. In this paper we consider a special case of this option of contact interaction when, the elastic body of finite size, subjected to external forces, is based on a rigid half-space. Contact occurs through a dedicated contact surface, which in general can change their sizes.Developed to solve this problem, a numerical algorithm is a further adaptation and development of the approaches described in [1]. The paper shows results of solving the model problem of the elasticity theory with and without taking friction into account. In the latter case, were additionally obtained numerical data characterizing the convergence of the solution.

Multisplitting for linear, least squares and nonlinear problems

Energy Technology Data Exchange (ETDEWEB)

Renaut, R.

1996-12-31

In earlier work, presented at the 1994 Iterative Methods meeting, a multisplitting (MS) method of block relaxation type was utilized for the solution of the least squares problem, and nonlinear unconstrained problems. This talk will focus on recent developments of the general approach and represents joint work both with Andreas Frommer, University of Wupertal for the linear problems and with Hans Mittelmann, Arizona State University for the nonlinear problems.

Multigrid Reduction in Time for Nonlinear Parabolic Problems

Energy Technology Data Exchange (ETDEWEB)

Falgout, R. D. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Manteuffel, T. A. [Univ. of Colorado, Boulder, CO (United States); O' Neill, B. [Univ. of Colorado, Boulder, CO (United States); Schroder, J. B. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2016-01-04

The need for parallel-in-time is being driven by changes in computer architectures, where future speed-ups will be available through greater concurrency, but not faster clock speeds, which are stagnant.This leads to a bottleneck for sequential time marching schemes, because they lack parallelism in the time dimension. Multigrid Reduction in Time (MGRIT) is an iterative procedure that allows for temporal parallelism by utilizing multigrid reduction techniques and a multilevel hierarchy of coarse time grids. MGRIT has been shown to be effective for linear problems, with speedups of up to 50 times. The goal of this work is the efficient solution of nonlinear problems with MGRIT, where efficient is defined as achieving similar performance when compared to a corresponding linear problem. As our benchmark, we use the p-Laplacian, where p = 4 corresponds to a well-known nonlinear diffusion equation and p = 2 corresponds to our benchmark linear diffusion problem. When considering linear problems and implicit methods, the use of optimal spatial solvers such as spatial multigrid imply that the cost of one time step evaluation is fixed across temporal levels, which have a large variation in time step sizes. This is not the case for nonlinear problems, where the work required increases dramatically on coarser time grids, where relatively large time steps lead to worse conditioned nonlinear solves and increased nonlinear iteration counts per time step evaluation. This is the key difficulty explored by this paper. We show that by using a variety of strategies, most importantly, spatial coarsening and an alternate initial guess to the nonlinear time-step solver, we can reduce the work per time step evaluation over all temporal levels to a range similar with the corresponding linear problem. This allows for parallel scaling behavior comparable to the corresponding linear problem.

Nonlinear acceleration of transport criticality problems

International Nuclear Information System (INIS)

Park, H.; Knoll, D.A.; Newman, C.K.

2011-01-01

We present a nonlinear acceleration algorithm for the transport criticality problem. The algorithm combines the well-known nonlinear diffusion acceleration (NDA) with a recently developed, Newton-based, nonlinear criticality acceleration (NCA) algorithm. The algorithm first employs the NDA to reduce the system to scalar flux, then the NCA is applied to the resulting drift-diffusion system. We apply a nonlinear elimination technique to eliminate the eigenvalue from the Jacobian matrix. Numerical results show that the algorithm reduces the CPU time a factor of 400 in a very diffusive system, and a factor of 5 in a non-diffusive system. (author)

Three-Field Modelling of Nonlinear Nonsmooth Boundary Value Problems and Stability of Differential Mixed Variational Inequalities

Directory of Open Access Journals (Sweden)

J. Gwinner

2013-01-01

Full Text Available The purpose of this paper is twofold. Firstly we consider nonlinear nonsmooth elliptic boundary value problems, and also related parabolic initial boundary value problems that model in a simplified way steady-state unilateral contact with Tresca friction in solid mechanics, respectively, stem from nonlinear transient heat conduction with unilateral boundary conditions. Here a recent duality approach, that augments the classical Babuška-Brezzi saddle point formulation for mixed variational problems to twofold saddle point formulations, is extended to the nonsmooth problems under consideration. This approach leads to variational inequalities of mixed form for three coupled fields as unknowns and to related differential mixed variational inequalities in the time-dependent case. Secondly we are concerned with the stability of the solution set of a general class of differential mixed variational inequalities. Here we present a novel upper set convergence result with respect to perturbations in the data, including perturbations of the associated nonlinear maps, the nonsmooth convex functionals, and the convex constraint set. We employ epiconvergence for the convergence of the functionals and Mosco convergence for set convergence. We impose weak convergence assumptions on the perturbed maps using the monotonicity method of Browder and Minty.

Non-linear frequency and amplitude modulation of a nano-contact spin torque oscillator

OpenAIRE

Muduli, P. K.; Pogoryelov, Ye.; Bonetti, S.; Consolo, G.; Mancoff, Fred; Åkerman, Johan

2009-01-01

We study the current controlled modulation of a nano-contact spin torque oscillator. Three principally different cases of frequency non-linearity ($d^{2}f/dI^{2}_{dc}$ being zero, positive, and negative) are investigated. Standard non-linear frequency modulation theory is able to accurately describe the frequency shifts during modulation. However, the power of the modulated sidebands only agrees with calculations based on a recent theory of combined non-linear frequency and amplitude modulation.

The intermittent contact impact problem in piping systems of nuclear reactor

International Nuclear Information System (INIS)

Martin, A.; Ricard, A.; Millard, A.

1981-09-01

The intermittent contact problem is important in many pipe whip studies, specially as to the safety of nuclear reactors. The impact concept adopted is that of instantaneous impact, so that at the time of impact the two impacting structures instantaneously acquire the same velocity in the impact direction. Energy is dissipated by some mechanism whose spatial and temporal scale is small compared to these scales in the discrete model. This dissipation is associated with local plastic deformation. Different solutions are presented for solving this problem. The first one is a generalization of the modal superposition method, when the nonlinearities of the structure are only due to impact between structural components; the other ones are included in a step by step time history and can take in account geometrical non linearities and of behavior. Some industrial applications in nuclear technology are presented

Quaternionic contact Einstein structures and the quaternionic contact Yamabe problem

CERN Document Server

Ivanov, Stefan; Vassilev, Dimiter

2014-01-01

A partial solution of the quaternionic contact Yamabe problem on the quaternionic sphere is given. It is shown that the torsion of the Biquard connection vanishes exactly when the trace-free part of the horizontal Ricci tensor of the Biquard connection is zero and this occurs precisely on 3-Sasakian manifolds. All conformal transformations sending the standard flat torsion-free quaternionic contact structure on the quaternionic Heisenberg group to a quaternionic contact structure with vanishing torsion of the Biquard connection are explicitly described. A "3-Hamiltonian form" of infinitesimal conformal automorphisms of quaternionic contact structures is presented.

Combined algorithms in nonlinear problems of magnetostatics

International Nuclear Information System (INIS)

Gregus, M.; Khoromskij, B.N.; Mazurkevich, G.E.; Zhidkov, E.P.

1988-01-01

To solve boundary problems of magnetostatics in unbounded two- and three-dimensional regions, we construct combined algorithms based on a combination of the method of boundary integral equations with the grid methods. We study the question of substantiation of the combined method of nonlinear magnetostatic problem without the preliminary discretization of equations and give some results on the convergence of iterative processes that arise in non-linear cases. We also discuss economical iterative processes and algorithms that solve boundary integral equations on certain surfaces. Finally, examples of numerical solutions of magnetostatic problems that arose when modelling the fields of electrophysical installations are given too. 14 refs.; 2 figs.; 1 tab

NONLINEAR DYNAMICS OF A ROTOR WITH CANTILEVERED DISK RESTING ON ANGULAR CONTACT BALL BEARINGS

Directory of Open Access Journals (Sweden)

S. Filipkovskyi

2016-06-01

Full Text Available The mathematical model of nonlinear oscillations of the rotor resting on angular contact ball bearings is developed. The disc is fixed on the console end of the shaft. The deflection of the shaft, and the elastic deformation of the bearings have the same order. Analysis of free oscillations is carried out, using nonlinear normal modes. The modes and backbone curves of rotor nonlinear oscillations are calculated. The system has soft characteristics.

Jacobian projection reduced-order models for dynamic systems with contact nonlinearities

Science.gov (United States)

Gastaldi, Chiara; Zucca, Stefano; Epureanu, Bogdan I.

2018-02-01

In structural dynamics, the prediction of the response of systems with localized nonlinearities, such as friction dampers, is of particular interest. This task becomes especially cumbersome when high-resolution finite element models are used. While state-of-the-art techniques such as Craig-Bampton component mode synthesis are employed to generate reduced order models, the interface (nonlinear) degrees of freedom must still be solved in-full. For this reason, a new generation of specialized techniques capable of reducing linear and nonlinear degrees of freedom alike is emerging. This paper proposes a new technique that exploits spatial correlations in the dynamics to compute a reduction basis. The basis is composed of a set of vectors obtained using the Jacobian of partial derivatives of the contact forces with respect to nodal displacements. These basis vectors correspond to specifically chosen boundary conditions at the contacts over one cycle of vibration. The technique is shown to be effective in the reduction of several models studied using multiple harmonics with a coupled static solution. In addition, this paper addresses another challenge common to all reduction techniques: it presents and validates a novel a posteriori error estimate capable of evaluating the quality of the reduced-order solution without involving a comparison with the full-order solution.

Higher-order techniques for some problems of nonlinear control

Directory of Open Access Journals (Sweden)

Sarychev Andrey V.

2002-01-01

Full Text Available A natural first step when dealing with a nonlinear problem is an application of some version of linearization principle. This includes the well known linearization principles for controllability, observability and stability and also first-order optimality conditions such as Lagrange multipliers rule or Pontryagin's maximum principle. In many interesting and important problems of nonlinear control the linearization principle fails to provide a solution. In the present paper we provide some examples of how higher-order methods of differential geometric control theory can be used for the study nonlinear control systems in such cases. The presentation includes: nonlinear systems with impulsive and distribution-like inputs; second-order optimality conditions for bang–bang extremals of optimal control problems; methods of high-order averaging for studying stability and stabilization of time-variant control systems.

Advanced Research Workshop on Nonlinear Hyperbolic Problems

CERN Document Server

Serre, Denis; Raviart, Pierre-Arnaud

1987-01-01

The field of nonlinear hyperbolic problems has been expanding very fast over the past few years, and has applications - actual and potential - in aerodynamics, multifluid flows, combustion, detonics amongst other. The difficulties that arise in application are of theoretical as well as numerical nature. In fact, the papers in this volume of proceedings deal to a greater extent with theoretical problems emerging in the resolution of nonlinear hyperbolic systems than with numerical methods. The volume provides an excellent up-to-date review of the current research trends in this area.

Nonlinear singular perturbation problems of arbitrary real orders

International Nuclear Information System (INIS)

Bijura, Angelina M.

2003-10-01

Higher order asymptotic solutions of singularly perturbed nonlinear fractional integral and derivatives of order 1/2 are investigated. It is particularly shown that whilst certain asymptotic expansions are applied successfully to linear equations and particular nonlinear problems, the standard formal asymptotic expansion is appropriate for the general class of nonlinear equations. This theory is then generalised to the general equation (of order β, 0 < β < 1). (author)

A simple numerical model of a geometrically nonlinear Timoshenko beam

NARCIS (Netherlands)

Keijdener, C.; Metrikine, A.

2015-01-01

In the original problem for which this model was developed, onedimensional flexible objects interact through a non-linear contact model. Due to the non-linear nature of the contact model, a numerical time-domain approach was adopted. One of the goals was to see if the coupling between axial and

Solving Large Scale Nonlinear Eigenvalue Problem in Next-Generation Accelerator Design

Energy Technology Data Exchange (ETDEWEB)

Liao, Ben-Shan; Bai, Zhaojun; /UC, Davis; Lee, Lie-Quan; Ko, Kwok; /SLAC

2006-09-28

A number of numerical methods, including inverse iteration, method of successive linear problem and nonlinear Arnoldi algorithm, are studied in this paper to solve a large scale nonlinear eigenvalue problem arising from finite element analysis of resonant frequencies and external Q{sub e} values of a waveguide loaded cavity in the next-generation accelerator design. They present a nonlinear Rayleigh-Ritz iterative projection algorithm, NRRIT in short and demonstrate that it is the most promising approach for a model scale cavity design. The NRRIT algorithm is an extension of the nonlinear Arnoldi algorithm due to Voss. Computational challenges of solving such a nonlinear eigenvalue problem for a full scale cavity design are outlined.

Analytical Solutions to Non-linear Mechanical Oscillation Problems

DEFF Research Database (Denmark)

Kaliji, H. D.; Ghadimi, M.; Barari, Amin

2011-01-01

In this paper, the Max-Min Method is utilized for solving the nonlinear oscillation problems. The proposed approach is applied to three systems with complex nonlinear terms in their motion equations. By means of this method, the dynamic behavior of oscillation systems can be easily approximated u...

Actuation of atomic force microscopy microcantilevers using contact acoustic nonlinearities

Energy Technology Data Exchange (ETDEWEB)

Torello, D.; Degertekin, F. Levent, E-mail: levent.degertekin@me.gatech.edu [George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332 (United States)

2013-11-15

A new method of actuating atomic force microscopy (AFM) cantilevers is proposed in which a high frequency (>5 MHz) wave modulated by a lower frequency (∼300 kHz) wave passes through a contact acoustic nonlinearity at the contact interface between the actuator and the cantilever chip. The nonlinearity converts the high frequency, modulated signal to a low frequency drive signal suitable for actuation of tapping-mode AFM probes. The higher harmonic content of this signal is filtered out mechanically by the cantilever transfer function, providing for clean output. A custom probe holder was designed and constructed using rapid prototyping technologies and off-the-shelf components and was interfaced with an Asylum Research MFP-3D AFM, which was then used to evaluate the performance characteristics with respect to standard hardware and linear actuation techniques. Using a carrier frequency of 14.19 MHz, it was observed that the cantilever output was cleaner with this actuation technique and added no significant noise to the system. This setup, without any optimization, was determined to have an actuation bandwidth on the order of 10 MHz, suitable for high speed imaging applications. Using this method, an image was taken that demonstrates the viability of the technique and is compared favorably to images taken with a standard AFM setup.

A solution to nonlinearity problems

International Nuclear Information System (INIS)

Neuffer, D.V.

1989-01-01

New methods of correcting dynamic nonlinearities resulting from the multipole content of a synchrotron or transport line are presented. In a simplest form, correction elements are places at the center (C) of the accelerator half-cells as well as near the focusing (F) and defocusing (D) quadrupoles. In a first approximation, the corrector strengths follow Simpson's Rule, forming an accurate quasi-local canceling approximation to the nonlinearity. The F, C, and D correctors may also be used to obtain precise control of the horizontal, coupled, and vertical motion. Correction by three or more orders of magnitude can be obtained, and simple solutions to a fundamental problem in beam transport have been obtained. 13 refs., 1 fig., 1 tab

A Linearized Relaxing Algorithm for the Specific Nonlinear Optimization Problem

Directory of Open Access Journals (Sweden)

Mio Horai

2016-01-01

Full Text Available We propose a new method for the specific nonlinear and nonconvex global optimization problem by using a linear relaxation technique. To simplify the specific nonlinear and nonconvex optimization problem, we transform the problem to the lower linear relaxation form, and we solve the linear relaxation optimization problem by the Branch and Bound Algorithm. Under some reasonable assumptions, the global convergence of the algorithm is certified for the problem. Numerical results show that this method is more efficient than the previous methods.

Cross-constrained problems for nonlinear Schrodinger equation with harmonic potential

Directory of Open Access Journals (Sweden)

Runzhang Xu

2012-11-01

Full Text Available This article studies a nonlinear Schodinger equation with harmonic potential by constructing different cross-constrained problems. By comparing the different cross-constrained problems, we derive different sharp criterion and different invariant manifolds that separate the global solutions and blowup solutions. Moreover, we conclude that some manifolds are empty due to the essence of the cross-constrained problems. Besides, we compare the three cross-constrained problems and the three depths of the potential wells. In this way, we explain the gaps in [J. Shu and J. Zhang, Nonlinear Shrodinger equation with harmonic potential, Journal of Mathematical Physics, 47, 063503 (2006], which was pointed out in [R. Xu and Y. Liu, Remarks on nonlinear Schrodinger equation with harmonic potential, Journal of Mathematical Physics, 49, 043512 (2008].

A two-dimensional finite element method for analysis of solid body contact problems in fuel rod mechanics

International Nuclear Information System (INIS)

Nissen, K.L.

1988-06-01

Two computer codes for the analysis of fuel rod behavior have been developed. Fuel rod mechanics is treated by a two-dimensional, axisymmetric finite element method. The program KONTAKT is used for detailed examinations on fuel rod sections, whereas the second program METHOD2D allows instationary calculations of whole fuel rods. The mechanical contact of fuel and cladding during heating of the fuel rod is very important for it's integrity. Both computer codes use a Newton-Raphson iteration for the solution of the nonlinear solid body contact problem. A constitutive equation is applied for the dependency of contact pressure on normal approach of the surfaces which are assumed to be rough. If friction is present on the contacting surfaces, Coulomb's friction law is used. Code validation is done by comparison with known analytical solutions for special problems. Results of the contact algorithm for an elastic ball pressing against a rigid surface are confronted with Hertzian theory. Influences of fuel-pellet geometry as well as influences of discretisation of displacements and stresses of a single fuel pellet are studied. Contact of fuel and cladding is calculated for a fuel rod section with two fuel pellets. The influence of friction forces between fuel and cladding on their axial expansion is demonstrated. By calculation of deformations and temperatures during an instationary fuel rod experiment of the CABRI-series the feasibility of two-dimensional finite element analysis of whole fuel rods is shown. (orig.) [de

Nonlinear dynamic model for skidding behavior of angular contact ball bearings

Science.gov (United States)

Han, Qinkai; Chu, Fulei

2015-10-01

A three-dimensional nonlinear dynamic model is proposed to predict the skidding behavior of angular contact ball bearings under combined load condition. The centrifugal and gyroscopic effects induced by ball rotation and revolution, Hertz contact between the ball and inner/outer races, discontinuous contact between the ball and cage and elastohydrodynamic lubrication are considered in the model. Through comparisons with the tested results of the reference, the dynamic model is verified. Based upon these, variations of ball slipping speed with time and space are discussed for the bearing under combined load condition. It is shown that radial load leads to the fluctuations in the slipping velocity of the ball contacting with inner/outer races, especially for the ball in load-decreasing regions. Adding the radial load would significantly increase the amplitude and range of slipping velocity, indicating that the skidding becomes more serious. As the ball still withstands contact load in the load-decreasing region, large slipping velocity would increase the temperature of both bearing and lubricant oil, intensify the wear and then might shorten the bearing service life. Therefore, the radial load should be considered carefully in the design and monitoring of rotating machinery.

Renormgroup symmetries in problems of nonlinear geometrical optics

International Nuclear Information System (INIS)

Kovalev, V.F.

1996-01-01

Utilization and further development of the previously announced approach [1,2] enables one to construct renormgroup symmetries for a boundary value problem for the system of equations which describes propagation of a powerful radiation in a nonlinear medium in geometrical optics approximation. With the help of renormgroup symmetries new rigorous and approximate analytical solutions of nonlinear geometrical optics equations are obtained. Explicit analytical expressions are presented that characterize spatial evolution of laser beam which has an arbitrary intensity dependence at the boundary of the nonlinear medium. (author)

Non-linear finite element analysis in structural mechanics

CERN Document Server

Rust, Wilhelm

2015-01-01

This monograph describes the numerical analysis of non-linearities in structural mechanics, i.e. large rotations, large strain (geometric non-linearities), non-linear material behaviour, in particular elasto-plasticity as well as time-dependent behaviour, and contact. Based on that, the book treats stability problems and limit-load analyses, as well as non-linear equations of a large number of variables. Moreover, the author presents a wide range of problem sets and their solutions. The target audience primarily comprises advanced undergraduate and graduate students of mechanical and civil engineering, but the book may also be beneficial for practising engineers in industry.

Renormalization-group approach to nonlinear radiation-transport problems

International Nuclear Information System (INIS)

Chapline, G.F.

1980-01-01

A Monte Carlo method is derived for solving nonlinear radiation-transport problems that allows one to average over the effects of many photon absorptions and emissions at frequencies where the opacity is large. This method should allow one to treat radiation-transport problems with large optical depths, e.g., line-transport problems, with little increase in computational effort over that which is required for optically thin problems

Numerical methods for solution of some nonlinear problems of mathematical physics

International Nuclear Information System (INIS)

Zhidkov, E.P.

1981-01-01

The continuous analog of the Newton method and its application to some nonlinear problems of mathematical physics using a computer is considered. It is shown that the application of this method in JINR to the wide range of nonlinear problems has shown its universality and high efficiency [ru

A frictional contact problem for an electro-viscoelastic body

Directory of Open Access Journals (Sweden)

Mircea Sofonea

2007-12-01

Full Text Available A mathematical model which describes the quasistatic frictional contact between a piezoelectric body and a deformable conductive foundation is studied. A nonlinear electro-viscoelastic constitutive law is used to model the piezoelectric material. Contact is described with the normal compliance condition, a version of Coulomb's law of dry friction, and a regularized electrical conductivity condition. A variational formulation of the model, in the form of a coupled system for the displacements and the electric potential, is derived. The existence of a unique weak solution of the model is established under a smallness assumption on the surface conductance. The proof is based on arguments of evolutionary variational inequalities and fixed points of operators.

A parametric FE modeling of brake for non-linear analysis

Energy Technology Data Exchange (ETDEWEB)

Ahmed,Ibrahim; Fatouh, Yasser [Automotive and Tractors Technology Department, Faculty of Industrial Education, Helwan University, Cairo (Egypt); Aly, Wael [Refrigeration and Air-Conditioning Technology Department, Faculty of Industrial Education, Helwan University, Cairo (Egypt)

2013-07-01

A parametric modeling of a drum brake based on 3-D Finite Element Methods (FEM) for non-contact analysis is presented. Many parameters are examined during this study such as the effect of drum-lining interface stiffness, coefficient of friction, and line pressure on the interface contact. Firstly, the modal analysis of the drum brake is also studied to get the natural frequency and instability of the drum to facilitate transforming the modal elements to non-contact elements. It is shown that the Unsymmetric solver of the modal analysis is efficient enough to solve this linear problem after transforming the non-linear behavior of the contact between the drum and the lining to a linear behavior. A SOLID45 which is a linear element is used in the modal analysis and then transferred to non-linear elements which are Targe170 and Conta173 that represent the drum and lining for contact analysis study. The contact analysis problems are highly non-linear and require significant computer resources to solve it, however, the contact problem give two significant difficulties. Firstly, the region of contact is not known based on the boundary conditions such as line pressure, and drum and friction material specs. Secondly, these contact problems need to take the friction into consideration. Finally, it showed a good distribution of the nodal reaction forces on the slotted lining contact surface and existing of the slot in the middle of the lining can help in wear removal due to the friction between the lining and the drum. Accurate contact stiffness can give a good representation for the pressure distribution between the lining and the drum. However, a full contact of the front part of the slotted lining could occur in case of 20, 40, 60 and 80 bar of piston pressure and a partially contact between the drum and lining can occur in the rear part of the slotted lining.

Morozov-type discrepancy principle for nonlinear ill-posed problems ...

Indian Academy of Sciences (India)

For proving the existence of a regularization parameter under a Morozov-type discrepancy principle for Tikhonov regularization of nonlinear ill-posed problems, it is required to impose additional nonlinearity assumptions on the forward operator. Lipschitz continuity of the Freéchet derivative and requirement of the Lipschitz ...

Morozov-type discrepancy principle for nonlinear ill-posed problems ...

Indian Academy of Sciences (India)

2016-08-26

Aug 26, 2016 ... For proving the existence of a regularization parameter under a Morozov-type discrepancy principle for Tikhonov regularization of nonlinear ill-posed problems, it is required to impose additional nonlinearity assumptions on the forward operator. Lipschitz continuity of the Freéchet derivative and requirement ...

Global Optimization of Nonlinear Blend-Scheduling Problems

Directory of Open Access Journals (Sweden)

Pedro A. Castillo Castillo

2017-04-01

Full Text Available The scheduling of gasoline-blending operations is an important problem in the oil refining industry. This problem not only exhibits the combinatorial nature that is intrinsic to scheduling problems, but also non-convex nonlinear behavior, due to the blending of various materials with different quality properties. In this work, a global optimization algorithm is proposed to solve a previously published continuous-time mixed-integer nonlinear scheduling model for gasoline blending. The model includes blend recipe optimization, the distribution problem, and several important operational features and constraints. The algorithm employs piecewise McCormick relaxation (PMCR and normalized multiparametric disaggregation technique (NMDT to compute estimates of the global optimum. These techniques partition the domain of one of the variables in a bilinear term and generate convex relaxations for each partition. By increasing the number of partitions and reducing the domain of the variables, the algorithm is able to refine the estimates of the global solution. The algorithm is compared to two commercial global solvers and two heuristic methods by solving four examples from the literature. Results show that the proposed global optimization algorithm performs on par with commercial solvers but is not as fast as heuristic approaches.

Bonus algorithm for large scale stochastic nonlinear programming problems

CERN Document Server

Diwekar, Urmila

2015-01-01

This book presents the details of the BONUS algorithm and its real world applications in areas like sensor placement in large scale drinking water networks, sensor placement in advanced power systems, water management in power systems, and capacity expansion of energy systems. A generalized method for stochastic nonlinear programming based on a sampling based approach for uncertainty analysis and statistical reweighting to obtain probability information is demonstrated in this book. Stochastic optimization problems are difficult to solve since they involve dealing with optimization and uncertainty loops. There are two fundamental approaches used to solve such problems. The first being the decomposition techniques and the second method identifies problem specific structures and transforms the problem into a deterministic nonlinear programming problem. These techniques have significant limitations on either the objective function type or the underlying distributions for the uncertain variables. Moreover, these ...

Initial boundary value problems of nonlinear wave equations in an exterior domain

International Nuclear Information System (INIS)

Chen Yunmei.

1987-06-01

In this paper, we investigate the existence and uniqueness of the global solutions to the initial boundary value problems of nonlinear wave equations in an exterior domain. When the space dimension n >= 3, the unique global solution of the above problem is obtained for small initial data, even if the nonlinear term is fully nonlinear and contains the unknown function itself. (author). 10 refs

Some New Results in Astrophysical Problems of Nonlinear Theory of Radiative Transfer

Science.gov (United States)

Pikichyan, H. V.

2017-07-01

In the interpretation of the observed astrophysical spectra, a decisive role is related to nonlinear problems of radiative transfer, because the processes of multiple interactions of matter of cosmic medium with the exciting intense radiation ubiquitously occur in astrophysical objects, and in their vicinities. Whereas, the intensity of the exciting radiation changes the physical properties of the original medium, and itself was modified, simultaneously, in a self-consistent manner under its influence. In the present report, we show that the consistent application of the principle of invariance in the nonlinear problem of bilateral external illumination of a scattering/absorbing one-dimensional anisotropic medium of finite geometrical thickness allows for simplifications that were previously considered as a prerogative only of linear problems. The nonlinear problem is analyzed through the three methods of the principle of invariance: (i) an adding of layers, (ii) its limiting form, described by differential equations of invariant imbedding, and (iii) a transition to the, so-called, functional equations of the "Ambartsumyan's complete invariance". Thereby, as an alternative to the Boltzmann equation, a new type of equations, so-called "kinetic equations of equivalence", are obtained. By the introduction of new functions - the so-called "linear images" of solution of nonlinear problem of radiative transfer, the linear structure of the solution of the nonlinear problem under study is further revealed. Linear images allow to convert naturally the statistical characteristics of random walk of a "single quantum" or their "beam of unit intensity", as well as widely known "probabilistic interpretation of phenomena of transfer", to the field of nonlinear problems. The structure of the equations obtained for determination of linear images is typical of linear problems.

SEACAS Theory Manuals: Part 1. Problem Formulation in Nonlinear Solid Mechancis

Energy Technology Data Exchange (ETDEWEB)

Attaway, S.W.; Laursen, T.A.; Zadoks, R.I.

1998-08-01

This report gives an introduction to the basic concepts and principles involved in the formulation of nonlinear problems in solid mechanics. By way of motivation, the discussion begins with a survey of some of the important sources of nonlinearity in solid mechanics applications, using wherever possible simple one dimensional idealizations to demonstrate the physical concepts. This discussion is then generalized by presenting generic statements of initial/boundary value problems in solid mechanics, using linear elasticity as a template and encompassing such ideas as strong and weak forms of boundary value problems, boundary and initial conditions, and dynamic and quasistatic idealizations. The notational framework used for the linearized problem is then extended to account for finite deformation of possibly inelastic solids, providing the context for the descriptions of nonlinear continuum mechanics, constitutive modeling, and finite element technology given in three companion reports.

Linear differential equations to solve nonlinear mechanical problems: A novel approach

OpenAIRE

Nair, C. Radhakrishnan

2004-01-01

Often a non-linear mechanical problem is formulated as a non-linear differential equation. A new method is introduced to find out new solutions of non-linear differential equations if one of the solutions of a given non-linear differential equation is known. Using the known solution of the non-linear differential equation, linear differential equations are set up. The solutions of these linear differential equations are found using standard techniques. Then the solutions of the linear differe...

Implicit solvers for large-scale nonlinear problems

International Nuclear Information System (INIS)

Keyes, David E; Reynolds, Daniel R; Woodward, Carol S

2006-01-01

Computational scientists are grappling with increasingly complex, multi-rate applications that couple such physical phenomena as fluid dynamics, electromagnetics, radiation transport, chemical and nuclear reactions, and wave and material propagation in inhomogeneous media. Parallel computers with large storage capacities are paving the way for high-resolution simulations of coupled problems; however, hardware improvements alone will not prove enough to enable simulations based on brute-force algorithmic approaches. To accurately capture nonlinear couplings between dynamically relevant phenomena, often while stepping over rapid adjustments to quasi-equilibria, simulation scientists are increasingly turning to implicit formulations that require a discrete nonlinear system to be solved for each time step or steady state solution. Recent advances in iterative methods have made fully implicit formulations a viable option for solution of these large-scale problems. In this paper, we overview one of the most effective iterative methods, Newton-Krylov, for nonlinear systems and point to software packages with its implementation. We illustrate the method with an example from magnetically confined plasma fusion and briefly survey other areas in which implicit methods have bestowed important advantages, such as allowing high-order temporal integration and providing a pathway to sensitivity analyses and optimization. Lastly, we overview algorithm extensions under development motivated by current SciDAC applications

Response of Non-Linear Shock Absorbers-Boundary Value Problem Analysis

Science.gov (United States)

Rahman, M. A.; Ahmed, U.; Uddin, M. S.

2013-08-01

A nonlinear boundary value problem of two degrees-of-freedom (DOF) untuned vibration damper systems using nonlinear springs and dampers has been numerically studied. As far as untuned damper is concerned, sixteen different combinations of linear and nonlinear springs and dampers have been comprehensively analyzed taking into account transient terms. For different cases, a comparative study is made for response versus time for different spring and damper types at three important frequency ratios: one at r = 1, one at r > 1 and one at r <1. The response of the system is changed because of the spring and damper nonlinearities; the change is different for different cases. Accordingly, an initially stable absorber may become unstable with time and vice versa. The analysis also shows that higher nonlinearity terms make the system more unstable. Numerical simulation includes transient vibrations. Although problems are much more complicated compared to those for a tuned absorber, a comparison of the results generated by the present numerical scheme with the exact one shows quite a reasonable agreement

Non-linear analytic and coanalytic problems (Lp-theory, Clifford analysis, examples)

International Nuclear Information System (INIS)

Dubinskii, Yu A; Osipenko, A S

2000-01-01

Two kinds of new mathematical model of variational type are put forward: non-linear analytic and coanalytic problems. The formulation of these non-linear boundary-value problems is based on a decomposition of the complete scale of Sobolev spaces into the 'orthogonal' sum of analytic and coanalytic subspaces. A similar decomposition is considered in the framework of Clifford analysis. Explicit examples are presented

Bending of a nonlinear beam reposing on an unilateral foundation

Directory of Open Access Journals (Sweden)

Machalová J.

2011-06-01

Full Text Available This article is going to deal with bending of a nonlinear beam whose mathematical model was proposed by D. Y. Gao in (Gao, D. Y., Nonlinear elastic beam theory with application in contact problems and variational approaches,Mech. Research Communication, 23 (1 1996. The model is based on the Euler-Bernoulli hypothesis and under assumption of nonzero lateral stress component enables moderately large deflections but with small strains. This is here extended by the unilateralWinkler foundation. The attribution unilateral means that the foundation is not connected with the beam. For this problem we demonstrate a mathematical formulation resulting from its natural decomposition which leads to a saddle-point problem with a proper Lagrangian. Next we are concerned with methods of solution for our problem by means of the finite element method as the paper (Gao, D. Y., Nonlinear elastic beam theory with application in contact problems and variational approaches, Mech. Research Communication, 23 (1 1996 has no mention of it. The main alternatives are here the solution of a system of nonlinear nondifferentiable equations or finding of a saddle point through the use of the augmented Lagrangian method. This is illustrated by an example in the final part of the article.

Morozov-type discrepancy principle for nonlinear ill-posed problems ...

Indian Academy of Sciences (India)

[3] Engl H W, Kunisch K and Neubauer A, Convergence rates for Tikhonov regularization of nonliner problems, Inverse Problems 5 (1989) 523–540. [4] Hanke M, Neubauer A and Scherzer O, A convergence analysis of Landweber iteration for nonlinear ill-posed problems, Numer. Math. 72 (1995) 21–37. [5] Hofmann B and ...

Nonlinear problems in theoretical physics

International Nuclear Information System (INIS)

Ranada, A.F.

1979-01-01

This volume contains the lecture notes and review talks delivered at the 9th GIFT international seminar on theoretical physics on the general subject 'Nonlinear Problems in Theoretical Physics'. Mist contributions deal with recent developments in the theory of the spectral transformation and solitons, but there are also articles from the field of transport theory and plasma physics and an unconventional view of classical and quantum electrodynamics. All contributions to this volume will appear under their corresponding subject categories. (HJ)

Influence of Stress Shape Function on Analysis of Contact Problem Using Hybrid Photoelasticity

International Nuclear Information System (INIS)

Shin, Dongchul; Hawong, Jaisug

2013-01-01

In this research, a study on stress shape functions was conducted to analyze the contact stress problem by using a hybrid photoelasticity. Because the contact stress problem is generally solved as a half-plane problem, the relationship between two analytical stress functions, which are compositions of the Airy stress function, was similar to one of the crack problem. However, this relationship in itself could not be used to solve the contact stress problem (especially one with singular points). Therefore, to analyze the contact stress problem more correctly, stress shape functions based on the condition of two contact end points had to be considered in the form of these two analytical stress functions. The four types of stress shape functions were related to the stress singularities at the two contact end points. Among them, the primary two types used for the analysis of an O-ring were selected, and their validities were verified in this work

On a non-linear pseudodifferential boundary value problem

International Nuclear Information System (INIS)

Nguyen Minh Chuong.

1989-12-01

A pseudodifferential boundary value problem for operators with symbols taking values in Sobolev spaces and with non-linear right-hand side was studied. Existence and uniqueness theorems were proved. (author). 11 refs

Simultaneous multigrid techniques for nonlinear eigenvalue problems: Solutions of the nonlinear Schrödinger-Poisson eigenvalue problem in two and three dimensions

Science.gov (United States)

Costiner, Sorin; Ta'asan, Shlomo

1995-07-01

Algorithms for nonlinear eigenvalue problems (EP's) often require solving self-consistently a large number of EP's. Convergence difficulties may occur if the solution is not sought in an appropriate region, if global constraints have to be satisfied, or if close or equal eigenvalues are present. Multigrid (MG) algorithms for nonlinear problems and for EP's obtained from discretizations of partial differential EP have often been shown to be more efficient than single level algorithms. This paper presents MG techniques and a MG algorithm for nonlinear Schrödinger Poisson EP's. The algorithm overcomes the above mentioned difficulties combining the following techniques: a MG simultaneous treatment of the eigenvectors and nonlinearity, and with the global constrains; MG stable subspace continuation techniques for the treatment of nonlinearity; and a MG projection coupled with backrotations for separation of solutions. These techniques keep the solutions in an appropriate region, where the algorithm converges fast, and reduce the large number of self-consistent iterations to only a few or one MG simultaneous iteration. The MG projection makes it possible to efficiently overcome difficulties related to clusters of close and equal eigenvalues. Computational examples for the nonlinear Schrödinger-Poisson EP in two and three dimensions, presenting special computational difficulties that are due to the nonlinearity and to the equal and closely clustered eigenvalues are demonstrated. For these cases, the algorithm requires O(qN) operations for the calculation of q eigenvectors of size N and for the corresponding eigenvalues. One MG simultaneous cycle per fine level was performed. The total computational cost is equivalent to only a few Gauss-Seidel relaxations per eigenvector. An asymptotic convergence rate of 0.15 per MG cycle is attained.

Simulation of unilateral contact problems departing from the classical boundary problems

International Nuclear Information System (INIS)

Frey, S.L.; Sampaio, R.; Gama, R.M.S. da.

1989-08-01

A numerical algorithm is proposed for simulating unilateral contact problems under the classical elasticity point of view. This simple algorithm may be employed by engineers with a minimum knowledge on classical elasticity. (A.C.A.S.) [pt

Nonlinear problems in fluid dynamics and inverse scattering: Nonlinear waves and inverse scattering

Science.gov (United States)

Ablowitz, Mark J.

1994-12-01

Research investigations involving the fundamental understanding and applications of nonlinear wave motion and related studies of inverse scattering and numerical computation have been carried out and a number of significant results have been obtained. A class of nonlinear wave equations which can be solved by the inverse scattering transform (IST) have been studied, including the Kadaomtsev-Petviashvili (KP) equation, the Davey-Stewartson equation, and the 2+1 Toda system. The solutions obtained by IST correspond to the Cauchy initial value problem with decaying initial data. We have also solved two important systems via the IST method: a 'Volterra' system in 2+1 dimensions and a new one dimensional nonlinear equation which we refer to as the Toda differential-delay equation. Research in computational chaos in moderate to long time numerical simulations continues.

Adomian decomposition method for nonlinear Sturm-Liouville problems

Directory of Open Access Journals (Sweden)

Sennur Somali

2007-09-01

Full Text Available In this paper the Adomian decomposition method is applied to the nonlinear Sturm-Liouville problem-y" + y(tp=λy(t, y(t > 0, t ∈ I = (0, 1, y(0 = y(1 = 0, where p > 1 is a constant and λ > 0 is an eigenvalue parameter. Also, the eigenvalues and the behavior of eigenfuctions of the problem are demonstrated.

Selected Problems in Nonlinear Dynamics and Sociophysics

Science.gov (United States)

Westley, Alexandra Renee

This Ph.D. dissertation focuses on a collection of problems on the dynamical behavior of nonlinear many-body systems, drawn from two substantially different areas. First, the dynamical behavior seen in strongly nonlinear lattices such as in the Fermi-Pasta-Ulam-Tsingou (FPUT) system (part I) and second, time evolution behavior of interacting living objects which can be broadly considered as sociophysics systems (part II). The studies on FPUT-like systems will comprise of five chapters, dedicated to the properties of solitary and anti-solitary waves in the system, how localized nonlinear excitations decay and spread throughout these lattices, how two colliding solitary waves can precipitate highly localized and stable excitations, a possible alternative way to view these localized excitations through Duffing oscillators, and finally an exploration of parametric resonance in an FPUT-like lattice. Part II consists of two problems in the context of sociophysics. I use molecular dynamics inspired simulations to study the size and the stability of social groups of chimpanzees (such as those seen in central Africa) and compare the results with existing observations on the stability of chimpanzee societies. Secondly, I use an agent-based model to simulate land battles between an intelligent army and an insurgency when both have access to equally powerful weaponry. The study considers genetic algorithm based adaptive strategies to infer the strategies needed for the intelligent army to win the battles.

Nonlinear finite element analyses: advances and challenges in dental applications.

Science.gov (United States)

Wakabayashi, N; Ona, M; Suzuki, T; Igarashi, Y

2008-07-01

To discuss the development and current status of application of nonlinear finite element method (FEM) in dentistry. The literature was searched for original research articles with keywords such as nonlinear, finite element analysis, and tooth/dental/implant. References were selected manually or searched from the PUBMED and MEDLINE databases through November 2007. The nonlinear problems analyzed in FEM studies were reviewed and categorized into: (A) nonlinear simulations of the periodontal ligament (PDL), (B) plastic and viscoelastic behaviors of dental materials, (C) contact phenomena in tooth-to-tooth contact, (D) contact phenomena within prosthodontic structures, and (E) interfacial mechanics between the tooth and the restoration. The FEM in dentistry recently focused on simulation of realistic intra-oral conditions such as the nonlinear stress-strain relationship in the periodontal tissues and the contact phenomena in teeth, which could hardly be solved by the linear static model. The definition of contact area critically affects the reliability of the contact analyses, especially for implant-abutment complexes. To predict the failure risk of a bonded tooth-restoration interface, it is essential to assess the normal and shear stresses relative to the interface. The inclusion of viscoelasticity and plastic deformation to the program to account for the time-dependent, thermal sensitive, and largely deformable nature of dental materials would enhance its application. Further improvement of the nonlinear FEM solutions should be encouraged to widen the range of applications in dental and oral health science.

Analytic approximations to nonlinear boundary value problems modeling beam-type nano-electromechanical systems

Energy Technology Data Exchange (ETDEWEB)

Zou, Li [Dalian Univ. of Technology, Dalian City (China). State Key Lab. of Structural Analysis for Industrial Equipment; Liang, Songxin; Li, Yawei [Dalian Univ. of Technology, Dalian City (China). School of Mathematical Sciences; Jeffrey, David J. [Univ. of Western Ontario, London (Canada). Dept. of Applied Mathematics

2017-06-01

Nonlinear boundary value problems arise frequently in physical and mechanical sciences. An effective analytic approach with two parameters is first proposed for solving nonlinear boundary value problems. It is demonstrated that solutions given by the two-parameter method are more accurate than solutions given by the Adomian decomposition method (ADM). It is further demonstrated that solutions given by the ADM can also be recovered from the solutions given by the two-parameter method. The effectiveness of this method is demonstrated by solving some nonlinear boundary value problems modeling beam-type nano-electromechanical systems.

A new integrability theory for certain nonlinear physical problems

International Nuclear Information System (INIS)

Berger, M.S.

1993-01-01

A new mathematically sound integrability theory for certain nonlinear problems defined by ordinary or partial differential equations is defined. The new theory works in an arbitrary finite number of space dimensions. Moreover, if a system is integrable in the new sense described here, it has a remarkable stability property that distinguishes if from any previously known integrability ideas. The new theory proceeds by establishing a ''global normal form'' for the problem at hand. This normal form holds subject to canonical coordinate transformations, extending such classical ideas by using new nonlinear methods of infinite dimensional functional analysis. The global normal form in question is related to the mathematical theory of singularities of mappings of H. Whitney and R. Thom extended globally and form finite to infinite dimensions. Thus bifurcation phenomena are naturally included in the new integrability theory. Typical examples include the classically nonintegrable Riccati equation, certain non-Euclidean mean field theories, certain parabolic reaction diffusion equations and the hyperbolic nonlinear telegrapher's equation. (Author)

Mathematical Modeling of Contact Problems of Elasticity Theory with Unilateral Discrete Contact

Directory of Open Access Journals (Sweden)

I. V. Stankevich

2015-01-01

Full Text Available Development and operation of modern machinery and latest technology require reliable estimates of the strength characteristics of the critical elements of structures and technological equipment under the impact of high-intensity thermomechanical loading, accompanied, as a rule, by complex contact interaction. Mathematical modeling of stress-strain state of such parts and components in the contact area, based on adequate mathematical models, modern numerical methods and efficient algorithms that implement the direct determination of displacement fields, strains and stresses, is the main tool that allows fast acquisition of data required for the calculations of strength and durability. The paper considers an algorithm for constructing the numerical solution of the contact problem of elasticity theory in relation to the body, which has an obvious one-sided discrete contact interaction with an elastic half-space. The proposed algorithm is specially designed to have a correction of the tangential forces at discrete contact points, allowing us to achieve sufficiently accurate implementation of the adopted law of friction. The algorithm is embedded in a general finite element technology, with which the application code is generated. Numerical study of discrete unilateral contact interaction of an elastic plate and a rigid half-space showed a high efficiency of the developed algorithm and the application code that implements it.

Computational Methods for Nonlinear Dynamics Problems in Solid and Structural Mechanics: Models of Dynamic Frictional Phenomena in Metallic Structures.

Science.gov (United States)

1986-03-31

Martins, J.A.C. and Campos , L.T. [1986], "Existence and Local Uniqueness of Solutions to Contact Problems in Elasticity with Nonlinear Friction...noisy and ttoubl esome vibt.t4ons. If the sound generated by the friction-induced oscillations of Rviolin strings may be the delight of all music lovers...formulation. See 0den and Martins - [1985] and Rabier, Martins, Oden and Campos [1986]. - It is now simple to show, in a 6o’uman manner, that, for

Non-steady state modelling of wheel-rail contact problem

Science.gov (United States)

Guiral, A.; Alonso, A.; Baeza, L.; Giménez, J. G.

2013-01-01

Among all the algorithms to solve the wheel-rail contact problem, Kalker's FastSim has become the most useful computation tool since it combines a low computational cost and enough precision for most of the typical railway dynamics problems. However, some types of dynamic problems require the use of a non-steady state analysis. Alonso and Giménez developed a non-stationary method based on FastSim, which provides both, sufficiently accurate results and a low computational cost. However, it presents some limitations; the method is developed for one time-dependent creepage and its accuracy for varying normal forces has not been checked. This article presents the required changes in order to deal with both problems and compares its results with those given by Kalker's Variational Method for rolling contact.

High efficiency algorithm for 3D transient thermo-elasto-plastic contact problem in reactor pressure vessel sealing system

International Nuclear Information System (INIS)

Xu Mingyu; Lin Tengjiao; Li Runfang; Du Xuesong; Li Shuian; Yang Yu

2005-01-01

There are some complex operating cases such as high temperature and high pressure during the operating process of nuclear reactor pressure vessel. It is necessary to carry out mechanical analysis and experimental investigation for its sealing ability. On the basis of the self-developed program for 3-D transient sealing analysis for nuclear reactor pressure vessel, some specific measures are presented to enhance the calculation efficiency in several aspects such as the non-linear solution of elasto-plastic problem, the mixed solution algorithm for contact problem as well as contract heat transfer problem and linear equation set solver. The 3-D transient sealing analysis program is amended and complemented, with which the sealing analysis result of the pressure vessel model can be obtained. The calculation results have good regularity and the calculation efficiency is twice more than before. (authors)

Some nonlinear problems in the manipulation of beams

International Nuclear Information System (INIS)

Sessler, A.M.

1990-01-01

An overview is given of nonlinear problems that arise in the manipulation of beams. Beams can be made of material particles or photons, can be intense or dilute, can be energetic or not, and they can be propagating in vacuum or in a medium. The nonlinear aspects of the motion are different in each case, and this diversity of behavior is categorized. Many examples are given, which serves to illustrate the categorization and, furthermore, display the richness of behavior encountered in the physics of beams. 25 refs., 5 figs

A gradient stable scheme for a phase field model for the moving contact line problem

KAUST Repository

Gao, Min

2012-02-01

In this paper, an efficient numerical scheme is designed for a phase field model for the moving contact line problem, which consists of a coupled system of the Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition [1,2,4]. The nonlinear version of the scheme is semi-implicit in time and is based on a convex splitting of the Cahn-Hilliard free energy (including the boundary energy) together with a projection method for the Navier-Stokes equations. We show, under certain conditions, the scheme has the total energy decaying property and is unconditionally stable. The linearized scheme is easy to implement and introduces only mild CFL time constraint. Numerical tests are carried out to verify the accuracy and stability of the scheme. The behavior of the solution near the contact line is examined. It is verified that, when the interface intersects with the boundary, the consistent splitting scheme [21,22] for the Navier Stokes equations has the better accuracy for pressure. © 2011 Elsevier Inc.

Non-linear analytic and coanalytic problems ( L_p-theory, Clifford analysis, examples)

Science.gov (United States)

Dubinskii, Yu A.; Osipenko, A. S.

2000-02-01

Two kinds of new mathematical model of variational type are put forward: non-linear analytic and coanalytic problems. The formulation of these non-linear boundary-value problems is based on a decomposition of the complete scale of Sobolev spaces into the "orthogonal" sum of analytic and coanalytic subspaces. A similar decomposition is considered in the framework of Clifford analysis. Explicit examples are presented.

The Cauchy problem for non-linear Klein-Gordon equations

International Nuclear Information System (INIS)

Simon, J.C.H.; Taflin, E.

1993-01-01

We consider in R n+1 , n≥2, the non-linear Klein-Gordon equation. We prove for such an equation that there is neighbourhood of zero in a Hilbert space of initial conditions for which the Cauchy problem has global solutions and on which there is asymptotic completeness. The inverse of the wave operator linearizes the non-linear equation. If, moreover, the equation is manifestly Poincare covariant then the non-linear representation of the Poincare-Lie algebra, associated with the non-linear Klein-Gordon equation is integrated to a non-linear representation of the Poincare group on an invariant neighbourhood of zero in the Hilbert space. This representation is linearized by the inverse of the wave operator. The Hilbert space is, in both cases, the closure of the space of the differentiable vectors for the linear representation of the Poincare group, associated with the Klein-Gordon equation, with respect to a norm defined by the representation of the enveloping algebra. (orig.)

Prediction of bead area contact load at the tire-wheel interface using NASTRAN

Science.gov (United States)

Chen, C. H. S.

1982-01-01

The theoretical prediction of the bead area contact load at the tire wheel interface using NASTRAN is reported. The application of the linear code to a basically nonlinear problem results in excessive deformation of the structure and the tire-wheel contact conditions become impossible to achieve. A psuedo-nonlinear approach was adopted in which the moduli of the cord reinforced composite are increased so that the computed key deformations matched that of the experiment. Numerical results presented are discussed.

Intrinsic nonlinearity and method of disturbed observations in inverse problems of celestial mechanics

Science.gov (United States)

Avdyushev, Victor A.

2017-12-01

Orbit determination from a small sample of observations over a very short observed orbital arc is a strongly nonlinear inverse problem. In such problems an evaluation of orbital uncertainty due to random observation errors is greatly complicated, since linear estimations conventionally used are no longer acceptable for describing the uncertainty even as a rough approximation. Nevertheless, if an inverse problem is weakly intrinsically nonlinear, then one can resort to the so-called method of disturbed observations (aka observational Monte Carlo). Previously, we showed that the weaker the intrinsic nonlinearity, the more efficient the method, i.e. the more accurate it enables one to simulate stochastically the orbital uncertainty, while it is strictly exact only when the problem is intrinsically linear. However, as we ascertained experimentally, its efficiency was found to be higher than that of other stochastic methods widely applied in practice. In the present paper we investigate the intrinsic nonlinearity in complicated inverse problems of Celestial Mechanics when orbits are determined from little informative samples of observations, which typically occurs for recently discovered asteroids. To inquire into the question, we introduce an index of intrinsic nonlinearity. In asteroid problems it evinces that the intrinsic nonlinearity can be strong enough to affect appreciably probabilistic estimates, especially at the very short observed orbital arcs that the asteroids travel on for about a hundredth of their orbital periods and less. As it is known from regression analysis, the source of intrinsic nonlinearity is the nonflatness of the estimation subspace specified by a dynamical model in the observation space. Our numerical results indicate that when determining asteroid orbits it is actually very slight. However, in the parametric space the effect of intrinsic nonlinearity is exaggerated mainly by the ill-conditioning of the inverse problem. Even so, as for the

On a Highly Nonlinear Self-Obstacle Optimal Control Problem

Energy Technology Data Exchange (ETDEWEB)

Di Donato, Daniela, E-mail: daniela.didonato@unitn.it [University of Trento, Department of Mathematics (Italy); Mugnai, Dimitri, E-mail: dimitri.mugnai@unipg.it [Università di Perugia, Dipartimento di Matematica e Informatica (Italy)

2015-10-15

We consider a non-quadratic optimal control problem associated to a nonlinear elliptic variational inequality, where the obstacle is the control itself. We show that, fixed a desired profile, there exists an optimal solution which is not far from it. Detailed characterizations of the optimal solution are given, also in terms of approximating problems.

Experimental analysis of nonlinear problems in solid mechanics

International Nuclear Information System (INIS)

1982-01-01

The booklet presents abstracts of papers from the Euromech Colloqium No. 152 held from Sept. 20th to 24th, 1982 in Wuppertal, Federal Republic of Germany. All the papers are dealing with Experimental Analysis of Nonlinear Problems in Solid Mechanics. (RW)

Lectures on nonlinear evolution equations initial value problems

CERN Document Server

Racke, Reinhard

2015-01-01

This book mainly serves as an elementary, self-contained introduction to several important aspects of the theory of global solutions to initial value problems for nonlinear evolution equations. The book employs the classical method of continuation of local solutions with the help of a priori estimates obtained for small data. The existence and uniqueness of small, smooth solutions that are defined for all values of the time parameter are investigated. Moreover, the asymptotic behavior of the solutions is described as time tends to infinity. The methods for nonlinear wave equations are discussed in detail. Other examples include the equations of elasticity, heat equations, the equations of thermoelasticity, Schrödinger equations, Klein-Gordon equations, Maxwell equations and plate equations. To emphasize the importance of studying the conditions under which small data problems offer global solutions, some blow-up results are briefly described. Moreover, the prospects for corresponding initial-boundary value p...

New Exact Penalty Functions for Nonlinear Constrained Optimization Problems

Directory of Open Access Journals (Sweden)

Bingzhuang Liu

2014-01-01

Full Text Available For two kinds of nonlinear constrained optimization problems, we propose two simple penalty functions, respectively, by augmenting the dimension of the primal problem with a variable that controls the weight of the penalty terms. Both of the penalty functions enjoy improved smoothness. Under mild conditions, it can be proved that our penalty functions are both exact in the sense that local minimizers of the associated penalty problem are precisely the local minimizers of the original constrained problem.

Galerkin approximations of nonlinear optimal control problems in Hilbert spaces

Directory of Open Access Journals (Sweden)

Mickael D. Chekroun

2017-07-01

Full Text Available Nonlinear optimal control problems in Hilbert spaces are considered for which we derive approximation theorems for Galerkin approximations. Approximation theorems are available in the literature. The originality of our approach relies on the identification of a set of natural assumptions that allows us to deal with a broad class of nonlinear evolution equations and cost functionals for which we derive convergence of the value functions associated with the optimal control problem of the Galerkin approximations. This convergence result holds for a broad class of nonlinear control strategies as well. In particular, we show that the framework applies to the optimal control of semilinear heat equations posed on a general compact manifold without boundary. The framework is then shown to apply to geoengineering and mitigation of greenhouse gas emissions formulated here in terms of optimal control of energy balance climate models posed on the sphere $\\mathbb{S}^2$.

Analytical vs. Simulation Solution Techniques for Pulse Problems in Non-linear Stochastic Dynamics

DEFF Research Database (Denmark)

Iwankiewicz, R.; Nielsen, Søren R. K.

Advantages and disadvantages of available analytical and simulation techniques for pulse problems in non-linear stochastic dynamics are discussed. First, random pulse problems, both those which do and do not lead to Markov theory, are presented. Next, the analytical and analytically-numerical tec......Advantages and disadvantages of available analytical and simulation techniques for pulse problems in non-linear stochastic dynamics are discussed. First, random pulse problems, both those which do and do not lead to Markov theory, are presented. Next, the analytical and analytically...

A quadratic approximation-based algorithm for the solution of multiparametric mixed-integer nonlinear programming problems

KAUST Repository

Domí nguez, Luis F.; Pistikopoulos, Efstratios N.

2012-01-01

An algorithm for the solution of convex multiparametric mixed-integer nonlinear programming problems arising in process engineering problems under uncertainty is introduced. The proposed algorithm iterates between a multiparametric nonlinear

Automatic contact algorithm in ppercase[dyna3d] for crashworthiness and impact problems

International Nuclear Information System (INIS)

Whirley, Robert G.; Engelmann, Bruce E.

1994-01-01

This paper presents a new approach for the automatic definition and treatment of mechanical contact in explicit non-linear finite element analysis. Automatic contact offers the benefits of significantly reduced model construction time and fewer opportunities for user error, but faces significant challenges in reliability and computational costs. Key aspects of the proposed new method include automatic identification of adjacent and opposite surfaces in the global search phase, and the use of a well-defined surface normal which allows a consistent treatment of shell intersection and corner contact conditions without adhoc rules. The paper concludes with three examples which illustrate the performance of the newly proposed algorithm in the public ppercase[dyna3d] code. ((orig.))

On the solvability of initial boundary value problems for nonlinear ...

African Journals Online (AJOL)

In this paper, we study the initial boundary value problems for a non-linear time dependent Schrödinger equation with Dirichlet and Neumann boundary conditions, respectively. We prove the existence and uniqueness of solutions of the initial boundary value problems by using Galerkin's method. Keywords: Initial boundary ...

First international conference on nonlinear problems in aviation and aerospace

International Nuclear Information System (INIS)

Sivasundaram, S.

1994-01-01

The International Conference on Nonlinear Problems in Aviation and Aerospace was held at Embry-Riddle Aeronautical University, Daytona Beach, Florida on May 9-11, 1996. This conference was sponsored by the International Federation of Nonlinear Analysts, International Federation of Information Processing, and Embry-Riddle Aeronautical University. Over one hundred engineers, scientists, and mathematicians from seventeen countries attended. These proceedings include keynote addresses, invited lectures, and contributed papers presented during the conference

Some problems on nonlinear hyperbolic equations and applications

CERN Document Server

Peng, YueJun

2010-01-01

This volume is composed of two parts: Mathematical and Numerical Analysis for Strongly Nonlinear Plasma Models and Exact Controllability and Observability for Quasilinear Hyperbolic Systems and Applications. It presents recent progress and results obtained in the domains related to both subjects without attaching much importance to the details of proofs but rather to difficulties encountered, to open problems and possible ways to be exploited. It will be very useful for promoting further study on some important problems in the future.

From a Nonlinear, Nonconvex Variational Problem to a Linear, Convex Formulation

International Nuclear Information System (INIS)

Egozcue, J.; Meziat, R.; Pedregal, P.

2002-01-01

We propose a general approach to deal with nonlinear, nonconvex variational problems based on a reformulation of the problem resulting in an optimization problem with linear cost functional and convex constraints. As a first step we explicitly explore these ideas to some one-dimensional variational problems and obtain specific conclusions of an analytical and numerical nature

On the real-time estimation of the wheel-rail contact force by means of a new nonlinear estimator design model

Science.gov (United States)

Strano, Salvatore; Terzo, Mario

2018-05-01

The dynamics of the railway vehicles is strongly influenced by the interaction between the wheel and the rail. This kind of contact is affected by several conditioning factors such as vehicle speed, wear, adhesion level and, moreover, it is nonlinear. As a consequence, the modelling and the observation of this kind of phenomenon are complex tasks but, at the same time, they constitute a fundamental step for the estimation of the adhesion level or for the vehicle condition monitoring. This paper presents a novel technique for the real time estimation of the wheel-rail contact forces based on an estimator design model that takes into account the nonlinearities of the interaction by means of a fitting model functional to reproduce the contact mechanics in a wide range of slip and to be easily integrated in a complete model based estimator for railway vehicle.

Nonlinear reaction-diffusion equations with delay: some theorems, test problems, exact and numerical solutions

Science.gov (United States)

Polyanin, A. D.; Sorokin, V. G.

2017-12-01

The paper deals with nonlinear reaction-diffusion equations with one or several delays. We formulate theorems that allow constructing exact solutions for some classes of these equations, which depend on several arbitrary functions. Examples of application of these theorems for obtaining new exact solutions in elementary functions are provided. We state basic principles of construction, selection, and use of test problems for nonlinear partial differential equations with delay. Some test problems which can be suitable for estimating accuracy of approximate analytical and numerical methods of solving reaction-diffusion equations with delay are presented. Some examples of numerical solutions of nonlinear test problems with delay are considered.

On Newton-Raphson formulation and algorithm for displacement based structural dynamics problem with quadratic damping nonlinearity

Directory of Open Access Journals (Sweden)

Koh Kim Jie

2017-01-01

Full Text Available Quadratic damping nonlinearity is challenging for displacement based structural dynamics problem as the problem is nonlinear in time derivative of the primitive variable. For such nonlinearity, the formulation of tangent stiffness matrix is not lucid in the literature. Consequently, ambiguity related to kinematics update arises when implementing the time integration-iterative algorithm. In present work, an Euler-Bernoulli beam vibration problem with quadratic damping nonlinearity is addressed as the main source of quadratic damping nonlinearity arises from drag force estimation, which is generally valid only for slender structures. Employing Newton-Raphson formulation, tangent stiffness components associated with quadratic damping nonlinearity requires velocity input for evaluation purpose. For this reason, two mathematically equivalent algorithm structures with different kinematics arrangement are tested. Both algorithm structures result in the same accuracy and convergence characteristic of solution.

Arbitrary Lagrangian-Eulerian method for non-linear problems of geomechanics

International Nuclear Information System (INIS)

Nazem, M; Carter, J P; Airey, D W

2010-01-01

In many geotechnical problems it is vital to consider the geometrical non-linearity caused by large deformation in order to capture a more realistic model of the true behaviour. The solutions so obtained should then be more accurate and reliable, which should ultimately lead to cheaper and safer design. The Arbitrary Lagrangian-Eulerian (ALE) method originated from fluid mechanics, but has now been well established for solving large deformation problems in geomechanics. This paper provides an overview of the ALE method and its challenges in tackling problems involving non-linearities due to material behaviour, large deformation, changing boundary conditions and time-dependency, including material rate effects and inertia effects in dynamic loading applications. Important aspects of ALE implementation into a finite element framework will also be discussed. This method is then employed to solve some interesting and challenging geotechnical problems such as the dynamic bearing capacity of footings on soft soils, consolidation of a soil layer under a footing, and the modelling of dynamic penetration of objects into soil layers.

Consensus problem in directed networks of multi-agents via nonlinear protocols

International Nuclear Information System (INIS)

Liu Xiwei; Chen Tianping; Lu Wenlian

2009-01-01

In this Letter, the consensus problem via distributed nonlinear protocols for directed networks is investigated. Its dynamical behaviors are described by ordinary differential equations (ODEs). Based on graph theory, matrix theory and the Lyapunov direct method, some sufficient conditions of nonlinear protocols guaranteeing asymptotical or exponential consensus are presented and rigorously proved. The main contribution of this work is that for nonlinearly coupled networks, we generalize the results for undirected networks to directed networks. Consensus under pinning control technique is also developed here. Simulations are also given to show the validity of the theories.

Inverse Boundary Value Problem for Non-linear Hyperbolic Partial Differential Equations

OpenAIRE

Nakamura, Gen; Vashisth, Manmohan

2017-01-01

In this article we are concerned with an inverse boundary value problem for a non-linear wave equation of divergence form with space dimension $n\\geq 3$. This non-linear wave equation has a trivial solution, i.e. zero solution. By linearizing this equation at the trivial solution, we have the usual linear isotropic wave equation with the speed $\\sqrt{\\gamma(x)}$ at each point $x$ in a given spacial domain. For any small solution $u=u(t,x)$ of this non-linear equation, we have the linear isotr...

Asymptotics for inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation

Energy Technology Data Exchange (ETDEWEB)

Kaikina, Elena I., E-mail: ekaikina@matmor.unam.mx [Centro de Ciencias Matemáticas, UNAM Campus Morelia, AP 61-3 (Xangari), Morelia CP 58089, Michoacán (Mexico)

2013-11-15

We consider the inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation, formulated on a half-line. We study traditionally important problems of the theory of nonlinear partial differential equations, such as global in time existence of solutions to the initial-boundary value problem and the asymptotic behavior of solutions for large time.

Asymptotics for inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation

International Nuclear Information System (INIS)

Kaikina, Elena I.

2013-01-01

We consider the inhomogeneous Dirichlet initial-boundary value problem for the nonlinear Schrödinger equation, formulated on a half-line. We study traditionally important problems of the theory of nonlinear partial differential equations, such as global in time existence of solutions to the initial-boundary value problem and the asymptotic behavior of solutions for large time

History-Dependent Problems with Applications to Contact Models for Elastic Beams

International Nuclear Information System (INIS)

Bartosz, Krzysztof; Kalita, Piotr; Migórski, Stanisław; Ochal, Anna; Sofonea, Mircea

2016-01-01

We prove an existence and uniqueness result for a class of subdifferential inclusions which involve a history-dependent operator. Then we specialize this result in the study of a class of history-dependent hemivariational inequalities. Problems of such kind arise in a large number of mathematical models which describe quasistatic processes of contact. To provide an example we consider an elastic beam in contact with a reactive obstacle. The contact is modeled with a new and nonstandard condition which involves both the subdifferential of a nonconvex and nonsmooth function and a Volterra-type integral term. We derive a variational formulation of the problem which is in the form of a history-dependent hemivariational inequality for the displacement field. Then, we use our abstract result to prove its unique weak solvability. Finally, we consider a numerical approximation of the model, solve effectively the approximate problems and provide numerical simulations

History-Dependent Problems with Applications to Contact Models for Elastic Beams

Energy Technology Data Exchange (ETDEWEB)

Bartosz, Krzysztof; Kalita, Piotr; Migórski, Stanisław; Ochal, Anna, E-mail: ochal@ii.uj.edu.pl [Jagiellonian University, Faculty of Mathematics and Computer Science (Poland); Sofonea, Mircea [Université de Perpignan Via Domitia, Laboratoire de Mathématiques et Physique (France)

2016-02-15

We prove an existence and uniqueness result for a class of subdifferential inclusions which involve a history-dependent operator. Then we specialize this result in the study of a class of history-dependent hemivariational inequalities. Problems of such kind arise in a large number of mathematical models which describe quasistatic processes of contact. To provide an example we consider an elastic beam in contact with a reactive obstacle. The contact is modeled with a new and nonstandard condition which involves both the subdifferential of a nonconvex and nonsmooth function and a Volterra-type integral term. We derive a variational formulation of the problem which is in the form of a history-dependent hemivariational inequality for the displacement field. Then, we use our abstract result to prove its unique weak solvability. Finally, we consider a numerical approximation of the model, solve effectively the approximate problems and provide numerical simulations.

Application of HPEM to investigate the response and stability of nonlinear problems in vibration

DEFF Research Database (Denmark)

Mohammadi, M.H.; Mohammadi, A.; Kimiaeifar, A.

2010-01-01

In this work, a powerful analytical method, called He's Parameter Expanding Methods (HPEM) is used to obtain the exact solution of nonlinear problems in nonlinear vibration. In this work, the governing equation is obtained by using Lagrange method, then the nonlinear governing equation is solved...

On discrete maximum principles for nonlinear elliptic problems

Czech Academy of Sciences Publication Activity Database

Karátson, J.; Korotov, S.; Křížek, Michal

2007-01-01

Roč. 76, č. 1 (2007), s. 99-108 ISSN 0378-4754 R&D Projects: GA MŠk 1P05ME749; GA AV ČR IAA1019201 Institutional research plan: CEZ:AV0Z10190503 Keywords : nonlinear elliptic problem * mixed boundary conditions * finite element method Subject RIV: BA - General Mathematics Impact factor: 0.738, year: 2007

Assessment of Two Analytical Methods in Solving the Linear and Nonlinear Elastic Beam Deformation Problems

DEFF Research Database (Denmark)

Barari, Amin; Ganjavi, B.; Jeloudar, M. Ghanbari

2010-01-01

and fluid mechanics. Design/methodology/approach – Two new but powerful analytical methods, namely, He's VIM and HPM, are introduced to solve some boundary value problems in structural engineering and fluid mechanics. Findings – Analytical solutions often fit under classical perturbation methods. However......, as with other analytical techniques, certain limitations restrict the wide application of perturbation methods, most important of which is the dependence of these methods on the existence of a small parameter in the equation. Disappointingly, the majority of nonlinear problems have no small parameter at all......Purpose – In the last two decades with the rapid development of nonlinear science, there has appeared ever-increasing interest of scientists and engineers in the analytical techniques for nonlinear problems. This paper considers linear and nonlinear systems that are not only regarded as general...

On some nonlinear problems arising in the physics of ionized gases

International Nuclear Information System (INIS)

Hilhorst-Goldman, D.

1981-01-01

The author reports results obtained by rigorous analysis of a nonlinear differential equation for the electron density nsub(e) in a specific type of electrical discharge. The problem is essentially two-dimensional. She discusses in particular the escape of electrons to infinity above a critical temperature and the boundary layer exhibited by nsub(e) near zero temperature. A singular boundary value problem arising in a pre-breakdown gas discharge is discussed. A Coulomb gas is considered in a special experimental situation: the pre-breakdown gas discharge between two electrodes. The equation for the negative charge density can be formulated as a nonlinear parabolic equation degenerate at the origin. The existence and uniqueness of the solution are proved as well as the asymptotic stability of its unique steady state. Some results are also given about the rate of convergence. The variational characterisation of the limit solution of a singular perturbation problem and variational analysis of a perturbed free boundary problem are considered. (Auth./C.F.)

Non-steady state modeling of wheel-rail contact problem

OpenAIRE

Guiral, A.; Alonso, A.; Baeza González, Luis Miguel; Giménez, J.G.

2013-01-01

Among all the algorithms to solve the wheel–rail contact problem, Kalker's FastSim has become the most useful computation tool since it combines a low computational cost and enough precision for most of the typical railway dynamics problems. However, some types of dynamic problems require the use of a non-steady state analysis. Alonso and Giménez developed a non-stationary method based on FastSim, which provides both, sufficiently accurate results and a low computational cost. However, it pre...

Inverse periodic problem for the discrete approximation of the Schroedinger nonlinear equation

International Nuclear Information System (INIS)

Bogolyubov, N.N.; Prikarpatskij, A.K.; AN Ukrainskoj SSR, Lvov. Inst. Prikladnykh Problem Mekhaniki i Matematiki)

1982-01-01

The problem of numerical solution of the Schroedinger nonlinear equation (1) iPSIsub(t) = PSIsub(xx)+-2(PSI)sup(2)PSI. The numerical solution of nonlinear differential equation supposes its discrete approximation is required for the realization of the computer calculation process. Tor the equation (1) there exists the following discrete approximation by variable x(2) iPSIsub(n, t) = (PSIsub(n+1)-2PSIsub(n)+PSIsub(n-1))/(Δx)sup(2)+-(PSIsub(n))sup(2)(PSIsub(n+1)+PSIsub(n-1)), n=0, +-1, +-2... where PSIsub(n)(+) is the corresponding value of PSI(x, t) function in the node and divisions with the equilibrium step Δx. The main problem is obtaining analytically exact solutions of the equations (2). The analysis of the equation system (2) is performed on the base of the discrete analogue of the periodic variant of the inverse scattering problem method developed with the aid of nonlinear equations of the Korteweg-de Vries type. Obtained in explicit form are analytical solutions of the equations system (2). The solutions are expressed through the Riemann THETA-function [ru

Multigrid techniques for nonlinear eigenvalue probems: Solutions of a nonlinear Schroedinger eigenvalue problem in 2D and 3D

Science.gov (United States)

Costiner, Sorin; Taasan, Shlomo

1994-01-01

This paper presents multigrid (MG) techniques for nonlinear eigenvalue problems (EP) and emphasizes an MG algorithm for a nonlinear Schrodinger EP. The algorithm overcomes the mentioned difficulties combining the following techniques: an MG projection coupled with backrotations for separation of solutions and treatment of difficulties related to clusters of close and equal eigenvalues; MG subspace continuation techniques for treatment of the nonlinearity; an MG simultaneous treatment of the eigenvectors at the same time with the nonlinearity and with the global constraints. The simultaneous MG techniques reduce the large number of self consistent iterations to only a few or one MG simultaneous iteration and keep the solutions in a right neighborhood where the algorithm converges fast.

Modeling of fatigue crack induced nonlinear ultrasonics using a highly parallelized explicit local interaction simulation approach

Science.gov (United States)

Shen, Yanfeng; Cesnik, Carlos E. S.

2016-04-01

This paper presents a parallelized modeling technique for the efficient simulation of nonlinear ultrasonics introduced by the wave interaction with fatigue cracks. The elastodynamic wave equations with contact effects are formulated using an explicit Local Interaction Simulation Approach (LISA). The LISA formulation is extended to capture the contact-impact phenomena during the wave damage interaction based on the penalty method. A Coulomb friction model is integrated into the computation procedure to capture the stick-slip contact shear motion. The LISA procedure is coded using the Compute Unified Device Architecture (CUDA), which enables the highly parallelized supercomputing on powerful graphic cards. Both the explicit contact formulation and the parallel feature facilitates LISA's superb computational efficiency over the conventional finite element method (FEM). The theoretical formulations based on the penalty method is introduced and a guideline for the proper choice of the contact stiffness is given. The convergence behavior of the solution under various contact stiffness values is examined. A numerical benchmark problem is used to investigate the new LISA formulation and results are compared with a conventional contact finite element solution. Various nonlinear ultrasonic phenomena are successfully captured using this contact LISA formulation, including the generation of nonlinear higher harmonic responses. Nonlinear mode conversion of guided waves at fatigue cracks is also studied.

A new solution method for wheel/rail rolling contact.

Science.gov (United States)

Yang, Jian; Song, Hua; Fu, Lihua; Wang, Meng; Li, Wei

2016-01-01

To solve the problem of wheel/rail rolling contact of nonlinear steady-state curving, a three-dimensional transient finite element (FE) model is developed by the explicit software ANSYS/LS-DYNA. To improve the solving speed and efficiency, an explicit-explicit order solution method is put forward based on analysis of the features of implicit and explicit algorithm. The solution method was first applied to calculate the pre-loading of wheel/rail rolling contact with explicit algorithm, and then the results became the initial conditions in solving the dynamic process of wheel/rail rolling contact with explicit algorithm as well. Simultaneously, the common implicit-explicit order solution method is used to solve the FE model. Results show that the explicit-explicit order solution method has faster operation speed and higher efficiency than the implicit-explicit order solution method while the solution accuracy is almost the same. Hence, the explicit-explicit order solution method is more suitable for the wheel/rail rolling contact model with large scale and high nonlinearity.

Techniques for contact and contact with friction problems; Tecnicas para problemas de contacto y contacto con friccion

Energy Technology Data Exchange (ETDEWEB)

Velandia Arana, Gonzalo [Instituto de Investigaciones Electricas, Cuernavaca (Mexico)

1990-12-31

Different numerical techniques are presented based in the finite element method to obtain numerical solutions to contact and contact with friction problems between solid bodies, and compared between each other. [Espanol] Se presentan diferentes tecnicas numericas basadas en el metodo de elementos finitos para la obtencion de soluciones numericas de problemas de contacto y contacto con friccion entre cuerpos solidos, y se comparan entre si.

Techniques for contact and contact with friction problems; Tecnicas para problemas de contacto y contacto con friccion

Energy Technology Data Exchange (ETDEWEB)

Velandia Arana, Gonzalo [Instituto de Investigaciones Electricas, Cuernavaca (Mexico)

1989-12-31

Different numerical techniques are presented based in the finite element method to obtain numerical solutions to contact and contact with friction problems between solid bodies, and compared between each other. [Espanol] Se presentan diferentes tecnicas numericas basadas en el metodo de elementos finitos para la obtencion de soluciones numericas de problemas de contacto y contacto con friccion entre cuerpos solidos, y se comparan entre si.

Application of nonlinear Krylov acceleration to radiative transfer problems

International Nuclear Information System (INIS)

Till, A. T.; Adams, M. L.; Morel, J. E.

2013-01-01

The iterative solution technique used for radiative transfer is normally nested, with outer thermal iterations and inner transport iterations. We implement a nonlinear Krylov acceleration (NKA) method in the PDT code for radiative transfer problems that breaks nesting, resulting in more thermal iterations but significantly fewer total inner transport iterations. Using the metric of total inner transport iterations, we investigate a crooked-pipe-like problem and a pseudo-shock-tube problem. Using only sweep preconditioning, we compare NKA against a typical inner / outer method employing GMRES / Newton and find NKA to be comparable or superior. Finally, we demonstrate the efficacy of applying diffusion-based preconditioning to grey problems in conjunction with NKA. (authors)

Analytical Solution of Nonlinear Problems in Classical Dynamics by Means of Lagrange-Ham

DEFF Research Database (Denmark)

Kimiaeifar, Amin; Mahdavi, S. H; Rabbani, A.

2011-01-01

In this work, a powerful analytical method, called Homotopy Analysis Methods (HAM) is coupled with Lagrange method to obtain the exact solution for nonlinear problems in classic dynamics. In this work, the governing equations are obtained by using Lagrange method, and then the nonlinear governing...

The solution of a coupled system of nonlinear physical problems using the homotopy analysis method

International Nuclear Information System (INIS)

El-Wakil, S A; Abdou, M A

2010-01-01

In this article, the homotopy analysis method (HAM) has been applied to solve coupled nonlinear evolution equations in physics. The validity of this method has been successfully demonstrated by applying it to two nonlinear evolution equations, namely coupled nonlinear diffusion reaction equations and the (2+1)-dimensional Nizhnik-Novikov Veselov system. The results obtained by this method show good agreement with the ones obtained by other methods. The proposed method is a powerful and easy to use analytic tool for nonlinear problems and does not need small parameters in the equations. The HAM solutions contain an auxiliary parameter that provides a convenient way of controlling the convergence region of series solutions. The results obtained here reveal that the proposed method is very effective and simple for solving nonlinear evolution equations. The basic ideas of this approach can be widely employed to solve other strongly nonlinear problems.

Numerical algorithms for contact problems in linear elastostatics

International Nuclear Information System (INIS)

Barbosa, H.J.C.; Feijoo, R.A.

1984-01-01

In this work contact problems in linear elasticity are analysed by means of Finite Elements and Mathematical Programming Techniques. The principle of virtual work leads in this case to a variational inequality which in turn is equivalent, for Hookean materials and infinitesimal strains, to the minimization of the total potential energy over the set of all admissible virtual displacements. The use of Gauss-Seidel algorithm with relaxation and projection and also Lemke's algorithm and Uzawa's algorithm for solving the minimization problem is discussed. Finally numerical examples are presented. (Author) [pt

COYOTE: a finite element computer program for nonlinear heat conduction problems

International Nuclear Information System (INIS)

Gartling, D.K.

1978-06-01

COYOTE is a finite element computer program designed for the solution of two-dimensional, nonlinear heat conduction problems. The theoretical and mathematical basis used to develop the code is described. Program capabilities and complete user instructions are presented. Several example problems are described in detail to demonstrate the use of the program

On topological derivatives for contact problems in elasticity

Czech Academy of Sciences Publication Activity Database

Giusti, S.M.; Sokolowski, S.; Stebel, Jan

2015-01-01

Roč. 165, č. 1 (2015), s. 279-294 ISSN 0022-3239 R&D Projects: GA ČR GA201/09/0917 Institutional support: RVO:67985840 Keywords : topological derivative * static frictionless contact problem * asymptotic analysis Subject RIV: BA - General Mathematics Impact factor: 1.160, year: 2015 http://link.springer.com/article/10.1007%2Fs10957-014-0594-7

Nonlinear triple-point problems on time scales

Directory of Open Access Journals (Sweden)

Douglas R. Anderson

2004-04-01

Full Text Available We establish the existence of multiple positive solutions to the nonlinear second-order triple-point boundary-value problem on time scales, $$displaylines{ u^{Delta abla}(t+h(tf(t,u(t=0, cr u(a=alpha u(b+delta u^Delta(a,quad eta u(c+gamma u^Delta(c=0 }$$ for $tin[a,c]subsetmathbb{T}$, where $mathbb{T}$ is a time scale, $eta, gamma, deltage 0$ with $Beta+gamma>0$, $0

Efficient numerical methods for the large-scale, parallel solution of elastoplastic contact problems

KAUST Repository

Frohne, Jö rg; Heister, Timo; Bangerth, Wolfgang

2015-01-01

© 2016 John Wiley & Sons, Ltd. Quasi-static elastoplastic contact problems are ubiquitous in many industrial processes and other contexts, and their numerical simulation is consequently of great interest in accurately describing and optimizing production processes. The key component in these simulations is the solution of a single load step of a time iteration. From a mathematical perspective, the problems to be solved in each time step are characterized by the difficulties of variational inequalities for both the plastic behavior and the contact problem. Computationally, they also often lead to very large problems. In this paper, we present and evaluate a complete set of methods that are (1) designed to work well together and (2) allow for the efficient solution of such problems. In particular, we use adaptive finite element meshes with linear and quadratic elements, a Newton linearization of the plasticity, active set methods for the contact problem, and multigrid-preconditioned linear solvers. Through a sequence of numerical experiments, we show the performance of these methods. This includes highly accurate solutions of a three-dimensional benchmark problem and scaling our methods in parallel to 1024 cores and more than a billion unknowns.

Efficient numerical methods for the large-scale, parallel solution of elastoplastic contact problems

KAUST Repository

Frohne, Jörg

2015-08-06

© 2016 John Wiley & Sons, Ltd. Quasi-static elastoplastic contact problems are ubiquitous in many industrial processes and other contexts, and their numerical simulation is consequently of great interest in accurately describing and optimizing production processes. The key component in these simulations is the solution of a single load step of a time iteration. From a mathematical perspective, the problems to be solved in each time step are characterized by the difficulties of variational inequalities for both the plastic behavior and the contact problem. Computationally, they also often lead to very large problems. In this paper, we present and evaluate a complete set of methods that are (1) designed to work well together and (2) allow for the efficient solution of such problems. In particular, we use adaptive finite element meshes with linear and quadratic elements, a Newton linearization of the plasticity, active set methods for the contact problem, and multigrid-preconditioned linear solvers. Through a sequence of numerical experiments, we show the performance of these methods. This includes highly accurate solutions of a three-dimensional benchmark problem and scaling our methods in parallel to 1024 cores and more than a billion unknowns.

TENSOLVE: A software package for solving systems of nonlinear equations and nonlinear least squares problems using tensor methods

Energy Technology Data Exchange (ETDEWEB)

Bouaricha, A. [Argonne National Lab., IL (United States). Mathematics and Computer Science Div.; Schnabel, R.B. [Colorado Univ., Boulder, CO (United States). Dept. of Computer Science

1996-12-31

This paper describes a modular software package for solving systems of nonlinear equations and nonlinear least squares problems, using a new class of methods called tensor methods. It is intended for small to medium-sized problems, say with up to 100 equations and unknowns, in cases where it is reasonable to calculate the Jacobian matrix or approximate it by finite differences at each iteration. The software allows the user to select between a tensor method and a standard method based upon a linear model. The tensor method models F({ital x}) by a quadratic model, where the second-order term is chosen so that the model is hardly more expensive to form, store, or solve than the standard linear model. Moreover, the software provides two different global strategies, a line search and a two- dimensional trust region approach. Test results indicate that, in general, tensor methods are significantly more efficient and robust than standard methods on small and medium-sized problems in iterations and function evaluations.

Multiple solutions for inhomogeneous nonlinear elliptic problems arising in astrophyiscs

Directory of Open Access Journals (Sweden)

Marco Calahorrano

2004-04-01

Full Text Available Using variational methods we prove the existence and multiplicity of solutions for some nonlinear inhomogeneous elliptic problems on a bounded domain in $mathbb{R}^n$, with $ngeq 2$ and a smooth boundary, and when the domain is $mathbb{R}_+^n$

Lavrentiev regularization method for nonlinear ill-posed problems

International Nuclear Information System (INIS)

Kinh, Nguyen Van

2002-10-01

In this paper we shall be concerned with Lavientiev regularization method to reconstruct solutions x 0 of non ill-posed problems F(x)=y o , where instead of y 0 noisy data y δ is an element of X with absolut(y δ -y 0 ) ≤ δ are given and F:X→X is an accretive nonlinear operator from a real reflexive Banach space X into itself. In this regularization method solutions x α δ are obtained by solving the singularly perturbed nonlinear operator equation F(x)+α(x-x*)=y δ with some initial guess x*. Assuming certain conditions concerning the operator F and the smoothness of the element x*-x 0 we derive stability estimates which show that the accuracy of the regularized solutions is order optimal provided that the regularization parameter α has been chosen properly. (author)

Iterative solution of a nonlinear system arising in phase change problems

International Nuclear Information System (INIS)

Williams, M.A.

1987-01-01

We consider several iterative methods for solving the nonlinear system arising from an enthalpy formulation of a phase change problem. We present the formulation of the problem. Implicit discretization of the governing equations results in a mildly nonlinear system at each time step. We discuss solving this system using Jacobi, Gauss-Seidel, and SOR iterations and a new modified preconditioned conjugate gradient (MPCG) algorithm. The new MPCG algorithm and its properties are discussed in detail. Numerical results are presented comparing the performance of the SOR algorithm and the MPCG algorithm with 1-step SSOR preconditioning. The MPCG algorithm exhibits a superlinear rate of convergence. The SOR algorithm exhibits a linear rate of convergence. Thus, the MPCG algorithm requires fewer iterations to converge than the SOR algorithm. However in most cases, the SOR algorithm requires less total computation time than the MPCG algorithm. Hence, the SOR algorithm appears to be more appropriate for the class of problems considered. 27 refs., 11 figs

Systems of general nonlinear set-valued mixed variational inequalities problems in Hilbert spaces

Directory of Open Access Journals (Sweden)

Cho Yeol

2011-01-01

Full Text Available Abstract In this paper, the existing theorems and methods for finding solutions of systems of general nonlinear set-valued mixed variational inequalities problems in Hilbert spaces are studied. To overcome the difficulties, due to the presence of a proper convex lower semi-continuous function, φ and a mapping g, which appeared in the considered problem, we have used some applications of the resolvent operator technique. We would like to point out that although many authors have proved results for finding solutions of the systems of nonlinear set-valued (mixed variational inequalities problems, it is clear that it cannot be directly applied to the problems that we have considered in this paper because of φ and g. 2000 AMS Subject Classification: 47H05; 47H09; 47J25; 65J15.

Periodic Contact and Crack Problems in Plane Elasticity

DEFF Research Database (Denmark)

Krenk, Steen

1976-01-01

By use of singular integral equations it is demonstrated how some periodic contact and crack problems can be solved in closed form. The integral equation in question is the same which is encountered when dealing with mixed boundary conditions on a circle. As analytical evaluation of the solution ...... may be quite complicated attention is drawn to a numerical quadrature method...

New preconditioned conjugate gradient algorithms for nonlinear unconstrained optimization problems

International Nuclear Information System (INIS)

Al-Bayati, A.; Al-Asadi, N.

1997-01-01

This paper presents two new predilection conjugate gradient algorithms for nonlinear unconstrained optimization problems and examines their computational performance. Computational experience shows that the new proposed algorithms generally imp lone the efficiency of Nazareth's [13] preconditioned conjugate gradient algorithm. (authors). 16 refs., 1 tab

Research on the Problem of Spur Gear Teeth Contact in the Car Gear Box

Directory of Open Access Journals (Sweden)

Viktor Skrickij

2011-04-01

Full Text Available The article presents research on the problem of two gear contact in the car gearbox. Contact stiffness is evaluated for the whole period of mesh. Also, contact stresses are evaluated in the contact place. The presented method can be used for calculating spur gear.Article in Lithuanian

Modal analysis of graphene-based structures for large deformations, contact and material nonlinearities

Science.gov (United States)

Ghaffari, Reza; Sauer, Roger A.

2018-06-01

The nonlinear frequencies of pre-stressed graphene-based structures, such as flat graphene sheets and carbon nanotubes, are calculated. These structures are modeled with a nonlinear hyperelastic shell model. The model is calibrated with quantum mechanics data and is valid for high strains. Analytical solutions of the natural frequencies of various plates are obtained for the Canham bending model by assuming infinitesimal strains. These solutions are used for the verification of the numerical results. The performance of the model is illustrated by means of several examples. Modal analysis is performed for square plates under pure dilatation or uniaxial stretch, circular plates under pure dilatation or under the effects of an adhesive substrate, and carbon nanotubes under uniaxial compression or stretch. The adhesive substrate is modeled with van der Waals interaction (based on the Lennard-Jones potential) and a coarse grained contact model. It is shown that the analytical natural frequencies underestimate the real ones, and this should be considered in the design of devices based on graphene structures.

Bayesian nonlinear regression for large small problems

KAUST Repository

Chakraborty, Sounak; Ghosh, Malay; Mallick, Bani K.

2012-01-01

Statistical modeling and inference problems with sample sizes substantially smaller than the number of available covariates are challenging. This is known as large p small n problem. Furthermore, the problem is more complicated when we have multiple correlated responses. We develop multivariate nonlinear regression models in this setup for accurate prediction. In this paper, we introduce a full Bayesian support vector regression model with Vapnik's ε-insensitive loss function, based on reproducing kernel Hilbert spaces (RKHS) under the multivariate correlated response setup. This provides a full probabilistic description of support vector machine (SVM) rather than an algorithm for fitting purposes. We have also introduced a multivariate version of the relevance vector machine (RVM). Instead of the original treatment of the RVM relying on the use of type II maximum likelihood estimates of the hyper-parameters, we put a prior on the hyper-parameters and use Markov chain Monte Carlo technique for computation. We have also proposed an empirical Bayes method for our RVM and SVM. Our methods are illustrated with a prediction problem in the near-infrared (NIR) spectroscopy. A simulation study is also undertaken to check the prediction accuracy of our models. © 2012 Elsevier Inc.

Bayesian nonlinear regression for large small problems

KAUST Repository

Chakraborty, Sounak

2012-07-01

Statistical modeling and inference problems with sample sizes substantially smaller than the number of available covariates are challenging. This is known as large p small n problem. Furthermore, the problem is more complicated when we have multiple correlated responses. We develop multivariate nonlinear regression models in this setup for accurate prediction. In this paper, we introduce a full Bayesian support vector regression model with Vapnik\\'s ε-insensitive loss function, based on reproducing kernel Hilbert spaces (RKHS) under the multivariate correlated response setup. This provides a full probabilistic description of support vector machine (SVM) rather than an algorithm for fitting purposes. We have also introduced a multivariate version of the relevance vector machine (RVM). Instead of the original treatment of the RVM relying on the use of type II maximum likelihood estimates of the hyper-parameters, we put a prior on the hyper-parameters and use Markov chain Monte Carlo technique for computation. We have also proposed an empirical Bayes method for our RVM and SVM. Our methods are illustrated with a prediction problem in the near-infrared (NIR) spectroscopy. A simulation study is also undertaken to check the prediction accuracy of our models. © 2012 Elsevier Inc.

A Smooth Newton Method for Nonlinear Programming Problems with Inequality Constraints

Directory of Open Access Journals (Sweden)

Vasile Moraru

2012-02-01

Full Text Available The paper presents a reformulation of the Karush-Kuhn-Tucker (KKT system associated nonlinear programming problem into an equivalent system of smooth equations. Classical Newton method is applied to solve the system of equations. The superlinear convergence of the primal sequence, generated by proposed method, is proved. The preliminary numerical results with a problems test set are presented.

Testing and Modeling of Contact Problems in Resistance Welding

DEFF Research Database (Denmark)

Song, Quanfeng

together two or three cylindrical parts as well as disc-ring pairs of dissimilar metals. The tests have demonstrated the effectiveness of the model. A theoretical and experimental study is performed on the contact resistance aiming at a more reliable model for numerical simulation of resistance welding......As a part of the efforts towards a professional and reliable numerical tool for resistance welding engineers, this Ph.D. project is dedicated to refining the numerical models related to the interface behavior. An FE algorithm for the contact problems in resistance welding has been developed...... for the formulation, and the interfaces are treated in a symmetric pattern. The frictional sliding contact is also solved employing the constant friction model. The algorithm is incorporated into the finite element code. Verification is carried out in some numerical tests as well as experiments such as upsetting...

Nonlinear programming for classification problems in machine learning

Science.gov (United States)

Astorino, Annabella; Fuduli, Antonio; Gaudioso, Manlio

2016-10-01

We survey some nonlinear models for classification problems arising in machine learning. In the last years this field has become more and more relevant due to a lot of practical applications, such as text and web classification, object recognition in machine vision, gene expression profile analysis, DNA and protein analysis, medical diagnosis, customer profiling etc. Classification deals with separation of sets by means of appropriate separation surfaces, which is generally obtained by solving a numerical optimization model. While linear separability is the basis of the most popular approach to classification, the Support Vector Machine (SVM), in the recent years using nonlinear separating surfaces has received some attention. The objective of this work is to recall some of such proposals, mainly in terms of the numerical optimization models. In particular we tackle the polyhedral, ellipsoidal, spherical and conical separation approaches and, for some of them, we also consider the semisupervised versions.

On a mixed problem for a coupled nonlinear system

Directory of Open Access Journals (Sweden)

Marcondes R. Clark

1997-03-01

Full Text Available In this article we prove the existence and uniqueness of solutions to the mixed problem associated with the nonlinear system $$ u_{tt}-M(int_Omega |abla u|^2dxDelta u+|u|^ ho u+heta =f $$ $$ heta _t -Delta heta +u_{t}=g $$ where $M$ is a positive real function, and $f$ and $g$ are known real functions.

Green's function-stochastic methods framework for probing nonlinear evolution problems: Burger's equation, the nonlinear Schroedinger's equation, and hydrodynamic organization of near-molecular-scale vorticity

International Nuclear Information System (INIS)

Keanini, R.G.

2011-01-01

Research highlights: → Systematic approach for physically probing nonlinear and random evolution problems. → Evolution of vortex sheets corresponds to evolution of an Ornstein-Uhlenbeck process. → Organization of near-molecular scale vorticity mediated by hydrodynamic modes. → Framework allows calculation of vorticity evolution within random strain fields. - Abstract: A framework which combines Green's function (GF) methods and techniques from the theory of stochastic processes is proposed for tackling nonlinear evolution problems. The framework, established by a series of easy-to-derive equivalences between Green's function and stochastic representative solutions of linear drift-diffusion problems, provides a flexible structure within which nonlinear evolution problems can be analyzed and physically probed. As a preliminary test bed, two canonical, nonlinear evolution problems - Burgers' equation and the nonlinear Schroedinger's equation - are first treated. In the first case, the framework provides a rigorous, probabilistic derivation of the well known Cole-Hopf ansatz. Likewise, in the second, the machinery allows systematic recovery of a known soliton solution. The framework is then applied to a fairly extensive exploration of physical features underlying evolution of randomly stretched and advected Burger's vortex sheets. Here, the governing vorticity equation corresponds to the Fokker-Planck equation of an Ornstein-Uhlenbeck process, a correspondence that motivates an investigation of sub-sheet vorticity evolution and organization. Under the assumption that weak hydrodynamic fluctuations organize disordered, near-molecular-scale, sub-sheet vorticity, it is shown that these modes consist of two weakly damped counter-propagating cross-sheet acoustic modes, a diffusive cross-sheet shear mode, and a diffusive cross-sheet entropy mode. Once a consistent picture of in-sheet vorticity evolution is established, a number of analytical results, describing the

Global gradient estimates for divergence-type elliptic problems involving general nonlinear operators

Science.gov (United States)

Cho, Yumi

2018-05-01

We study nonlinear elliptic problems with nonstandard growth and ellipticity related to an N-function. We establish global Calderón-Zygmund estimates of the weak solutions in the framework of Orlicz spaces over bounded non-smooth domains. Moreover, we prove a global regularity result for asymptotically regular problems which are getting close to the regular problems considered, when the gradient variable goes to infinity.

Piecewise-linear and bilinear approaches to nonlinear differential equations approximation problem of computational structural mechanics

OpenAIRE

Leibov Roman

2017-01-01

This paper presents a bilinear approach to nonlinear differential equations system approximation problem. Sometimes the nonlinear differential equations right-hand sides linearization is extremely difficult or even impossible. Then piecewise-linear approximation of nonlinear differential equations can be used. The bilinear differential equations allow to improve piecewise-linear differential equations behavior and reduce errors on the border of different linear differential equations systems ...

Quasi-stability of a vector trajectorial problem with non-linear partial criteria

Directory of Open Access Journals (Sweden)

Vladimir A. Emelichev

2003-10-01

Full Text Available Multi-objective (vector combinatorial problem of finding the Pareto set with four kinds of non-linear partial criteria is considered. Necessary and sufficient conditions of that kind of stability of the problem (quasi-stability are obtained. The problem is a discrete analogue of the lower semicontinuity by Hausdorff of the optimal mapping. Mathematics Subject Classification 2000: 90C10, 90C05, 90C29, 90C31.

Determination of Nonlinear Stiffness Coefficients for Finite Element Models with Application to the Random Vibration Problem

Science.gov (United States)

Muravyov, Alexander A.

1999-01-01

In this paper, a method for obtaining nonlinear stiffness coefficients in modal coordinates for geometrically nonlinear finite-element models is developed. The method requires application of a finite-element program with a geometrically non- linear static capability. The MSC/NASTRAN code is employed for this purpose. The equations of motion of a MDOF system are formulated in modal coordinates. A set of linear eigenvectors is used to approximate the solution of the nonlinear problem. The random vibration problem of the MDOF nonlinear system is then considered. The solutions obtained by application of two different versions of a stochastic linearization technique are compared with linear and exact (analytical) solutions in terms of root-mean-square (RMS) displacements and strains for a beam structure.

A frictional contact problem with damage and adhesion for an electro elastic-viscoplastic body

Directory of Open Access Journals (Sweden)

Adel Aissaoui

2014-01-01

Full Text Available We consider a quasistatic frictional contact problem for an electro elastic-viscopalastic body with damage and adhestion. The contact is modelled with normal compliance. The adhesion of the contact surfaces is taken into account and modelled by a surface variable. We derive variational formulation for the model which is in the form of a system involving the displacement field, the electric potential field, the damage field and the adhesion field. We prove the existence of a unique weak solution to the problem. The proof is based on arguments of time-dependent variational inequalities, parabolic inequalities, differential equations and fixed point.

A Class of time-fractional hemivariational inequalities with application to frictional contact problem

Science.gov (United States)

Zeng, Shengda; Migórski, Stanisław

2018-03-01

In this paper a class of elliptic hemivariational inequalities involving the time-fractional order integral operator is investigated. Exploiting the Rothe method and using the surjectivity of multivalued pseudomonotone operators, a result on existence of solution to the problem is established. Then, this abstract result is applied to provide a theorem on the weak solvability of a fractional viscoelastic contact problem. The process is quasistatic and the constitutive relation is modeled with the fractional Kelvin-Voigt law. The friction and contact conditions are described by the Clarke generalized gradient of nonconvex and nonsmooth functionals. The variational formulation of this problem leads to a fractional hemivariational inequality.

Nonlinear Preconditioning and its Application in Multicomponent Problems

KAUST Repository

Liu, Lulu

2015-12-07

The Multiplicative Schwarz Preconditioned Inexact Newton (MSPIN) algorithm is presented as a complement to Additive Schwarz Preconditioned Inexact Newton (ASPIN). At an algebraic level, ASPIN and MSPIN are variants of the same strategy to improve the convergence of systems with unbalanced nonlinearities; however, they have natural complementarity in practice. MSPIN is naturally based on partitioning of degrees of freedom in a nonlinear PDE system by field type rather than by subdomain, where a modest factor of concurrency can be sacrificed for physically motivated convergence robustness. ASPIN, originally introduced for decompositions into subdomains, is natural for high concurrency and reduction of global synchronization. The ASPIN framework, as an option for the outermost solver, successfully handles strong nonlinearities in computational fluid dynamics, but is barely explored for the highly nonlinear models of complex multiphase flow with capillarity, heterogeneity, and complex geometry. In this dissertation, the fully implicit ASPIN method is demonstrated for a finite volume discretization based on incompressible two-phase reservoir simulators in the presence of capillary forces and gravity. Numerical experiments show that the number of global nonlinear iterations is not only scalable with respect to the number of processors, but also significantly reduced compared with the standard inexact Newton method with a backtracking technique. Moreover, the ASPIN method, in contrast with the IMPES method, saves overall execution time because of the savings in timestep size. We consider the additive and multiplicative types of inexact Newton algorithms in the field-split context, and we augment the classical convergence theory of ASPIN for the multiplicative case. Moreover, we provide the convergence analysis of the MSPIN algorithm. Under suitable assumptions, it is shown that MSPIN is locally convergent, and desired superlinear or even quadratic convergence can be

Fermat collocation method for the solutions of nonlinear system of second order boundary value problems

Directory of Open Access Journals (Sweden)

Salih Yalcinbas

2016-01-01

Full Text Available In this study, a numerical approach is proposed to obtain approximate solutions of nonlinear system of second order boundary value problem. This technique is essentially based on the truncated Fermat series and its matrix representations with collocation points. Using the matrix method, we reduce the problem system of nonlinear algebraic equations. Numerical examples are also given to demonstrate the validity and applicability of the presented technique. The method is easy to implement and produces accurate results.

Domain decomposition based iterative methods for nonlinear elliptic finite element problems

Energy Technology Data Exchange (ETDEWEB)

Cai, X.C. [Univ. of Colorado, Boulder, CO (United States)

1994-12-31

The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.

Numerical solution of non-linear diffusion problems

International Nuclear Information System (INIS)

Carmen, A. del; Ferreri, J.C.

1998-01-01

This paper presents a method for the numerical solution of non-linear diffusion problems using finite-differences in moving grids. Due to the presence of steep fronts in the solution domain and to the presence of advective terms originating in the grid movement, an implicit TVD scheme, first order in time and second order in space has been developed. Some algebraic details of the derivation are given. Results are shown for the pure advection of a scalar as a test case and an example dealing with the slow spreading of viscous fluids over plane surfaces. The agreement between numerical and analytical solutions is excellent. (author). 8 refs., 3 figs

Iterative Runge–Kutta-type methods for nonlinear ill-posed problems

International Nuclear Information System (INIS)

Böckmann, C; Pornsawad, P

2008-01-01

We present a regularization method for solving nonlinear ill-posed problems by applying the family of Runge–Kutta methods to an initial value problem, in particular, to the asymptotical regularization method. We prove that the developed iterative regularization method converges to a solution under certain conditions and with a general stopping rule. Some particular iterative regularization methods are numerically implemented. Numerical results of the examples show that the developed Runge–Kutta-type regularization methods yield stable solutions and that particular implicit methods are very efficient in saving iteration steps

On the quantum inverse problem for a new type of nonlinear Schroedinger equation for Alfven waves in plasma

International Nuclear Information System (INIS)

Sen, S.; Roy Chowdhury, A.

1989-06-01

The nonlinear Alfven waves are governed by the Vector Derivative nonlinear Schroedinger (VDNLS) equation, which for parallel or quasi parallel propagation reduces to the Derivative Nonlinear Schroedinger (DNLS) equation for the circularly polarized waves. We have formulated the Quantum Inverse problem for a new type of Nonlinear Schroedinger Equation which has many properties similar to the usual NLS problem but the structure of classical and quantum R matrix are distinctly different. The commutation rules of the scattering data are obtained and the Algebraic Bethe Ansatz is formulated to derive the eigenvalue equation for the energy of the excited states. 10 refs

An Algorithmic Comparison of the Hyper-Reduction and the Discrete Empirical Interpolation Method for a Nonlinear Thermal Problem

Directory of Open Access Journals (Sweden)

Felix Fritzen

2018-02-01

Full Text Available A novel algorithmic discussion of the methodological and numerical differences of competing parametric model reduction techniques for nonlinear problems is presented. First, the Galerkin reduced basis (RB formulation is presented, which fails at providing significant gains with respect to the computational efficiency for nonlinear problems. Renowned methods for the reduction of the computing time of nonlinear reduced order models are the Hyper-Reduction and the (Discrete Empirical Interpolation Method (EIM, DEIM. An algorithmic description and a methodological comparison of both methods are provided. The accuracy of the predictions of the hyper-reduced model and the (DEIM in comparison to the Galerkin RB is investigated. All three approaches are applied to a simple uncertainty quantification of a planar nonlinear thermal conduction problem. The results are compared to computationally intense finite element simulations.

Polyharmonic boundary value problems positivity preserving and nonlinear higher order elliptic equations in bounded domains

CERN Document Server

Gazzola, Filippo; Sweers, Guido

2010-01-01

This monograph covers higher order linear and nonlinear elliptic boundary value problems in bounded domains, mainly with the biharmonic or poly-harmonic operator as leading principal part. Underlying models and, in particular, the role of different boundary conditions are explained in detail. As for linear problems, after a brief summary of the existence theory and Lp and Schauder estimates, the focus is on positivity or - since, in contrast to second order equations, a general form of a comparison principle does not exist - on “near positivity.” The required kernel estimates are also presented in detail. As for nonlinear problems, several techniques well-known from second order equations cannot be utilized and have to be replaced by new and different methods. Subcritical, critical and supercritical nonlinearities are discussed and various existence and nonexistence results are proved. The interplay with the positivity topic from the first part is emphasized and, moreover, a far-reaching Gidas-Ni-Nirenbe...

On the Cauchy problem for nonlinear Schrödinger equations with rotation

KAUST Repository

Antonelli, Paolo; Marahrens, Daniel; Sparber, Christof

2011-01-01

We consider the Cauchy problem for (energy-subcritical) nonlinear Schrödinger equations with sub-quadratic external potentials and an additional angular momentum rotation term. This equation is a well-known model for superuid quantum gases in rotating traps. We prove global existence (in the energy space) for defocusing nonlinearities without any restriction on the rotation frequency, generalizing earlier results given in [11, 12]. Moreover, we find that the rotation term has a considerable in fiuence in proving finite time blow-up in the focusing case.

On the Cauchy problem for nonlinear Schrödinger equations with rotation

KAUST Repository

Antonelli, Paolo

2011-10-01

We consider the Cauchy problem for (energy-subcritical) nonlinear Schrödinger equations with sub-quadratic external potentials and an additional angular momentum rotation term. This equation is a well-known model for superuid quantum gases in rotating traps. We prove global existence (in the energy space) for defocusing nonlinearities without any restriction on the rotation frequency, generalizing earlier results given in [11, 12]. Moreover, we find that the rotation term has a considerable in fiuence in proving finite time blow-up in the focusing case.

Inequivalence of interior and exterior dynamical problems

International Nuclear Information System (INIS)

Santilli, R.M.

1991-09-01

We begin a series of notes with the review of the historical distinction by Lagrange, Hamilton, Jacobi and other Founding Fathers of analytic dynamics, between the exteriordynamical problem, consisting of motion in vacuum under action-at-a-distance interactions, and the interior dynamical problem, consisting of motion within a resistive medium with the additional presence of contact, nonlinear, nonlocal and nonhamiltonian internal forces. After recalling some of the historical reasons that led to the contemporary, virtually complete restriction of research to the exterior problem, we show that the interior dynamical problem cannot be reduced to the exterior one. This establishes the open character of the central objective of these notes: the identification of the space-time symmetries and relativities that are applicable to interior, nonlinear, nonlocal and nonhamiltonian systems. (author). 29 refs

Existence of bounded solutions of Neumann problem for a nonlinear degenerate elliptic equation

Directory of Open Access Journals (Sweden)

Salvatore Bonafede

2017-10-01

Full Text Available We prove the existence of bounded solutions of Neumann problem for nonlinear degenerate elliptic equations of second order in divergence form. We also study some properties as the Phragmen-Lindelof property and the asymptotic behavior of the solutions of Dirichlet problem associated to our equation in an unbounded domain.

Criminal Justice Contact, Stressors, and Obesity-Related Health Problems Among Black Adults in the USA.

Science.gov (United States)

Archibald, Paul C; Parker, Lauren; Thorpe, Roland

2018-04-01

Criminal justice contact-defined as lifetime arrest, parole, or incarceration, seems to exacerbate chronic conditions, and those who are most likely to have had contact with the criminal justice system, such as Black adults, often already have pre-existing disproportionately high rates of stress and chronic conditions due to the social determinants of health that affect underrepresented minorities. Findings from this study suggest that there is a mechanism that links the stressors among Black adults manifested by such factors as family, financial, neighborhood, and personal problems with criminal justice contact to obesity-related health status. Using the National Survey of American Life (NSAL), modified Poisson regression analyses were used to determine the association between criminal justice contact, stressors, and obesity-related health problems among a national sample of Black adults (n = 5008). In the full model, the odds of experiencing obesity-related health problems for Black adults who had criminal justice contact was reduced (PR, 1.23 to 1.14) and not statistically significant. Black adults who reported experiencing family stressors (PR, 1.21; 95% CI, 1.08, 1.36), financial stressors (PR, 1.30; 95% CI, 1.16, 1.47), and personal stressors (PR, 1.16; 95% CI, 1.02, 1.31) were statistically significant and higher than those who reported not experiencing any of these stressors; neighborhood stressors was not statistically significant. The evidence suggests a relationship between the stressors associated with criminal justice contact and obesity-related health status. These findings emphasize the need to further explore the family, financial, and personal stressors for Black adults with criminal justice contact in order to further our understanding of their obesity-related health problems.ᅟ.

Active-Set Reduced-Space Methods with Nonlinear Elimination for Two-Phase Flow Problems in Porous Media

KAUST Repository

Yang, Haijian

2016-07-26

Fully implicit methods are drawing more attention in scientific and engineering applications due to the allowance of large time steps in extreme-scale simulations. When using a fully implicit method to solve two-phase flow problems in porous media, one major challenge is the solution of the resultant nonlinear system at each time step. To solve such nonlinear systems, traditional nonlinear iterative methods, such as the class of the Newton methods, often fail to achieve the desired convergent rate due to the high nonlinearity of the system and/or the violation of the boundedness requirement of the saturation. In the paper, we reformulate the two-phase model as a variational inequality that naturally ensures the physical feasibility of the saturation variable. The variational inequality is then solved by an active-set reduced-space method with a nonlinear elimination preconditioner to remove the high nonlinear components that often causes the failure of the nonlinear iteration for convergence. To validate the effectiveness of the proposed method, we compare it with the classical implicit pressure-explicit saturation method for two-phase flow problems with strong heterogeneity. The numerical results show that our nonlinear solver overcomes the often severe limits on the time step associated with existing methods, results in superior convergence performance, and achieves reduction in the total computing time by more than one order of magnitude.

Active-Set Reduced-Space Methods with Nonlinear Elimination for Two-Phase Flow Problems in Porous Media

KAUST Repository

Yang, Haijian; Yang, Chao; Sun, Shuyu

2016-01-01

Fully implicit methods are drawing more attention in scientific and engineering applications due to the allowance of large time steps in extreme-scale simulations. When using a fully implicit method to solve two-phase flow problems in porous media, one major challenge is the solution of the resultant nonlinear system at each time step. To solve such nonlinear systems, traditional nonlinear iterative methods, such as the class of the Newton methods, often fail to achieve the desired convergent rate due to the high nonlinearity of the system and/or the violation of the boundedness requirement of the saturation. In the paper, we reformulate the two-phase model as a variational inequality that naturally ensures the physical feasibility of the saturation variable. The variational inequality is then solved by an active-set reduced-space method with a nonlinear elimination preconditioner to remove the high nonlinear components that often causes the failure of the nonlinear iteration for convergence. To validate the effectiveness of the proposed method, we compare it with the classical implicit pressure-explicit saturation method for two-phase flow problems with strong heterogeneity. The numerical results show that our nonlinear solver overcomes the often severe limits on the time step associated with existing methods, results in superior convergence performance, and achieves reduction in the total computing time by more than one order of magnitude.

Spectral methods for a nonlinear initial value problem involving pseudo differential operators

International Nuclear Information System (INIS)

Pasciak, J.E.

1982-01-01

Spectral methods (Fourier methods) for approximating the solution of a nonlinear initial value problem involving pseudo differential operators are defined and analyzed. A semidiscrete approximation to the nonlinear equation based on an L 2 projection is described. The semidiscrete L 2 approximation is shown to be a priori stable and convergent under sufficient decay and smoothness assumptions on the initial data. It is shown that the semidiscrete method converges with infinite order, that is, higher order decay and smoothness assumptions imply higher order error bounds. Spectral schemes based on spacial collocation are also discussed

APPLICATION OF PARAMETER CONTINUATION METHOD FOR INVESTIGATION OF VIBROIMPACT SYSTEMS DYNAMIC BEHAVIOUR. PROBLEM STATE. SHORT SURVEY OF WORLD SCIENTIFIC LITERATURE

Directory of Open Access Journals (Sweden)

V.A. Bazhenov

2014-12-01

Full Text Available Authors in their works study vibroimpact system dynamic behaviour by numerical parametric continuation technique combined with shooting and Newton-Raphson’s methods. The technique is adapted to two-mass two-degree-of-freedom vibroimpact system under periodic excitation. Impact is simulated by nonlinear contact interaction force based on Hertz’s contact theory. Stability or instability of obtained periodic solutions is determined by monodromy matrix eigenvalues (multipliers based on Floquet’s theory. In the present paper we describe the state of problem of parameter continuation method using for nonlinear tasks solution. Also we give the short survey of numerous contemporary literature in English and Russian about parameter continuation method application for nonlinear problems. This method is applied for vibroimpact problem solving more rarely because of the difficulties connected with repeated impacts.

Initial-value problem for the Gardner equation applied to nonlinear internal waves

Science.gov (United States)

Rouvinskaya, Ekaterina; Kurkina, Oxana; Kurkin, Andrey; Talipova, Tatiana; Pelinovsky, Efim

2017-04-01

The Gardner equation is a fundamental mathematical model for the description of weakly nonlinear weakly dispersive internal waves, when cubic nonlinearity cannot be neglected. Within this model coefficients of quadratic and cubic nonlinearity can both be positive as well as negative, depending on background conditions of the medium, where waves propagate (sea water density stratification, shear flow profile) [Rouvinskaya et al., 2014, Kurkina et al., 2011, 2015]. For the investigation of weakly dispersive behavior in the framework of nondimensional Gardner equation with fixed (positive) sign of quadratic nonlinearity and positive or negative cubic nonlinearity {eq1} partial η/partial t+6η( {1± η} )partial η/partial x+partial ^3η/partial x^3=0, } the series of numerical experiments of initial-value problem was carried out for evolution of a bell-shaped impulse of negative polarity (opposite to the sign of quadratic nonlinear coefficient): {eq2} η(x,t=0)=-asech2 ( {x/x0 } ), for which amplitude a and width x0 was varied. Similar initial-value problem was considered in the paper [Trillo et al., 2016] for the Korteweg - de Vries equation. For the Gardner equation with different signs of cubic nonlinearity the initial-value problem for piece-wise constant initial condition was considered in detail in [Grimshaw et al., 2002, 2010]. It is widely known, for example, [Pelinovsky et al., 2007], that the Gardner equation (1) with negative cubic nonlinearity has a family of classic solitary wave solutions with only positive polarity,and with limiting amplitude equal to 1. Therefore evolution of impulses (2) of negative polarity (whose amplitudes a were varied from 0.1 to 3, and widths at the level of a/2 were equal to triple width of solitons with the same amplitude for a 1) was going on a universal scenario with the generation of nonlinear Airy wave. For the Gardner equation (1) with the positive cubic nonlinearity coefficient there exist two one-parametric families of

A quadratic approximation-based algorithm for the solution of multiparametric mixed-integer nonlinear programming problems

KAUST Repository

Domínguez, Luis F.

2012-06-25

An algorithm for the solution of convex multiparametric mixed-integer nonlinear programming problems arising in process engineering problems under uncertainty is introduced. The proposed algorithm iterates between a multiparametric nonlinear programming subproblem and a mixed-integer nonlinear programming subproblem to provide a series of parametric upper and lower bounds. The primal subproblem is formulated by fixing the integer variables and solved through a series of multiparametric quadratic programming (mp-QP) problems based on quadratic approximations of the objective function, while the deterministic master subproblem is formulated so as to provide feasible integer solutions for the next primal subproblem. To reduce the computational effort when infeasibilities are encountered at the vertices of the critical regions (CRs) generated by the primal subproblem, a simplicial approximation approach is used to obtain CRs that are feasible at each of their vertices. The algorithm terminates when there does not exist an integer solution that is better than the one previously used by the primal problem. Through a series of examples, the proposed algorithm is compared with a multiparametric mixed-integer outer approximation (mp-MIOA) algorithm to demonstrate its computational advantages. © 2012 American Institute of Chemical Engineers (AIChE).

A Simple FEM Formulation Applied to Nonlinear Problems of Impact with Thermomechanical Coupling

Directory of Open Access Journals (Sweden)

João Paulo de Barros Cavalcante

Full Text Available Abstract The thermal effects of problems involving deformable structures are essential to describe the behavior of materials in feasible terms. Verifying the transformation of mechanical energy into heat it is possible to predict the modifications of mechanical properties of materials due to its temperature changes. The current paper presents the numerical development of a finite element method suitable for nonlinear structures coupled with thermomechanical behavior; including impact problems. A simple and effective alternative formulation is presented, called FEM positional, to deal with the dynamic nonlinear systems. The developed numerical is based on the minimum potential energy written in terms of nodal positions instead of displacements. The effects of geometrical, material and thermal nonlinearities are considered. The thermodynamically consistent formulation is based on the laws of thermodynamics and the Helmholtz free-energy, used to describe the thermoelastic and the thermoplastic behaviors. The coupled thermomechanical model can result in secondary effects that cause redistributions of internal efforts, depending on the history of deformation and material properties. The numerical results of the proposed formulation are compared with examples found in the literature.

Adaptive discretizations for the choice of a Tikhonov regularization parameter in nonlinear inverse problems

International Nuclear Information System (INIS)

Kaltenbacher, Barbara; Kirchner, Alana; Vexler, Boris

2011-01-01

Parameter identification problems for partial differential equations usually lead to nonlinear inverse problems. A typical property of such problems is their instability, which requires regularization techniques, like, e.g., Tikhonov regularization. The main focus of this paper will be on efficient methods for determining a suitable regularization parameter by using adaptive finite element discretizations based on goal-oriented error estimators. A well-established method for the determination of a regularization parameter is the discrepancy principle where the residual norm, considered as a function i of the regularization parameter, should equal an appropriate multiple of the noise level. We suggest to solve the resulting scalar nonlinear equation by an inexact Newton method, where in each iteration step, a regularized problem is solved at a different discretization level. The proposed algorithm is an extension of the method suggested in Griesbaum A et al (2008 Inverse Problems 24 025025) for linear inverse problems, where goal-oriented error estimators for i and its derivative are used for adaptive refinement strategies in order to keep the discretization level as coarse as possible to save computational effort but fine enough to guarantee global convergence of the inexact Newton method. This concept leads to a highly efficient method for determining the Tikhonov regularization parameter for nonlinear ill-posed problems. Moreover, we prove that with the so-obtained regularization parameter and an also adaptively discretized Tikhonov minimizer, usual convergence and regularization results from the continuous setting can be recovered. As a matter of fact, it is shown that it suffices to use stationary points of the Tikhonov functional. The efficiency of the proposed method is demonstrated by means of numerical experiments. (paper)

Stokes phenomena and monodromy deformation problem for nonlinear Schrodinger equation

International Nuclear Information System (INIS)

Chowdury, A.R.; Naskar, M.

1986-01-01

Following Flaschka and Newell, the inverse problem for Painleve IV is formulated with the help of similarity variables. The Painleve IV arises as the eliminant of the two second-order ordinary differential equations originating from the nonlinear Schrodinger equation. Asymptotic expansions are obtained near the singularities at zero and infinity of the complex eigenvalue plane. The corresponding analysis then displays the Stokes phenomena. The monodromy matrices connecting the solution Y /sub j/ in the sector S /sub j/ to that in S /sub j+1/ are fixed in structure by the imposition of certain conditions. It is then shown that a deformation keeping the monodromy data fixed leads to the nonlinear Schrodinger equation. While Flaschka and Newell did not make any absolute determination of the Stokes parameters, the present approach yields the values of the Stokes parameters in an explicit way, which in turn can determine the matrix connecting the solutions near zero and infinity. Finally, it is shown that the integral equation originating from the analyticity and asymptotic nature of the problem leads to the similarity solution previously determined by Boiti and Pampinelli

A Kind of Nonlinear Programming Problem Based on Mixed Fuzzy Relation Equations Constraints

Science.gov (United States)

Li, Jinquan; Feng, Shuang; Mi, Honghai

In this work, a kind of nonlinear programming problem with non-differential objective function and under the constraints expressed by a system of mixed fuzzy relation equations is investigated. First, some properties of this kind of optimization problem are obtained. Then, a polynomial-time algorithm for this kind of optimization problem is proposed based on these properties. Furthermore, we show that this algorithm is optimal for the considered optimization problem in this paper. Finally, numerical examples are provided to illustrate our algorithms.

A finite element method with contact for tensile analysis in fuel rods

International Nuclear Information System (INIS)

Tanajura, C.A.S.; Galeao, A.C.N.R.

1987-01-01

Elements for mechanical analysis of fuel rod of a PWR type reactor, are presented. The rod, consists basically in a cylindrical coating of zircalloy which contains pilling of UO 2 pellets, is submitted to strong internal and external pressures, intense temperature gradients and neutron flux. These conditions lead several phenomena in the pellet (swelling, fracture, densification, creep) and in the cladding (embrittlement, corrosion, creep) which undergo deformations leading them to contact the restriction for the interpenetration is included in the problem without restriction by Lagrange multipliers. Considering a non-linear problem, due to the surface of contact to be not known a priori, the numerical solutions were obtained using the finite element method. (M.C.K.) [pt

Approximate stresses in 2-D flat elastic contact fretting problems

Science.gov (United States)

Urban, Michael Rene

Fatigue results from the cyclic loading of a solid body. If the body subject to fatigue is in contact with another body and relative sliding motion occurs between these two bodies, then rubbing surface damage can accelerate fatigue failure. The acceleration of fatigue failure is especially important if the relative motion between the two bodies results in surface damage without excessive surface removal via wear. The situation just described is referred to as fretting fatigue. Understanding of fretting fatigue is greatly enhanced if the stress state associated with fretting can be characterized. For Hertzian contact, this can readily be done. Unfortunately, simple stress formulae are not available for flat body contact. The primary result of the present research is the development of a new, reasonably accurate, approximate closed form expression for 2-dimensional contact stresses which has been verified using finite element modeling. This expression is also combined with fracture mechanics to provide a simple method of determining when a crack is long enough to no longer be affected by the contact stress field. Lower bounds on fatigue life can then easily be calculated using fracture mechanics. This closed form expression can also be used to calculate crack propagation within the contact stress field. The problem of determining the cycles required to generate an initial crack and what to choose as an initial crack size is unresolved as it is in non-fretting fatigue.

Contact Problem of Disk on Shaft Fixed by Induction Shrink Fit

Czech Academy of Sciences Publication Activity Database

Ulrych, B.; Kotlan, V.; Doležel, Ivo

2012-01-01

Roč. 88, 12B (2012), s. 32-34 ISSN 0033-2097 Institutional research plan: CEZ:AV0Z20760514 Keywords : induction shrink fit * contact problem * transfer of torque Subject RIV: BA - General Mathematics Impact factor: 0.244, year: 2011

Numerical methods for the design of large-scale nonlinear discrete ill-posed inverse problems

International Nuclear Information System (INIS)

Haber, E; Horesh, L; Tenorio, L

2010-01-01

Design of experiments for discrete ill-posed problems is a relatively new area of research. While there has been some limited work concerning the linear case, little has been done to study design criteria and numerical methods for ill-posed nonlinear problems. We present an algorithmic framework for nonlinear experimental design with an efficient numerical implementation. The data are modeled as indirect, noisy observations of the model collected via a set of plausible experiments. An inversion estimate based on these data is obtained by a weighted Tikhonov regularization whose weights control the contribution of the different experiments to the data misfit term. These weights are selected by minimization of an empirical estimate of the Bayes risk that is penalized to promote sparsity. This formulation entails a bilevel optimization problem that is solved using a simple descent method. We demonstrate the viability of our design with a problem in electromagnetic imaging based on direct current resistivity and magnetotelluric data

A dynamic unilateral contact problem with adhesion and friction in viscoelasticity

Science.gov (United States)

Cocou, Marius; Schryve, Mathieu; Raous, Michel

2010-08-01

The aim of this paper is to study an interaction law coupling recoverable adhesion, friction and unilateral contact between two viscoelastic bodies of Kelvin-Voigt type. A dynamic contact problem with adhesion and nonlocal friction is considered and its variational formulation is written as the coupling between an implicit variational inequality and a parabolic variational inequality describing the evolution of the intensity of adhesion. The existence and approximation of variational solutions are analysed, based on a penalty method, some abstract results and compactness properties. Finally, some numerical examples are presented.

Optimization for nonlinear inverse problem

International Nuclear Information System (INIS)

Boyadzhiev, G.; Brandmayr, E.; Pinat, T.; Panza, G.F.

2007-06-01

The nonlinear inversion of geophysical data in general does not yield a unique solution, but a single model, representing the investigated field, is preferred for an easy geological interpretation of the observations. The analyzed region is constituted by a number of sub-regions where the multi-valued nonlinear inversion is applied, which leads to a multi-valued solution. Therefore, combining the values of the solution in each sub-region, many acceptable models are obtained for the entire region and this complicates the geological interpretation of geophysical investigations. In this paper are presented new methodologies, capable to select one model, among all acceptable ones, that satisfies different criteria of smoothness in the explored space of solutions. In this work we focus on the non-linear inversion of surface waves dispersion curves, which gives structural models of shear-wave velocity versus depth, but the basic concepts have a general validity. (author)

A model-adaptivity method for the solution of Lennard-Jones based adhesive contact problems

Science.gov (United States)

Ben Dhia, Hachmi; Du, Shuimiao

2018-05-01

The surface micro-interaction model of Lennard-Jones (LJ) is used for adhesive contact problems (ACP). To address theoretical and numerical pitfalls of this model, a sequence of partitions of contact models is adaptively constructed to both extend and approximate the LJ model. It is formed by a combination of the LJ model with a sequence of shifted-Signorini (or, alternatively, -Linearized-LJ) models, indexed by a shift parameter field. For each model of this sequence, a weak formulation of the associated local ACP is developed. To track critical localized adhesive areas, a two-step strategy is developed: firstly, a macroscopic frictionless (as first approach) linear-elastic contact problem is solved once to detect contact separation zones. Secondly, at each shift-adaptive iteration, a micro-macro ACP is re-formulated and solved within the multiscale Arlequin framework, with significant reduction of computational costs. Comparison of our results with available analytical and numerical solutions shows the effectiveness of our global strategy.

Adaptive wavelet collocation methods for initial value boundary problems of nonlinear PDE's

Science.gov (United States)

Cai, Wei; Wang, Jian-Zhong

1993-01-01

We have designed a cubic spline wavelet decomposition for the Sobolev space H(sup 2)(sub 0)(I) where I is a bounded interval. Based on a special 'point-wise orthogonality' of the wavelet basis functions, a fast Discrete Wavelet Transform (DWT) is constructed. This DWT transform will map discrete samples of a function to its wavelet expansion coefficients in O(N log N) operations. Using this transform, we propose a collocation method for the initial value boundary problem of nonlinear PDE's. Then, we test the efficiency of the DWT transform and apply the collocation method to solve linear and nonlinear PDE's.

Frictional Heating During Sliding of two Semi-Spaces with Arbitrary Thermal Nonlinearity

Directory of Open Access Journals (Sweden)

Och Ewa

2014-12-01

Full Text Available Analytical and numerical solution for transient thermal problems of friction were presented for semi limited bodies made from thermosensitive materials in which coefficient of thermal conductivity and specific heat arbitrarily depend on the temperature (materials with arbitrary non-linearity. With the constant power of friction assumption and imperfect thermal contact linearization of nonlinear problems formulated initial-boundary thermal conductivity, using Kirchhoff transformation is partial. In order to complete linearization, method of successive approximations was used. On the basis of obtained solutions a numerical analysis of two friction systems in which one element is constant (cermet FMC-845 and another is variable (grey iron ChNMKh or aluminum-based composite alloy AL MMC was conducted

Parallel supercomputing: Advanced methods, algorithms, and software for large-scale linear and nonlinear problems

Energy Technology Data Exchange (ETDEWEB)

Carey, G.F.; Young, D.M.

1993-12-31

The program outlined here is directed to research on methods, algorithms, and software for distributed parallel supercomputers. Of particular interest are finite element methods and finite difference methods together with sparse iterative solution schemes for scientific and engineering computations of very large-scale systems. Both linear and nonlinear problems will be investigated. In the nonlinear case, applications with bifurcation to multiple solutions will be considered using continuation strategies. The parallelizable numerical methods of particular interest are a family of partitioning schemes embracing domain decomposition, element-by-element strategies, and multi-level techniques. The methods will be further developed incorporating parallel iterative solution algorithms with associated preconditioners in parallel computer software. The schemes will be implemented on distributed memory parallel architectures such as the CRAY MPP, Intel Paragon, the NCUBE3, and the Connection Machine. We will also consider other new architectures such as the Kendall-Square (KSQ) and proposed machines such as the TERA. The applications will focus on large-scale three-dimensional nonlinear flow and reservoir problems with strong convective transport contributions. These are legitimate grand challenge class computational fluid dynamics (CFD) problems of significant practical interest to DOE. The methods developed and algorithms will, however, be of wider interest.

Nonlinear Elliptic Differential Equations with Multivalued Nonlinearities

Indian Academy of Sciences (India)

In this paper we study nonlinear elliptic boundary value problems with monotone and nonmonotone multivalued nonlinearities. First we consider the case of monotone nonlinearities. In the first result we assume that the multivalued nonlinearity is defined on all R R . Assuming the existence of an upper and of a lower ...

A New Spectral Local Linearization Method for Nonlinear Boundary Layer Flow Problems

Directory of Open Access Journals (Sweden)

S. S. Motsa

2013-01-01

Full Text Available We propose a simple and efficient method for solving highly nonlinear systems of boundary layer flow problems with exponentially decaying profiles. The algorithm of the proposed method is based on an innovative idea of linearizing and decoupling the governing systems of equations and reducing them into a sequence of subsystems of differential equations which are solved using spectral collocation methods. The applicability of the proposed method, hereinafter referred to as the spectral local linearization method (SLLM, is tested on some well-known boundary layer flow equations. The numerical results presented in this investigation indicate that the proposed method, despite being easy to develop and numerically implement, is very robust in that it converges rapidly to yield accurate results and is more efficient in solving very large systems of nonlinear boundary value problems of the similarity variable boundary layer type. The accuracy and numerical stability of the SLLM can further be improved by using successive overrelaxation techniques.

Solvability of Static Contact Problems with Coulomb Friction for Orthotropic Material

Czech Academy of Sciences Publication Activity Database

Eck, C.; Jarušek, Jiří

2008-01-01

Roč. 93, č. 1 (2008), s. 93-104 ISSN 0374-3535 R&D Projects: GA AV ČR IAA1075402 Institutional research plan: CEZ:AV0Z10190503 Keywords : contact problem * Coulomb friction * orthotropic elasticity Subject RIV: BA - General Mathematics Impact factor: 1.277, year: 2008

Dynamic contact problem with adhesion and damage between thermo-electro-elasto-viscoplastic bodies

Science.gov (United States)

Hadj ammar, Tedjani; Saïdi, Abdelkader; Azeb Ahmed, Abdelaziz

2017-05-01

We study of a dynamic contact problem between two thermo-electro-elasto-viscoplastic bodies with damage and adhesion. The contact is frictionless and is modeled with normal compliance condition. We derive variational formulation for the model and prove an existence and uniqueness result of the weak solution. The proof is based on arguments of evolutionary variational inequalities, parabolic inequalities, differential equations, and fixed point theorem.

A frictional contact problem with wear involving elastic-viscoplastic materials with damage and thermal effects

Directory of Open Access Journals (Sweden)

Abdelmoumene Djabi

2015-05-01

Full Text Available We consider a mathematical problem for quasistatic contact between a thermo-elastic-viscoplastic body with damage and an obstacle. The contact is frictional and bilateral with a moving rigid foundation which results in the wear of the contacting surface. We employ the thermo-elasticviscoplastic with damage constitutive law for the material. The damage of the material caused by elastic deformations. The evolution of the damage is described by an inclusion of parabolic type. The problem is formulated as a coupled system of an elliptic variational inequality for the displacement, a parabolic variational inequality for the damage and the heat equation for the temperature. We establish a variational formulation for the model and we prove the existence of a unique weak solution to the problem. The proof is based on a classical existence and uniqueness result on parabolic inequalities, differential equations and fixed point arguments.

Special function solutions of a spectral problem for a nonlinear quantum oscillator

International Nuclear Information System (INIS)

Schulze-Halberg, A; Morris, J R

2012-01-01

We construct exact solutions of a spectral problem involving the Schrödinger equation for a nonlinear, one-parameter oscillator potential. In contrast to a previous analysis of the problem (Carinena et al 2007 Ann. Phys. 322 434–59), where solutions were given through a Rodrigues-type formula, our approach leads to closed-form representations of the solutions in terms of special functions, not containing any derivative operators. We show normalizability and orthogonality of our solutions, as well as correct reduction of the problem to the harmonic oscillator model, if the parameter in the potential gets close to zero. (paper)

Multi-level nonlinear diffusion acceleration method for multigroup transport k-Eigenvalue problems

International Nuclear Information System (INIS)

Anistratov, Dmitriy Y.

2011-01-01

The nonlinear diffusion acceleration (NDA) method is an efficient and flexible transport iterative scheme for solving reactor-physics problems. This paper presents a fast iterative algorithm for solving multigroup neutron transport eigenvalue problems in 1D slab geometry. The proposed method is defined by a multi-level system of equations that includes multigroup and effective one-group low-order NDA equations. The Eigenvalue is evaluated in the exact projected solution space of smallest dimensionality, namely, by solving the effective one- group eigenvalue transport problem. Numerical results that illustrate performance of the new algorithm are demonstrated. (author)

Lamé Parameter Estimation from Static Displacement Field Measurements in the Framework of Nonlinear Inverse Problems

DEFF Research Database (Denmark)

Hubmer, Simon; Sherina, Ekaterina; Neubauer, Andreas

2018-01-01

. The main result of this paper is the verification of a nonlinearity condition in an infinite dimensional Hilbert space context. This condition guarantees convergence of iterative regularization methods. Furthermore, numerical examples for recovery of the Lam´e parameters from displacement data simulating......We consider a problem of quantitative static elastography, the estimation of the Lam´e parameters from internal displacement field data. This problem is formulated as a nonlinear operator equation. To solve this equation, we investigate the Landweber iteration both analytically and numerically...... a static elastography experiment are presented....

Studies in nonlinear problems of energy

Energy Technology Data Exchange (ETDEWEB)

Matkowsky, B.J.

1992-07-01

Emphasis has been on combustion and flame propagation. The research program was on modeling, analysis and computation of combustion phenomena, with emphasis on transition from laminar to turbulent combustion. Nonlinear dynamics and pattern formation were investigated in the transition. Stability of combustion waves, and transitions to complex waves are described. Combustion waves possess large activation energies, so that chemical reactions are significant only in thin layers, or reaction zones. In limit of infinite activation energy, the zones shrink to moving surfaces, (fronts) which must be found during the analysis, so that (moving free boundary problems). The studies are carried out for limiting case with fronts, while the numerical studies are carried out for finite, though large, activation energy. Accurate resolution of the solution in the reaction zones is essential, otherwise false predictions of dynamics are possible. Since the the reaction zones move, adaptive pseudo-spectral methods were developed. The approach is based on a synergism of analytical and computational methods. The numerical computations build on and extend the analytical information. Furthermore, analytical solutions serve as benchmarks for testing the accuracy of the computation. Finally, ideas from analysis (singular perturbation theory) have induced new approaches to computations. The computational results suggest new analysis to be considered. Among the recent interesting results, was spatio-temporal chaos in combustion. One goal is extension of the adaptive pseudo-spectral methods to adaptive domain decomposition methods. Efforts have begun to develop such methods for problems with multiple reaction zones, corresponding to problems with more complex, and more realistic chemistry. Other topics included stochastics, oscillators, Rysteretic Josephson junctions, DC SQUID, Markov jumps, laser with saturable absorber, chemical physics, Brownian movement, combustion synthesis, etc.

Multigrid Methods for EHL Problems

Science.gov (United States)

Nurgat, Elyas; Berzins, Martin

1996-01-01

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of

A MODIFIED DECOMPOSITION METHOD FOR SOLVING NONLINEAR PROBLEM OF FLOW IN CONVERGING- DIVERGING CHANNEL

Directory of Open Access Journals (Sweden)

MOHAMED KEZZAR

2015-08-01

Full Text Available In this research, an efficient technique of computation considered as a modified decomposition method was proposed and then successfully applied for solving the nonlinear problem of the two dimensional flow of an incompressible viscous fluid between nonparallel plane walls. In fact this method gives the nonlinear term Nu and the solution of the studied problem as a power series. The proposed iterative procedure gives on the one hand a computationally efficient formulation with an acceleration of convergence rate and on the other hand finds the solution without any discretization, linearization or restrictive assumptions. The comparison of our results with those of numerical treatment and other earlier works shows clearly the higher accuracy and efficiency of the used Modified Decomposition Method.

Strongly nonlinear nonhomogeneous elliptic unilateral problems with L^1 data and no sign conditions

Directory of Open Access Journals (Sweden)

Elhoussine Azroul

2012-05-01

Full Text Available In this article, we prove the existence of solutions to unilateral problems involving nonlinear operators of the form: $$ Au+H(x,u,abla u=f $$ where $A$ is a Leray Lions operator from $W_0^{1,p(x}(Omega$ into its dual $W^{-1,p'(x}(Omega$ and $H(x,s,xi$ is the nonlinear term satisfying some growth condition but no sign condition. The right hand side $f$ belong to $L^1(Omega$.

Energy transfer in contact binary systems

International Nuclear Information System (INIS)

Robertson, J.A.

1980-01-01

A simple model for the transfer of energy by steady circulation within the envelope of a contact binary system is presented. The model describes the fully compressible, two-dimensional flow of a perfect gas within a rectangular region in a uniform gravitational field. The region is heated non-uniformly from below. Coriolis forces are neglected but the interaction of the circulation with convection is discussed briefly. Numerical solutions of the linearized equations of the problem are discussed in detail, and the results of some non-linear calculations are also presented. The influence of alternative boundary conditions is examined. (author)

Modelling the nonlinear behaviour of an underplatform damper test rig for turbine applications

Science.gov (United States)

Pesaresi, L.; Salles, L.; Jones, A.; Green, J. S.; Schwingshackl, C. W.

2017-02-01

Underplatform dampers (UPD) are commonly used in aircraft engines to mitigate the risk of high-cycle fatigue failure of turbine blades. The energy dissipated at the friction contact interface of the damper reduces the vibration amplitude significantly, and the couplings of the blades can also lead to significant shifts of the resonance frequencies of the bladed disk. The highly nonlinear behaviour of bladed discs constrained by UPDs requires an advanced modelling approach to ensure that the correct damper geometry is selected during the design of the turbine, and that no unexpected resonance frequencies and amplitudes will occur in operation. Approaches based on an explicit model of the damper in combination with multi-harmonic balance solvers have emerged as a promising way to predict the nonlinear behaviour of UPDs correctly, however rigorous experimental validations are required before approaches of this type can be used with confidence. In this study, a nonlinear analysis based on an updated explicit damper model having different levels of detail is performed, and the results are evaluated against a newly-developed UPD test rig. Detailed linear finite element models are used as input for the nonlinear analysis, allowing the inclusion of damper flexibility and inertia effects. The nonlinear friction interface between the blades and the damper is described with a dense grid of 3D friction contact elements which allow accurate capturing of the underlying nonlinear mechanism that drives the global nonlinear behaviour. The introduced explicit damper model showed a great dependence on the correct contact pressure distribution. The use of an accurate, measurement based, distribution, better matched the nonlinear dynamic behaviour of the test rig. Good agreement with the measured frequency response data could only be reached when the zero harmonic term (constant term) was included in the multi-harmonic expansion of the nonlinear problem, highlighting its importance

On the solutions of the dKP equation: the nonlinear Riemann Hilbert problem, longtime behaviour, implicit solutions and wave breaking

International Nuclear Information System (INIS)

Manakov, S V; Santini, P M

2008-01-01

We have recently solved the inverse scattering problem for one-parameter families of vector fields, and used this result to construct the formal solution of the Cauchy problem for a class of integrable nonlinear partial differential equations in multidimensions, including the second heavenly equation of Plebanski and the dispersionless Kadomtsev-Petviashvili (dKP) equation. We showed, in particular, that the associated inverse problems can be expressed in terms of nonlinear Riemann-Hilbert problems on the real axis. In this paper, we make use of the nonlinear Riemann-Hilbert problem of dKP (i) to construct the longtime behaviour of the solutions of its Cauchy problem; (ii) to characterize a class of implicit solutions; (iii) to elucidate the spectral mechanism causing the gradient catastrophe of localized solutions of dKP, at finite time as well as in the longtime regime, and the corresponding universal behaviours near breaking

On the solutions of the dKP equation: the nonlinear Riemann Hilbert problem, longtime behaviour, implicit solutions and wave breaking

Energy Technology Data Exchange (ETDEWEB)

Manakov, S V [Landau Institute for Theoretical Physics, Moscow (Russian Federation); Santini, P M [Dipartimento di Fisica, Universita di Roma ' La Sapienza' , and Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, Piazz.le Aldo Moro 2, I-00185 Rome (Italy)

2008-02-08

We have recently solved the inverse scattering problem for one-parameter families of vector fields, and used this result to construct the formal solution of the Cauchy problem for a class of integrable nonlinear partial differential equations in multidimensions, including the second heavenly equation of Plebanski and the dispersionless Kadomtsev-Petviashvili (dKP) equation. We showed, in particular, that the associated inverse problems can be expressed in terms of nonlinear Riemann-Hilbert problems on the real axis. In this paper, we make use of the nonlinear Riemann-Hilbert problem of dKP (i) to construct the longtime behaviour of the solutions of its Cauchy problem; (ii) to characterize a class of implicit solutions; (iii) to elucidate the spectral mechanism causing the gradient catastrophe of localized solutions of dKP, at finite time as well as in the longtime regime, and the corresponding universal behaviours near breaking.

Nonlinear Analysis and Modeling of Tires

Science.gov (United States)

Noor, Ahmed K.

1996-01-01

The objective of the study was to develop efficient modeling techniques and computational strategies for: (1) predicting the nonlinear response of tires subjected to inflation pressure, mechanical and thermal loads; (2) determining the footprint region, and analyzing the tire pavement contact problem, including the effect of friction; and (3) determining the sensitivity of the tire response (displacements, stresses, strain energy, contact pressures and contact area) to variations in the different material and geometric parameters. Two computational strategies were developed. In the first strategy the tire was modeled by using either a two-dimensional shear flexible mixed shell finite elements or a quasi-three-dimensional solid model. The contact conditions were incorporated into the formulation by using a perturbed Lagrangian approach. A number of model reduction techniques were applied to substantially reduce the number of degrees of freedom used in describing the response outside the contact region. The second strategy exploited the axial symmetry of the undeformed tire, and uses cylindrical coordinates in the development of three-dimensional elements for modeling each of the different parts of the tire cross section. Model reduction techniques are also used with this strategy.

Finite element solution of nonlinear eddy current problems with periodic excitation and its industrial applications.

Science.gov (United States)

Bíró, Oszkár; Koczka, Gergely; Preis, Kurt

2014-05-01

An efficient finite element method to take account of the nonlinearity of the magnetic materials when analyzing three-dimensional eddy current problems is presented in this paper. The problem is formulated in terms of vector and scalar potentials approximated by edge and node based finite element basis functions. The application of Galerkin techniques leads to a large, nonlinear system of ordinary differential equations in the time domain. The excitations are assumed to be time-periodic and the steady-state periodic solution is of interest only. This is represented either in the frequency domain as a finite Fourier series or in the time domain as a set of discrete time values within one period for each finite element degree of freedom. The former approach is the (continuous) harmonic balance method and, in the latter one, discrete Fourier transformation will be shown to lead to a discrete harmonic balance method. Due to the nonlinearity, all harmonics, both continuous and discrete, are coupled to each other. The harmonics would be decoupled if the problem were linear, therefore, a special nonlinear iteration technique, the fixed-point method is used to linearize the equations by selecting a time-independent permeability distribution, the so-called fixed-point permeability in each nonlinear iteration step. This leads to uncoupled harmonics within these steps. As industrial applications, analyses of large power transformers are presented. The first example is the computation of the electromagnetic field of a single-phase transformer in the time domain with the results compared to those obtained by traditional time-stepping techniques. In the second application, an advanced model of the same transformer is analyzed in the frequency domain by the harmonic balance method with the effect of the presence of higher harmonics on the losses investigated. Finally a third example tackles the case of direct current (DC) bias in the coils of a single-phase transformer.

State and parameter estimation in nonlinear systems as an optimal tracking problem

International Nuclear Information System (INIS)

Creveling, Daniel R.; Gill, Philip E.; Abarbanel, Henry D.I.

2008-01-01

In verifying and validating models of nonlinear processes it is important to incorporate information from observations in an efficient manner. Using the idea of synchronization of nonlinear dynamical systems, we present a framework for connecting a data signal with a model in a way that minimizes the required coupling yet allows the estimation of unknown parameters in the model. The need to evaluate unknown parameters in models of nonlinear physical, biophysical, and engineering systems occurs throughout the development of phenomenological or reduced models of dynamics. Our approach builds on existing work that uses synchronization as a tool for parameter estimation. We address some of the critical issues in that work and provide a practical framework for finding an accurate solution. In particular, we show the equivalence of this problem to that of tracking within an optimal control framework. This equivalence allows the application of powerful numerical methods that provide robust practical tools for model development and validation

Two-Polarisation Physical Model of Bowed Strings with Nonlinear Contact and Friction Forces, and Application to Gesture-Based Sound Synthesis

Directory of Open Access Journals (Sweden)

Charlotte Desvages

2016-05-01

Full Text Available Recent bowed string sound synthesis has relied on physical modelling techniques; the achievable realism and flexibility of gestural control are appealing, and the heavier computational cost becomes less significant as technology improves. A bowed string sound synthesis algorithm is designed, by simulating two-polarisation string motion, discretising the partial differential equations governing the string’s behaviour with the finite difference method. A globally energy balanced scheme is used, as a guarantee of numerical stability under highly nonlinear conditions. In one polarisation, a nonlinear contact model is used for the normal forces exerted by the dynamic bow hair, left hand fingers, and fingerboard. In the other polarisation, a force-velocity friction curve is used for the resulting tangential forces. The scheme update requires the solution of two nonlinear vector equations. The dynamic input parameters allow for simulating a wide range of gestures; some typical bow and left hand gestures are presented, along with synthetic sound and video demonstrations.

A demonstration of the improved efficiency of the canonical coordinates method using nonlinear combined heat and power economic dispatch problems

Science.gov (United States)

Chang, Hung-Chieh; Lin, Pei-Chun

2014-02-01

Economic dispatch is the short-term determination of the optimal output from a number of electricity generation facilities to meet the system load while providing power. As such, it represents one of the main optimization problems in the operation of electrical power systems. This article presents techniques to substantially improve the efficiency of the canonical coordinates method (CCM) algorithm when applied to nonlinear combined heat and power economic dispatch (CHPED) problems. The improvement is to eliminate the need to solve a system of nonlinear differential equations, which appears in the line search process in the CCM algorithm. The modified algorithm was tested and the analytical solution was verified using nonlinear CHPED optimization problems, thereby demonstrating the effectiveness of the algorithm. The CCM methods proved numerically stable and, in the case of nonlinear programs, produced solutions with unprecedented accuracy within a reasonable time.

Group-invariant solutions of nonlinear elastodynamic problems of plates and shells

International Nuclear Information System (INIS)

Dzhupanov, V.A.; Vassilev, V.M.; Dzhondzhorov, P.A.

1993-01-01

Plates and shells are basic structural components in nuclear reactors and their equipment. The prediction of the dynamic response of these components to fast transient loadings (e.g., loadings caused by earthquakes, missile impacts, etc.) is a quite important problem in the general context of the design, reliability and safety of nuclear power stations. Due to the extreme loading conditions a more adequate treatment of the foregoing problem should rest on a suitable nonlinear shell model, which would allow large deflections of the structures regarded to be taken into account. Such a model is provided in the nonlinear Donnell-Mushtari-Vlasov (DMV) theory. The governing system of equations of the DMV theory consists of two coupled nonlinear fourth order partial differential equations in three independent and two dependent variables. It is clear, as the case stands, that the obtaining solutions to this system directly, by using any of the general analytical or numerical techniques, would involve considerable difficulties. In the present paper, the invariance of the governing equations of DMV theory for plates and cylindrical shells relative to local Lie groups of local point transformations will be employed to get some advantages in connection with the aforementioned problem. First, the symmetry of a functional, corresponding to the governing equations of DMV theory for plates and cylindrical shells is studied. Next, the densities in the corresponding conservation laws are determined on the basis of Noether theorem. Finally, we study a class of invariant solutions of the governing equations. As is well known, group-invariant solutions are often intermediate asymptotics for a wider class of solutions of the corresponding equations. When such solutions are considered, the number of the independent variables can be reduced. For the class of invariant solutions studied here, the system of governing equations converts into a system of ordinary differential equations

Comparison between FEBio and Abaqus for biphasic contact problems.

Science.gov (United States)

Meng, Qingen; Jin, Zhongmin; Fisher, John; Wilcox, Ruth

2013-09-01

Articular cartilage plays an important role in the function of diarthrodial joints. Computational methods have been used to study the biphasic mechanics of cartilage, and Abaqus has been one of the most widely used commercial software packages for this purpose. A newly developed open-source finite element solver, FEBio, has been developed specifically for biomechanical applications. The aim of this study was to undertake a direct comparison between FEBio and Abaqus for some practical contact problems involving cartilage. Three model types, representing a porous flat-ended indentation test, a spherical-ended indentation test, and a conceptual natural joint contact model, were compared. In addition, a parameter sensitivity study was also performed for the spherical-ended indentation test to investigate the effects of changes in the input material properties on the model outputs, using both FEBio and Abaqus. Excellent agreement was found between FEBio and Abaqus for all of the model types and across the range of material properties that were investigated.

Existence results for boundary problems for uniformly elliptic and parabolic fully nonlinear equations

Directory of Open Access Journals (Sweden)

M. G. Crandall

1999-07-01

Full Text Available We study existence of continuous weak (viscosity solutions of Dirichlet and Cauchy-Dirichlet problems for fully nonlinear uniformly elliptic and parabolic equations. Two types of results are obtained in contexts where uniqueness of solutions fails or is unknown. For equations with merely measurable coefficients we prove solvability of the problem, while in the continuous case we construct maximal and minimal solutions. Necessary barriers on external cones are also constructed.

On the solvability of initial-value problems for nonlinear implicit difference equations

Directory of Open Access Journals (Sweden)

Ha Thi Ngoc Yen

2004-07-01

Full Text Available Our aim is twofold. First, we propose a natural definition of index for linear nonautonomous implicit difference equations, which is similar to that of linear differential-algebraic equations. Then we extend this index notion to a class of nonlinear implicit difference equations and prove some existence theorems for their initial-value problems.

Exact solution to the problem of nonlinear pulse propagation through random layered media and its connection with number triangles

International Nuclear Information System (INIS)

Sokolow, Adam; Sen, Surajit

2007-01-01

An energy pulse refers to a spatially compact energy bundle. In nonlinear pulse propagation, the nonlinearity of the relevant dynamical equations could lead to pulse propagation that is nondispersive or weakly dispersive in space and time. Nonlinear pulse propagation through layered media with widely varying pulse transmission properties is not wave-like and a problem of broad interest in many areas such as optics, geophysics, atmospheric physics and ocean sciences. We study nonlinear pulse propagation through a semi-infinite sequence of layers where the layers can have arbitrary energy transmission properties. By assuming that the layers are rigid, we are able to develop exact expressions for the backscattered energy received at the surface layer. The present study is likely to be relevant in the context of energy transport through soil and similar complex media. Our study reveals a surprising connection between the problem of pulse propagation and the number patterns in the well known Pascal's and Catalan's triangles and hence provides an analytic benchmark in a challenging problem of broad interest. We close with comments on the relationship between this study and the vast body of literature on the problem of wave localization in disordered systems

Admissible solutions for a class of nonlinear parabolic problem with non-negative data

Czech Academy of Sciences Publication Activity Database

Feireisl, Eduard; Petzeltová, Hana; Simondon, F.

2001-01-01

Roč. 131, č. 5 (2001), s. 857-883 ISSN 0308-2105 R&D Projects: GA AV ČR IAA1019703 Keywords : admissible solutions%nonlinear parabolic problem * admissible solutions * comparison principle * non-negative data Subject RIV: BA - General Mathematics Impact factor: 0.441, year: 2001

Shape Optimization in Three-Dimensional Contact Problems with Coulomb Friction

Czech Academy of Sciences Publication Activity Database

Beremlijski, P.; Haslinger, J.; Kočvara, Michal; Kučera, R.; Outrata, Jiří

2009-01-01

Roč. 20, č. 1 (2009), s. 416-444 ISSN 1052-6234 R&D Projects: GA AV ČR IAA100750802; GA AV ČR IAA1075402 Grant - others:European Commision(XE) FP6 - 30717; GA ČR(CZ) GA201/07/0294 Institutional research plan: CEZ:AV0Z10750506 Keywords : shape optimization * contact problems * Coulomb friction Subject RIV: BA - General Mathematics Impact factor: 1.429, year: 2009

Frictional contact problems for electro-viscoelastic materials with long-term memory, damage, and adhesion

Directory of Open Access Journals (Sweden)

Tedjani Hadj Ammar

2014-10-01

Full Text Available We consider a quasistatic contact problem between two electro-viscoelastic bodies with long-term memory and damage. The contact is frictional and is modelled with a version of normal compliance condition and the associated Coulomb's law of friction in which the adhesion of contact surfaces is taken into account. We derive a variational formulation for the model and prove an existence and uniqueness result of the weak solution. The proof is based on arguments of evolutionary variational inequalities, a classical existence and uniqueness result on parabolic inequalities, and Banach fixed point theorem.

On a solvability of hydro-mechanical problem based on contact problem with visco-plastic friction in Bingham rheology

Czech Academy of Sciences Publication Activity Database

Nedoma, Jiří; Tomášek, Luboš

2008-01-01

Roč. 218, č. 1 (2008), s. 116-124 ISSN 0377-0427 Institutional research plan: CEZ:AV0Z10300504 Keywords : visco-plasticity * Bingham rheology * contact problems with friction * variational inequalities * FEM * geomechanics * hydromechanics Subject RIV: BA - General Mathematics Impact factor: 1.048, year: 2008

Iterative Methods for Solving Nonlinear Parabolic Problem in Pension Saving Management

Science.gov (United States)

Koleva, M. N.

2011-11-01

In this work we consider a nonlinear parabolic equation, obtained from Riccati like transformation of the Hamilton-Jacobi-Bellman equation, arising in pension saving management. We discuss two numerical iterative methods for solving the model problem—fully implicit Picard method and mixed Picard-Newton method, which preserves the parabolic characteristics of the differential problem. Numerical experiments for comparison the accuracy and effectiveness of the algorithms are discussed. Finally, observations are given.

Numerical Modeling of Electrical Contact Conductance of Rough Bodies

Directory of Open Access Journals (Sweden)

M. V. Murashov

2015-01-01

Full Text Available Since the beginning of the 20th century to the present time, efforts have been made to develop a model of the electrical contact conductance. The development of micro- and nanotechnologies make contact conductance problem more essential. To conduct borrowing from a welldeveloped thermal contact conductance models on the basis of thermal and electrical conductivity analogy is often not possible due to a number of fundamental differences. While some 3Dmodels of rough bodies deformation have been developed in one way or another, a 3D-model of the electrical conductance through rough bodies contact is still not. A spatial model of electrical contact of rough bodies is proposed, allows one to calculate the electrical contact conductance as a function of the contact pressure. Representative elements of the bodies are parallelepipeds with deterministic roughness on the contacting surfaces. First the non-linear elastic-plastic deformation of rough surface under external pressure is solved using the finite element software ANSYS. Then the solution of electrostatic problem goes on the same finite element mesh. Aluminum AD1 is used as the material of the contacting bodies with properties that account for cold work hardening of the surface. The numerical model is built within the continuum mechanics and nanoscale effects are not taken into account. The electrical contact conductance was calculated on the basis of the concept of electrical resistance of the model as the sum of the electrical resistances of the contacting bodies and the contact itself. It was assumed that there is no air in the gap between the bodies. The dependence of the electrical contact conductance on the contact pressure is calculated as well as voltage and current density distributions in the contact bodies. It is determined that the multi-asperity contact mode, adequate to real roughness, is achieved at pressures higher than 3MPa, while results within the single contact spot are

TRUMP3-JR: a finite difference computer program for nonlinear heat conduction problems

International Nuclear Information System (INIS)

Ikushima, Takeshi

1984-02-01

Computer program TRUMP3-JR is a revised version of TRUMP3 which is a finite difference computer program used for the solution of multi-dimensional nonlinear heat conduction problems. Pre- and post-processings for input data generation and graphical representations of calculation results of TRUMP3 are avaiable in TRUMP3-JR. The calculation equations, program descriptions and user's instruction are presented. A sample problem is described to demonstrate the use of the program. (author)

Comparative Study of Evolutionary Multi-objective Optimization Algorithms for a Non-linear Greenhouse Climate Control Problem

DEFF Research Database (Denmark)

Ghoreishi, Newsha; Sørensen, Jan Corfixen; Jørgensen, Bo Nørregaard

2015-01-01

Non-trivial real world decision-making processes usually involve multiple parties having potentially conflicting interests over a set of issues. State-of-the-art multi-objective evolutionary algorithms (MOEA) are well known to solve this class of complex real-world problems. In this paper, we...... compare the performance of state-of-the-art multi-objective evolutionary algorithms to solve a non-linear multi-objective multi-issue optimisation problem found in Greenhouse climate control. The chosen algorithms in the study includes NSGAII, eNSGAII, eMOEA, PAES, PESAII and SPEAII. The performance...... of all aforementioned algorithms is assessed and compared using performance indicators to evaluate proximity, diversity and consistency. Our insights to this comparative study enhanced our understanding of MOEAs performance in order to solve a non-linear complex climate control problem. The empirical...

Nonlinear Dynamic Characteristics of the Railway Vehicle

Science.gov (United States)

Uyulan, Çağlar; Gokasan, Metin

2017-06-01

The nonlinear dynamic characteristics of a railway vehicle are checked into thoroughly by applying two different wheel-rail contact model: a heuristic nonlinear friction creepage model derived by using Kalker 's theory and Polach model including dead-zone clearance. This two models are matched with the quasi-static form of the LuGre model to obtain more realistic wheel-rail contact model. LuGre model parameters are determined using nonlinear optimization method, which it's objective is to minimize the error between the output of the Polach and Kalker model and quasi-static LuGre model for specific operating conditions. The symmetric/asymmetric bifurcation attitude and stable/unstable motion of the railway vehicle in the presence of nonlinearities which are yaw damping forces in the longitudinal suspension system are analyzed in great detail by changing the vehicle speed. Phase portraits of the lateral displacement of the leading wheelset of the railway vehicle are drawn below and on the critical speeds, where sub-critical Hopf bifurcation take place, for two wheel-rail contact model. Asymmetric periodic motions have been observed during the simulation in the lateral displacement of the wheelset under different vehicle speed range. The coexistence of multiple steady states cause bounces in the amplitude of vibrations, resulting instability problems of the railway vehicle. By using Lyapunov's indirect method, the critical hunting speeds are calculated with respect to the radius of the curved track parameter changes. Hunting, which is defined as the oscillation of the lateral displacement of wheelset with a large domain, is described by a limit cycle-type oscillation nature. The evaluated accuracy of the LuGre model adopted from Kalker's model results for prediction of critical speed is higher than the results of the LuGre model adopted from Polach's model. From the results of the analysis, the critical hunting speed must be resolved by investigating the track tests

On the the Contact Lens Problem: Modeling Rigid and Elastic Beams on Thin Films

Science.gov (United States)

Trinh, Philippe; Wilson, Stephen; Stone, Howard

2011-11-01

Generally, contact lenses are prescribed by the practitioner to fit each individual patient's eye, but these fitting-philosophies are based on empirical studies and a certain degree of trial-and-error. A badly fitted lens can cause a range of afflictions, which varies from mild dry-eye-discomfort, to more serious corneal diseases. Thus, at this heart of this problem, is the question of how a rigid or elastic plate interacts with the free-surface of a thin viscous film. In this talk, we present several mathematical models for the study of these plate-and-fluid problems. Asymptotic and numerical results are described, and we explain the role of elasticity, surface tension, viscosity, and pressure in determining the equilibrium solutions. Finally, we discuss the implications of our work on the contact lens problem, as well as on other coating processes which involve elastic substrates.

Fault detection for nonlinear systems - A standard problem approach

DEFF Research Database (Denmark)

Stoustrup, Jakob; Niemann, Hans Henrik

1998-01-01

The paper describes a general method for designing (nonlinear) fault detection and isolation (FDI) systems for nonlinear processes. For a rich class of nonlinear systems, a nonlinear FDI system can be designed using convex optimization procedures. The proposed method is a natural extension...

Nonlinear Eigenvalue Problems in Elliptic Variational Inequalities: a local study

International Nuclear Information System (INIS)

Conrad, F.; Brauner, C.; Issard-Roch, F.; Nicolaenko, B.

1985-01-01

The authors consider a class of Nonlinear Eigenvalue Problems (N.L.E.P.) associated with Elliptic Variational Inequalities (E.V.I.). First the authors introduce the main tools for a local study of branches of solutions; the authors extend the linearization process required in the case of equations. Next the authors prove the existence of arcs of solutions close to regular vs singular points, and determine their local behavior up to the first order. Finally, the authors discuss the connection between their regularity condition and some stability concept. 37 references, 6 figures

Numerical Simulation of Some Biomechanical Problems

Czech Academy of Sciences Publication Activity Database

Nedoma, Jiří; Klézl, Z.; Fousek, J.; Kestřánek, Zdeněk; Stehlík, J.

2003-01-01

Roč. 61, 3-6 (2003), s. 283-295 ISSN 0378-4754. [MODELLING 2001. IMACS Conference on Mathematical Modelling and Computational Methods in Mechanics, Physics , Biomechanics and Geodynamics /2./. Pilsen, 19.06.2001-25.06.2001] Institutional research plan: AV0Z1030915 Keywords : non-linear elasticity * contact problems * variational inequality * finite element method * wrist * spine * fracture Subject RIV: BA - General Mathematics Impact factor: 0.558, year: 2003

DISPL-1, 2. Order Nonlinear Partial Differential Equation System Solution for Kinetics Diffusion Problems

International Nuclear Information System (INIS)

Leaf, G.K.; Minkoff, M.

1982-01-01

1 - Description of problem or function: DISPL1 is a software package for solving second-order nonlinear systems of partial differential equations including parabolic, elliptic, hyperbolic, and some mixed types. The package is designed primarily for chemical kinetics- diffusion problems, although not limited to these problems. Fairly general nonlinear boundary conditions are allowed as well as inter- face conditions for problems in an inhomogeneous medium. The spatial domain is one- or two-dimensional with rectangular Cartesian, cylindrical, or spherical (in one dimension only) geometry. 2 - Method of solution: The numerical method is based on the use of Galerkin's procedure combined with the use of B-Splines (C.W.R. de-Boor's B-spline package) to generate a system of ordinary differential equations. These equations are solved by a sophisticated ODE software package which is a modified version of Hindmarsh's GEAR package, NESC Abstract 592. 3 - Restrictions on the complexity of the problem: The spatial domain must be rectangular with sides parallel to the coordinate geometry. Cross derivative terms are not permitted in the PDE. The order of the B-Splines is at most 12. Other parameters such as the number of mesh points in each coordinate direction, the number of PDE's etc. are set in a macro table used by the MORTRAn2 preprocessor in generating the object code

New implementation method for essential boundary condition to extended element-free Galerkin method. Application to nonlinear problem

International Nuclear Information System (INIS)

Saitoh, Ayumu; Matsui, Nobuyuki; Itoh, Taku; Kamitani, Atsushi; Nakamura, Hiroaki

2011-01-01

A new method has been proposed for implementing essential boundary conditions to the Element-Free Galerkin Method (EFGM) without using the Lagrange multiplier. Furthermore, the performance of the proposed method has been investigated for a nonlinear Poisson problem. The results of computations show that, as interpolation functions become closer to delta functions, the accuracy of the solution is improved on the boundary. In addition, the accuracy of the proposed method is higher than that of the conventional EFGM. Therefore, it might be concluded that the proposed method is useful for solving the nonlinear Poisson problem. (author)

Kalker's algorithm Fastsim solves tangential contact problems with slip-dependent friction and friction anisotropy

Science.gov (United States)

Piotrowski, J.

2010-07-01

This paper presents two extensions of Kalker's algorithm Fastsim of the simplified theory of rolling contact. The first extension is for solving tangential contact problems with the coefficient of friction depending on slip velocity. Two friction laws have been considered: with and without recuperation of the static friction. According to the tribological hypothesis for metallic bodies shear failure, the friction law without recuperation of static friction is more suitable for wheel and rail than the other one. Sample results present local quantities inside the contact area (division to slip and adhesion, traction) as well as global ones (creep forces as functions of creepages and rolling velocity). For the coefficient of friction diminishing with slip, the creep forces decay after reaching the maximum and they depend on the rolling velocity. The second extension is for solving tangential contact problems with friction anisotropy characterised by a convex set of the permissible tangential tractions. The effect of the anisotropy has been shown on examples of rolling without spin and in the presence of pure spin for the elliptical set. The friction anisotropy influences tangential tractions and creep forces. Sample results present local and global quantities. Both extensions have been described with the same language of formulation and they may be merged into one, joint algorithm.

Semi-analog Monte Carlo (SMC) method for time-dependent non-linear three-dimensional heterogeneous radiative transfer problems

International Nuclear Information System (INIS)

Yun, Sung Hwan

2004-02-01

Radiative transfer is a complex phenomenon in which radiation field interacts with material. This thermal radiative transfer phenomenon is composed of two equations which are the balance equation of photons and the material energy balance equation. The two equations involve non-linearity due to the temperature and that makes the radiative transfer equation more difficult to solve. During the last several years, there have been many efforts to solve the non-linear radiative transfer problems by Monte Carlo method. Among them, it is known that Semi-Analog Monte Carlo (SMC) method developed by Ahrens and Larsen is accurate regard-less of the time step size in low temperature region. But their works are limited to one-dimensional, low temperature problems. In this thesis, we suggest some method to remove their limitations in the SMC method and apply to the more realistic problems. An initially cold problem was solved over entire temperature region by using piecewise linear interpolation of the heat capacity, while heat capacity is still fitted as a cubic curve within the lowest temperature region. If we assume the heat capacity to be linear in each temperature region, the non-linearity still remains in the radiative transfer equations. We then introduce the first-order Taylor expansion to linearize the non-linear radiative transfer equations. During the linearization procedure, absorption-reemission phenomena may be described by a conventional reemission time sampling scheme which is similar to the repetitive sampling scheme in particle transport Monte Carlo method. But this scheme causes significant stochastic errors, which necessitates many histories. Thus, we present a new reemission time sampling scheme which reduces stochastic errors by storing the information of absorption times. The results of the comparison of the two schemes show that the new scheme has less stochastic errors. Therefore, the improved SMC method is able to solve more realistic problems with

Unique solvability of a non-linear non-local boundary-value problem for systems of non-linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Dilna, N.; Rontó, András

2010-01-01

Roč. 60, č. 3 (2010), s. 327-338 ISSN 0139-9918 R&D Projects: GA ČR(CZ) GA201/06/0254 Institutional research plan: CEZ:AV0Z10190503 Keywords : non-linear boundary value-problem * functional differential equation * non-local condition * unique solvability * differential inequality Subject RIV: BA - General Mathematics Impact factor: 0.316, year: 2010 http://link.springer.com/article/10.2478%2Fs12175-010-0015-9

On Nonlinear Inverse Problems of Heat Transfer with Radiation Boundary Conditions: Application to Dehydration of Gypsum Plasterboards Exposed to Fire

OpenAIRE

Belmiloudi, A.; Mahé, F.

2014-01-01

International audience; The paper investigates boundary optimal controls and parameter estimates to the well-posedness nonlinear model of dehydration of thermic problems. We summarize the general formulations for the boundary control for initial-boundary value problem for nonlinear partial differential equations modeling the heat transfer and derive necessary optimality conditions, including the adjoint equation, for the optimal set of parameters minimizing objective functions J. Numerical si...

Nonlinear time-series analysis of current signal in cathodic contact glow discharge electrolysis

International Nuclear Information System (INIS)

Allagui, Anis; Abdelkareem, Mohammad Ali; Rojas, Andrea Espinel; Bonny, Talal; Elwakil, Ahmed S.

2016-01-01

In the standard two-electrode configuration employed in electrolytic process, when the control dc voltage is brought to a critical value, the system undergoes a transition from conventional electrolysis to contact glow discharge electrolysis (CGDE), which has also been referred to as liquid-submerged micro-plasma, glow discharge plasma electrolysis, electrode effect, electrolytic plasma, etc. The light-emitting process is associated with the development of an irregular and erratic current time-series which has been arbitrarily labelled as “random,” and thus dissuaded further research in this direction. Here, we examine the current time-series signals measured in cathodic CGDE configuration in a concentrated KOH solution at different dc bias voltages greater than the critical voltage. We show that the signals are, in fact, not random according to the NIST SP. 800-22 test suite definition. We also demonstrate that post-processing low-pass filtered sequences requires less time than the native as-measured sequences, suggesting a superposition of low frequency chaotic fluctuations and high frequency behaviors (which may be produced by more than one possible source of entropy). Using an array of nonlinear time-series analyses for dynamical systems, i.e., the computation of largest Lyapunov exponents and correlation dimensions, and re-construction of phase portraits, we found that low-pass filtered datasets undergo a transition from quasi-periodic to chaotic to quasi-hyper-chaotic behavior, and back again to chaos when the voltage controlling-parameter is increased. The high frequency part of the signals is discussed in terms of highly nonlinear turbulent motion developed around the working electrode.

A Semismooth Newton Method for Nonlinear Parameter Identification Problems with Impulsive Noise

KAUST Repository

Clason, Christian

2012-01-01

This work is concerned with nonlinear parameter identification in partial differential equations subject to impulsive noise. To cope with the non-Gaussian nature of the noise, we consider a model with L 1 fitting. However, the nonsmoothness of the problem makes its efficient numerical solution challenging. By approximating this problem using a family of smoothed functionals, a semismooth Newton method becomes applicable. In particular, its superlinear convergence is proved under a second-order condition. The convergence of the solution to the approximating problem as the smoothing parameter goes to zero is shown. A strategy for adaptively selecting the regularization parameter based on a balancing principle is suggested. The efficiency of the method is illustrated on several benchmark inverse problems of recovering coefficients in elliptic differential equations, for which one- and two-dimensional numerical examples are presented. © by SIAM.

Nonlinear problems of the theory of heterogeneous slightly curved shells

Science.gov (United States)

Kantor, B. Y.

1973-01-01

An account if given of the variational method of the solution of physically and geometrically nonlinear problems of the theory of heterogeneous slightly curved shells. Examined are the bending and supercritical behavior of plates and conical and spherical cupolas of variable thickness in a temperature field, taking into account the dependence of the elastic parameters on temperature. The bending, stability in general and load-bearing capacity of flexible isotropic elastic-plastic shells with different criteria of plasticity, taking into account compressibility and hardening. The effect of the plastic heterogeneity caused by heat treatment, surface work hardening and irradiation by fast neutron flux is investigated. Some problems of the dynamic behavior of flexible shells are solved. Calculations are performed in high approximations. Considerable attention is given to the construction of a machine algorithm and to the checking of the convergence of iterative processes.

LDRD report nonlinear model reduction

Energy Technology Data Exchange (ETDEWEB)

Segalman, D.; Heinstein, M.

1997-09-01

The very general problem of model reduction of nonlinear systems was made tractable by focusing on the very large subclass consisting of linear subsystems connected by nonlinear interfaces. Such problems constitute a large part of the nonlinear structural problems encountered in addressing the Sandia missions. A synthesis approach to this class of problems was developed consisting of: detailed modeling of the interface mechanics; collapsing the interface simulation results into simple nonlinear interface models; constructing system models by assembling model approximations of the linear subsystems and the nonlinear interface models. These system models, though nonlinear, would have very few degrees of freedom. A paradigm problem, that of machine tool vibration, was selected for application of the reduction approach outlined above. Research results achieved along the way as well as the overall modeling of a specific machine tool have been very encouraging. In order to confirm the interface models resulting from simulation, it was necessary to develop techniques to deduce interface mechanics from experimental data collected from the overall nonlinear structure. A program to develop such techniques was also pursued with good success.

An artificial nonlinear diffusivity method for supersonic reacting flows with shocks

Science.gov (United States)

Fiorina, B.; Lele, S. K.

2007-03-01

A computational approach for modeling interactions between shocks waves, contact discontinuities and reactions zones with a high-order compact scheme is investigated. To prevent the formation of spurious oscillations around shocks, artificial nonlinear viscosity [A.W. Cook, W.H. Cabot, A high-wavenumber viscosity for high resolution numerical method, J. Comput. Phys. 195 (2004) 594-601] based on high-order derivative of the strain rate tensor is used. To capture temperature and species discontinuities a nonlinear diffusivity based on the entropy gradient is added. It is shown that the damping of 'wiggles' is controlled by the model constants and is largely independent of the mesh size and the shock strength. The same holds for the numerical shock thickness and allows a determination of the L2 error. In the shock tube problem, with fluids of different initial entropy separated by the diaphragm, an artificial diffusivity is required to accurately capture the contact surface. Finally, the method is applied to a shock wave propagating into a medium with non-uniform density/entropy and to a CJ detonation wave. Multi-dimensional formulation of the model is presented and is illustrated by a 2D oblique wave reflection from an inviscid wall, by a 2D supersonic blunt body flow and by a Mach reflection problem.

POSITIVE SOLUTIONS OF A NONLINEAR THREE-POINT EIGENVALUE PROBLEM WITH INTEGRAL BOUNDARY CONDITIONS

Directory of Open Access Journals (Sweden)

FAOUZI HADDOUCHI

2015-11-01

Full Text Available In this paper, we study the existence of positive solutions of a three-point integral boundary value problem (BVP for the following second-order differential equation u''(t + \\lambda a(tf(u(t = 0; 0 0 is a parameter, 0 <\\eta < 1, 0 <\\alpha < 1/{\\eta}. . By using the properties of the Green's function and Krasnoselskii's fixed point theorem on cones, the eigenvalue intervals of the nonlinear boundary value problem are considered, some sufficient conditions for the existence of at least one positive solutions are established.

Efficient Output Solution for Nonlinear Stochastic Optimal Control Problem with Model-Reality Differences

Directory of Open Access Journals (Sweden)

Sie Long Kek

2015-01-01

Full Text Available A computational approach is proposed for solving the discrete time nonlinear stochastic optimal control problem. Our aim is to obtain the optimal output solution of the original optimal control problem through solving the simplified model-based optimal control problem iteratively. In our approach, the adjusted parameters are introduced into the model used such that the differences between the real system and the model used can be computed. Particularly, system optimization and parameter estimation are integrated interactively. On the other hand, the output is measured from the real plant and is fed back into the parameter estimation problem to establish a matching scheme. During the calculation procedure, the iterative solution is updated in order to approximate the true optimal solution of the original optimal control problem despite model-reality differences. For illustration, a wastewater treatment problem is studied and the results show the efficiency of the approach proposed.

A parallel multi-domain solution methodology applied to nonlinear thermal transport problems in nuclear fuel pins

Energy Technology Data Exchange (ETDEWEB)

Philip, Bobby, E-mail: philipb@ornl.gov [Oak Ridge National Laboratory, One Bethel Valley Road, Oak Ridge, TN 37831 (United States); Berrill, Mark A.; Allu, Srikanth; Hamilton, Steven P.; Sampath, Rahul S.; Clarno, Kevin T. [Oak Ridge National Laboratory, One Bethel Valley Road, Oak Ridge, TN 37831 (United States); Dilts, Gary A. [Los Alamos National Laboratory, PO Box 1663, Los Alamos, NM 87545 (United States)

2015-04-01

This paper describes an efficient and nonlinearly consistent parallel solution methodology for solving coupled nonlinear thermal transport problems that occur in nuclear reactor applications over hundreds of individual 3D physical subdomains. Efficiency is obtained by leveraging knowledge of the physical domains, the physics on individual domains, and the couplings between them for preconditioning within a Jacobian Free Newton Krylov method. Details of the computational infrastructure that enabled this work, namely the open source Advanced Multi-Physics (AMP) package developed by the authors is described. Details of verification and validation experiments, and parallel performance analysis in weak and strong scaling studies demonstrating the achieved efficiency of the algorithm are presented. Furthermore, numerical experiments demonstrate that the preconditioner developed is independent of the number of fuel subdomains in a fuel rod, which is particularly important when simulating different types of fuel rods. Finally, we demonstrate the power of the coupling methodology by considering problems with couplings between surface and volume physics and coupling of nonlinear thermal transport in fuel rods to an external radiation transport code.

A parallel algorithm for transient solid dynamics simulations with contact detection

International Nuclear Information System (INIS)

Attaway, S.; Hendrickson, B.; Plimpton, S.; Gardner, D.; Vaughan, C.; Heinstein, M.; Peery, J.

1996-01-01

Solid dynamics simulations with Lagrangian finite elements are used to model a wide variety of problems, such as the calculation of impact damage to shipping containers for nuclear waste and the analysis of vehicular crashes. Using parallel computers for these simulations has been hindered by the difficulty of searching efficiently for material surface contacts in parallel. A new parallel algorithm for calculation of arbitrary material contacts in finite element simulations has been developed and implemented in the PRONTO3D transient solid dynamics code. This paper will explore some of the issues involved in developing efficient, portable, parallel finite element models for nonlinear transient solid dynamics simulations. The contact-detection problem poses interesting challenges for efficient implementation of a solid dynamics simulation on a parallel computer. The finite element mesh is typically partitioned so that each processor owns a localized region of the finite element mesh. This mesh partitioning is optimal for the finite element portion of the calculation since each processor must communicate only with the few connected neighboring processors that share boundaries with the decomposed mesh. However, contacts can occur between surfaces that may be owned by any two arbitrary processors. Hence, a global search across all processors is required at every time step to search for these contacts. Load-imbalance can become a problem since the finite element decomposition divides the volumetric mesh evenly across processors but typically leaves the surface elements unevenly distributed. In practice, these complications have been limiting factors in the performance and scalability of transient solid dynamics on massively parallel computers. In this paper the authors present a new parallel algorithm for contact detection that overcomes many of these limitations

Positive solutions for a nonlinear periodic boundary-value problem with a parameter

Directory of Open Access Journals (Sweden)

Jingliang Qiu

2012-08-01

Full Text Available Using topological degree theory with a partially ordered structure of space, sufficient conditions for the existence and multiplicity of positive solutions for a second-order nonlinear periodic boundary-value problem are established. Inspired by ideas in Guo and Lakshmikantham [6], we study the dependence of positive periodic solutions as a parameter approaches infinity, $$ lim_{lambdao +infty}|x_{lambda}|=+infty,quadhbox{or}quad lim_{lambdao+infty}|x_{lambda}|=0. $$

Numerical solution of large nonlinear boundary value problems by quadratic minimization techniques

International Nuclear Information System (INIS)

Glowinski, R.; Le Tallec, P.

1984-01-01

The objective of this paper is to describe the numerical treatment of large highly nonlinear two or three dimensional boundary value problems by quadratic minimization techniques. In all the different situations where these techniques were applied, the methodology remains the same and is organized as follows: 1) derive a variational formulation of the original boundary value problem, and approximate it by Galerkin methods; 2) transform this variational formulation into a quadratic minimization problem (least squares methods) or into a sequence of quadratic minimization problems (augmented lagrangian decomposition); 3) solve each quadratic minimization problem by a conjugate gradient method with preconditioning, the preconditioning matrix being sparse, positive definite, and fixed once for all in the iterative process. This paper will illustrate the methodology above on two different examples: the description of least squares solution methods and their application to the solution of the unsteady Navier-Stokes equations for incompressible viscous fluids; the description of augmented lagrangian decomposition techniques and their application to the solution of equilibrium problems in finite elasticity

Eigenvalue problems for degenerate nonlinear elliptic equations in anisotropic media

Directory of Open Access Journals (Sweden)

Vicenţiu RăDulescu

2005-06-01

Full Text Available We study nonlinear eigenvalue problems of the type −div(a(x∇u=g(λ,x,u in ℝN, where a(x is a degenerate nonnegative weight. We establish the existence of solutions and we obtain information on qualitative properties as multiplicity and location of solutions. Our approach is based on the critical point theory in Sobolev weighted spaces combined with a Caffarelli-Kohn-Nirenberg-type inequality. A specific minimax method is developed without making use of Palais-Smale condition.

Nonlinear behavior of three-terminal graphene junctions at room temperature

International Nuclear Information System (INIS)

Kim, Wonjae; Riikonen, Juha; Lipsanen, Harri; Pasanen, Pirjo

2012-01-01

We demonstrate nonlinear behavior in three-terminal T-branch graphene devices at room temperature. A rectified nonlinear output at the center branch is observed when the device is biased by a push–pull configuration. Nonlinearity is assumed to arise from a difference in charge transfer through the metal–graphene contact barrier between two contacts. The sign of the rectification can be altered by changing the carrier type using the back-gate voltage. (paper)

Nonlinear approaches in engineering applications advanced analysis of vehicle related technologies

CERN Document Server

Dai, Liming

2016-01-01

This book looks at the broad field of engineering science through the lens of nonlinear approaches. Examples focus on issues in vehicle technology, including vehicle dynamics, vehicle-road interaction, steering, and control for electric and hybrid vehicles. Also included are discussions on train and tram systems, aerial vehicles, robot-human interaction, and contact and scratch analysis at the micro/nanoscale. Chapters are based on invited contributions from world-class experts in the field who advance the future of engineering by discussing the development of more optimal, accurate, efficient, and cost and energy effective systems. This book is appropriate for researchers, students, and practicing engineers who are interested in the applications of nonlinear approaches to solving engineering and science problems.

JAC, 2-D Finite Element Method Program for Quasi Static Mechanics Problems by Nonlinear Conjugate Gradient (CG) Method

International Nuclear Information System (INIS)

Biffle, J.H.

1991-01-01

1 - Description of program or function: JAC is a two-dimensional finite element program for solving large deformation, temperature dependent, quasi-static mechanics problems with the nonlinear conjugate gradient (CG) technique. Either plane strain or axisymmetric geometry may be used with material descriptions which include temperature dependent elastic-plastic, temperature dependent secondary creep, and isothermal soil models. The nonlinear effects examined include material and geometric nonlinearities due to large rotations, large strains, and surface which slide relative to one another. JAC is vectorized to perform efficiently on the Cray1 computer. A restart capability is included. 2 - Method of solution: The nonlinear conjugate gradient method is employed in a two-dimensional plane strain or axisymmetric setting with various techniques for accelerating convergence. Sliding interface conditions are also implemented. A four-node Lagrangian uniform strain element is used with orthogonal hourglass viscosity to control the zero energy modes. Three sets of continuum equations are needed - kinematic statements, constitutive equations, and equations of equilibrium - to describe the deformed configuration of the body. 3 - Restrictions on the complexity of the problem - Maxima of: 10 load and solution control functions, 4 materials. The strain rate is assumed constant over a time interval. Current large rotation theory is applicable to a maximum shear strain of 1.0. JAC should be used with caution for large shear strains. Problem size is limited only by available memory

A comparison of two efficient nonlinear heat conduction methodologies using a two-dimensional time-dependent benchmark problem

International Nuclear Information System (INIS)

Wilson, G.L.; Rydin, R.A.; Orivuori, S.

1988-01-01

Two highly efficient nonlinear time-dependent heat conduction methodologies, the nonlinear time-dependent nodal integral technique (NTDNT) and IVOHEAT are compared using one- and two-dimensional time-dependent benchmark problems. The NTDNT is completely based on newly developed time-dependent nodal integral methods, whereas IVOHEAT is based on finite elements in space and Crank-Nicholson finite differences in time. IVOHEAT contains the geometric flexibility of the finite element approach, whereas the nodal integral method is constrained at present to Cartesian geometry. For test problems where both methods are equally applicable, the nodal integral method is approximately six times more efficient per dimension than IVOHEAT when a comparable overall accuracy is chosen. This translates to a factor of 200 for a three-dimensional problem having relatively homogeneous regions, and to a smaller advantage as the degree of heterogeneity increases

On the solvability of asymmetric quasilinear finite element approximate problems in nonlinear incompressible elasticity

International Nuclear Information System (INIS)

Ruas, V.

1982-09-01

A class of simplicial finite elements for solving incompressible elasticity problems in n-dimensional space, n=2 or 3, is presented. An asymmetric structure of the shape functions with respect to the centroid of the simplex, renders them particularly stable in the large strain case, in which the incompressibility condition is nonlinear. It is proved that under certain assembling conditions of the elements, there exists a solution to the corresponding discrete problems. Numerical examples illustrate the efficiency of the method. (Author) [pt

Solution of the nonlinear inverse scattering problem by T-matrix completion. I. Theory.

Science.gov (United States)

Levinson, Howard W; Markel, Vadim A

2016-10-01

We propose a conceptually different method for solving nonlinear inverse scattering problems (ISPs) such as are commonly encountered in tomographic ultrasound imaging, seismology, and other applications. The method is inspired by the theory of nonlocality of physical interactions and utilizes the relevant formalism. We formulate the ISP as a problem whose goal is to determine an unknown interaction potential V from external scattering data. Although we seek a local (diagonally dominated) V as the solution to the posed problem, we allow V to be nonlocal at the intermediate stages of iterations. This allows us to utilize the one-to-one correspondence between V and the T matrix of the problem. Here it is important to realize that not every T corresponds to a diagonal V and we, therefore, relax the usual condition of strict diagonality (locality) of V. An iterative algorithm is proposed in which we seek T that is (i) compatible with the measured scattering data and (ii) corresponds to an interaction potential V that is as diagonally dominated as possible. We refer to this algorithm as to the data-compatible T-matrix completion. This paper is Part I in a two-part series and contains theory only. Numerical examples of image reconstruction in a strongly nonlinear regime are given in Part II [H. W. Levinson and V. A. Markel, Phys. Rev. E 94, 043318 (2016)10.1103/PhysRevE.94.043318]. The method described in this paper is particularly well suited for very large data sets that become increasingly available with the use of modern measurement techniques and instrumentation.

Near-optimal alternative generation using modified hit-and-run sampling for non-linear, non-convex problems

Science.gov (United States)

Rosenberg, D. E.; Alafifi, A.

2016-12-01

Water resources systems analysis often focuses on finding optimal solutions. Yet an optimal solution is optimal only for the modelled issues and managers often seek near-optimal alternatives that address un-modelled objectives, preferences, limits, uncertainties, and other issues. Early on, Modelling to Generate Alternatives (MGA) formalized near-optimal as the region comprising the original problem constraints plus a new constraint that allowed performance within a specified tolerance of the optimal objective function value. MGA identified a few maximally-different alternatives from the near-optimal region. Subsequent work applied Markov Chain Monte Carlo (MCMC) sampling to generate a larger number of alternatives that span the near-optimal region of linear problems or select portions for non-linear problems. We extend the MCMC Hit-And-Run method to generate alternatives that span the full extent of the near-optimal region for non-linear, non-convex problems. First, start at a feasible hit point within the near-optimal region, then run a random distance in a random direction to a new hit point. Next, repeat until generating the desired number of alternatives. The key step at each iterate is to run a random distance along the line in the specified direction to a new hit point. If linear equity constraints exist, we construct an orthogonal basis and use a null space transformation to confine hits and runs to a lower-dimensional space. Linear inequity constraints define the convex bounds on the line that runs through the current hit point in the specified direction. We then use slice sampling to identify a new hit point along the line within bounds defined by the non-linear inequity constraints. This technique is computationally efficient compared to prior near-optimal alternative generation techniques such MGA, MCMC Metropolis-Hastings, evolutionary, or firefly algorithms because search at each iteration is confined to the hit line, the algorithm can move in one

A finite-element model for moving contact line problems in immiscible two-phase flow

Science.gov (United States)

Kucala, Alec

2017-11-01

Accurate modeling of moving contact line (MCL) problems is imperative in predicting capillary pressure vs. saturation curves, permeability, and preferential flow paths for a variety of applications, including geological carbon storage (GCS) and enhanced oil recovery (EOR). The macroscale movement of the contact line is dependent on the molecular interactions occurring at the three-phase interface, however most MCL problems require resolution at the meso- and macro-scale. A phenomenological model must be developed to account for the microscale interactions, as resolving both the macro- and micro-scale would render most problems computationally intractable. Here, a model for the moving contact line is presented as a weak forcing term in the Navier-Stokes equation and applied directly at the location of the three-phase interface point. The moving interface is tracked with the level set method and discretized using the conformal decomposition finite element method (CDFEM), allowing for the surface tension and the wetting model to be computed at the exact interface location. A variety of verification test cases for simple two- and three-dimensional geometries are presented to validate the current MCL model, which can exhibit grid independence when a proper scaling for the slip length is chosen. Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA-0003525.

Nonlinear wave equations

CERN Document Server

Li, Tatsien

2017-01-01

This book focuses on nonlinear wave equations, which are of considerable significance from both physical and theoretical perspectives. It also presents complete results on the lower bound estimates of lifespan (including the global existence), which are established for classical solutions to the Cauchy problem of nonlinear wave equations with small initial data in all possible space dimensions and with all possible integer powers of nonlinear terms. Further, the book proposes the global iteration method, which offers a unified and straightforward approach for treating these kinds of problems. Purely based on the properties of solut ions to the corresponding linear problems, the method simply applies the contraction mapping principle.

A Hartman–Nagumo inequality for the vector ordinary -Laplacian and applications to nonlinear boundary value problems

Directory of Open Access Journals (Sweden)

Ureña Antonio J

2002-01-01

Full Text Available A generalization of the well-known Hartman–Nagumo inequality to the case of the vector ordinary -Laplacian and classical degree theory provide existence results for some associated nonlinear boundary value problems.

Some recent advances in 3D crack and contact analysis of elastic solids with transverse isotropy and multifield coupling

Science.gov (United States)

Chen, Wei-Qiu

2015-10-01

Significant progress has been made in mixed boundary-value problems associated with three-dimensional (3D) crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials. These include material anisotropy and multifield coupling, two typical characteristics of most current multifunctional materials. In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V. I. Fabrikant, whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory. We are particularly interested in crack and contact problems with certain nonlinear features. Emphasis is also placed on the coupling between the temperature field (or the like) and other physical fields (e.g., elastic, electric, and magnetic fields). We further highlight the practical significance of 3D contact solutions, in particular in applications related to modern scanning probe microscopes.

Nonlinear Preconditioning and its Application in Multicomponent Problems

KAUST Repository

Liu, Lulu

2015-01-01

the convergence of systems with unbalanced nonlinearities; however, they have natural complementarity in practice. MSPIN is naturally based on partitioning of degrees of freedom in a nonlinear PDE system by field type rather than by subdomain, where a modest

Existence for viscoplastic contact with Coulomb friction problems

Directory of Open Access Journals (Sweden)

Amina Amassad

2002-01-01

frictional contact between an elastic-viscoplastic body and a rigid obstacle. We model the frictional contact both by a Tresca's friction law and a regularized Coulomb's law. We assume, in a first part, that the contact is bilateral and that no separation takes place. In a second part, we consider the Signorini unilateral contact conditions. Proofs are based on a time-discretization method, Banach and Schauder fixed point theorems.

Surface tension and contact angles: Molecular origins and associated microstructure

Science.gov (United States)

Davis, H. T.

1982-01-01

Gradient theory converts the molecular theory of inhomogeneous fluid into nonlinear boundary value problems for density and stress distributions in fluid interfaces, contact line regions, nuclei and microdroplets, and other fluid microstructures. The relationship between the basic patterns of fluid phase behavior and the occurrence and stability of fluid microstructures was clearly established by the theory. All the inputs of the theory have molecular expressions which are computable from simple models. On another level, the theory becomes a phenomenological framework in which the equation of state of homogeneous fluid and sets of influence parameters of inhomogeneous fluids are the inputs and the structures, stress tensions and contact angles of menisci are the outputs. These outputs, which find applications in the science and technology of drops and bubbles, are discussed.

Ensemble Kalman Filtering with Residual Nudging: An Extension to State Estimation Problems with Nonlinear Observation Operators

KAUST Repository

Luo, Xiaodong

2014-10-01

The ensemble Kalman filter (EnKF) is an efficient algorithm for many data assimilation problems. In certain circumstances, however, divergence of the EnKF might be spotted. In previous studies, the authors proposed an observation-space-based strategy, called residual nudging, to improve the stability of the EnKF when dealing with linear observation operators. The main idea behind residual nudging is to monitor and, if necessary, adjust the distances (misfits) between the real observations and the simulated ones of the state estimates, in the hope that by doing so one may be able to obtain better estimation accuracy. In the present study, residual nudging is extended and modified in order to handle nonlinear observation operators. Such extension and modification result in an iterative filtering framework that, under suitable conditions, is able to achieve the objective of residual nudging for data assimilation problems with nonlinear observation operators. The 40-dimensional Lorenz-96 model is used to illustrate the performance of the iterative filter. Numerical results show that, while a normal EnKF may diverge with nonlinear observation operators, the proposed iterative filter remains stable and leads to reasonable estimation accuracy under various experimental settings.

On the Cauchy problem for a Sobolev-type equation with quadratic non-linearity

International Nuclear Information System (INIS)

Aristov, Anatoly I

2011-01-01

We investigate the asymptotic behaviour as t→∞ of the solution of the Cauchy problem for a Sobolev-type equation with quadratic non-linearity and develop ideas used by I. A. Shishmarev and other authors in the study of classical and Sobolev-type equations. Conditions are found under which it is possible to consider the case of an arbitrary dimension of the spatial variable.

Positive solutions of nonlinear fractional boundary value problems with Dirichlet boundary conditions

Directory of Open Access Journals (Sweden)

Qingkai Kong

2012-02-01

Full Text Available In this paper, we study the existence and multiplicity of positive solutions of a class of nonlinear fractional boundary value problems with  Dirichlet boundary conditions. By applying the fixed point theory on cones we establish a series of criteria for the existence of one, two, any arbitrary finite number, and an infinite number of positive solutions. A criterion for the nonexistence of positive solutions is also derived. Several examples are given for demonstration.

Indefinitely preconditioned inexact Newton method for large sparse equality constrained non-linear programming problems

Czech Academy of Sciences Publication Activity Database

Lukšan, Ladislav; Vlček, Jan

1998-01-01

Roč. 5, č. 3 (1998), s. 219-247 ISSN 1070-5325 R&D Projects: GA ČR GA201/96/0918 Keywords : nonlinear programming * sparse problems * equality constraints * truncated Newton method * augmented Lagrangian function * indefinite systems * indefinite preconditioners * conjugate gradient method * residual smoothing Subject RIV: BA - General Mathematics Impact factor: 0.741, year: 1998

The focal boundary value problem for strongly singular higher-order nonlinear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Půža, B.

2015-01-01

Roč. 2015, January (2015), s. 17 ISSN 1687-2770 Institutional support: RVO:67985840 Keywords : higher order nonlinear functional-differential equations * two-point right-focal boundary value problem * strong singularity Subject RIV: BA - General Mathematics Impact factor: 0.642, year: 2015 http://link.springer.com/article/10.1186%2Fs13661-014-0277-1

An iterative kernel based method for fourth order nonlinear equation with nonlinear boundary condition

Science.gov (United States)

Azarnavid, Babak; Parand, Kourosh; Abbasbandy, Saeid

2018-06-01

This article discusses an iterative reproducing kernel method with respect to its effectiveness and capability of solving a fourth-order boundary value problem with nonlinear boundary conditions modeling beams on elastic foundations. Since there is no method of obtaining reproducing kernel which satisfies nonlinear boundary conditions, the standard reproducing kernel methods cannot be used directly to solve boundary value problems with nonlinear boundary conditions as there is no knowledge about the existence and uniqueness of the solution. The aim of this paper is, therefore, to construct an iterative method by the use of a combination of reproducing kernel Hilbert space method and a shooting-like technique to solve the mentioned problems. Error estimation for reproducing kernel Hilbert space methods for nonlinear boundary value problems have yet to be discussed in the literature. In this paper, we present error estimation for the reproducing kernel method to solve nonlinear boundary value problems probably for the first time. Some numerical results are given out to demonstrate the applicability of the method.

Semi-Smooth Newton Method for Solving 2D Contact Problems with Tresca and Coulomb Friction

Directory of Open Access Journals (Sweden)

Kristina Motyckova

2013-01-01

Full Text Available The contribution deals with contact problems for two elastic bodies with friction. After the description of the problem we present its discretization based on linear or bilinear finite elements. The semi--smooth Newton method is used to find the solution, from which we derive active sets algorithms. Finally, we arrive at the globally convergent dual implementation of the algorithms in terms of the Langrange multipliers for the Tresca problem. Numerical experiments conclude the paper.

Solution of problems with material nonlinearities with a coupled finite element/boundary element scheme using an iterative solver. Yucca Mountain Site Characterization Project

International Nuclear Information System (INIS)

Koteras, J.R.

1996-01-01

The prediction of stresses and displacements around tunnels buried deep within the earth is an important class of geomechanics problems. The material behavior immediately surrounding the tunnel is typically nonlinear. The surrounding mass, even if it is nonlinear, can usually be characterized by a simple linear elastic model. The finite element method is best suited for modeling nonlinear materials of limited volume, while the boundary element method is well suited for modeling large volumes of linear elastic material. A computational scheme that couples the finite element and boundary element methods would seem particularly useful for geomechanics problems. A variety of coupling schemes have been proposed, but they rely on direct solution methods. Direct solution techniques have large storage requirements that become cumbersome for large-scale three-dimensional problems. An alternative to direct solution methods is iterative solution techniques. A scheme has been developed for coupling the finite element and boundary element methods that uses an iterative solution method. This report shows that this coupling scheme is valid for problems where nonlinear material behavior occurs in the finite element region

Application of power series to the solution of the boundary value problem for a second order nonlinear differential equation

International Nuclear Information System (INIS)

Semenova, V.N.

2016-01-01

A boundary value problem for a nonlinear second order differential equation has been considered. A numerical method has been proposed to solve this problem using power series. Results of numerical experiments have been presented in the paper [ru

Optimization of lift gas allocation in a gas lifted oil field as non-linear optimization problem

Directory of Open Access Journals (Sweden)

Roshan Sharma

2012-01-01

Full Text Available Proper allocation and distribution of lift gas is necessary for maximizing total oil production from a field with gas lifted oil wells. When the supply of the lift gas is limited, the total available gas should be optimally distributed among the oil wells of the field such that the total production of oil from the field is maximized. This paper describes a non-linear optimization problem with constraints associated with the optimal distribution of the lift gas. A non-linear objective function is developed using a simple dynamic model of the oil field where the decision variables represent the lift gas flow rate set points of each oil well of the field. The lift gas optimization problem is solved using the emph'fmincon' solver found in MATLAB. As an alternative and for verification, hill climbing method is utilized for solving the optimization problem. Using both of these methods, it has been shown that after optimization, the total oil production is increased by about 4. For multiple oil wells sharing lift gas from a common source, a cascade control strategy along with a nonlinear steady state optimizer behaves as a self-optimizing control structure when the total supply of lift gas is assumed to be the only input disturbance present in the process. Simulation results show that repeated optimization performed after the first time optimization under the presence of the input disturbance has no effect in the total oil production.

Uniqueness of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations

International Nuclear Information System (INIS)

Tran Duc Van

1994-01-01

The notion of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations is presented, some uniqueness theorems and a stability result are established by the method based on the theory of differential inclusions. In particular, the answer to an open problem of S.N. Kruzhkov is given. (author). 10 refs, 1 fig

Analysis of a Dynamic Viscoelastic Contact Problem with Normal Compliance, Normal Damped Response, and Nonmonotone Slip Rate Dependent Friction

Directory of Open Access Journals (Sweden)

Mikaël Barboteu

2016-01-01

Full Text Available We consider a mathematical model which describes the dynamic evolution of a viscoelastic body in frictional contact with an obstacle. The contact is modelled with a combination of a normal compliance and a normal damped response law associated with a slip rate-dependent version of Coulomb’s law of dry friction. We derive a variational formulation and an existence and uniqueness result of the weak solution of the problem is presented. Next, we introduce a fully discrete approximation of the variational problem based on a finite element method and on an implicit time integration scheme. We study this fully discrete approximation schemes and bound the errors of the approximate solutions. Under regularity assumptions imposed on the exact solution, optimal order error estimates are derived for the fully discrete solution. Finally, after recalling the solution of the frictional contact problem, some numerical simulations are provided in order to illustrate both the behavior of the solution related to the frictional contact conditions and the theoretical error estimate result.

Tyre-road contact using a particle-envelope surface model

Science.gov (United States)

Pinnington, Roger J.

2013-12-01

Determination of the contact forces is the central problem in all aspects of road-tyre interaction: i.e. noise, energy loss and friction. A procedure to find the contact forces under a rolling tyre is presented in four stages. First, the contact stiffness of a uniform peak array from indentations in the rubber tread, and also tyre carcass deflection, is described by some new simplified expressions. Second, a routine divides a single surface profile into equal search intervals, in which the highest peaks are identified. These are used to obtain the parameters for the interval, i.e. the mean envelope and the mean interval. The process is repeated at geometrically decreasing search intervals until the level of the data resolution, thereby describing the profile by a set of envelopes. The ‘strip profile’ ultimately used to describe the surface, is obtained by selecting the highest points across the profiles of one stone's width. The third stage is to combine the strip profile envelopes with the contact stiffness expressions, yielding the nonlinear stiffness-displacement, and force-displacement relationships for the chosen road-tyre combination. Finally the contact pressure distribution from a steady-state rolling tyre model is applied to the strip profile, via the force-displacement relationship, giving the local tyre displacements on the road texture. This displacement pattern is shown to be proportional to the time and space varying contact pressure, which then is incorporated into a wave equation for rolling contact.

Spike-layer solutions to nonlinear fractional Schrodinger equations with almost optimal nonlinearities

Directory of Open Access Journals (Sweden)

Jinmyoung Seok

2015-07-01

Full Text Available In this article, we are interested in singularly perturbed nonlinear elliptic problems involving a fractional Laplacian. Under a class of nonlinearity which is believed to be almost optimal, we construct a positive solution which exhibits multiple spikes near any given local minimum components of an exterior potential of the problem.

Existence of weak solutions to a nonlinear reaction-diffusion system with singular sources

Directory of Open Access Journals (Sweden)

Ida de Bonis

2017-09-01

Full Text Available We discuss the existence of a class of weak solutions to a nonlinear parabolic system of reaction-diffusion type endowed with singular production terms by reaction. The singularity is due to a potential occurrence of quenching localized to the domain boundary. The kind of quenching we have in mind is due to a twofold contribution: (i the choice of boundary conditions, modeling in our case the contact with an infinite reservoir filled with ready-to-react chemicals and (ii the use of a particular nonlinear, non-Lipschitz structure of the reaction kinetics. Our working techniques use fine energy estimates for approximating non-singular problems and uniform control on the set where singularities are localizing.

Rigid Bodies in Contact

DEFF Research Database (Denmark)

Niebe, Sarah Maria

. A contact point determination method, based on boolean surface maps, is developed to handle collisions between tetrahedral meshes. The novel nonsmooth nonlinear conjugate gradient (NNCG) method is presented. The NNCG method is comparable in terms of accuracy to the state-of-the-art method, projected Gauss...

Coupled nonlinear oscillators

Energy Technology Data Exchange (ETDEWEB)

Chandra, J; Scott, A C

1983-01-01

Topics discussed include transitions in weakly coupled nonlinear oscillators, singularly perturbed delay-differential equations, and chaos in simple laser systems. Papers are presented on truncated Navier-Stokes equations in a two-dimensional torus, on frequency locking in Josephson point contacts, and on soliton excitations in Josephson tunnel junctions. Attention is also given to the nonlinear coupling of radiation pulses to absorbing anharmonic molecular media, to aspects of interrupted coarse-graining in stimulated excitation, and to a statistical analysis of long-term dynamic irregularity in an exactly soluble quantum mechanical model.

Analytical and Computational Modeling of Mechanical Waves in Microscale Granular Crystals: Nonlinearity and Rotational Dynamics

Science.gov (United States)

Wallen, Samuel P.

Granular media are one of the most common, yet least understood forms of matter on earth. The difficulties in understanding the physics of granular media stem from the fact that they are typically heterogeneous and highly disordered, and the grains interact via nonlinear contact forces. Historically, one approach to reducing these complexities and gaining new insight has been the study of granular crystals, which are ordered arrays of similarly-shaped particles (typically spheres) in Hertzian contact. Using this setting, past works explored the rich nonlinear dynamics stemming from contact forces, and proposed avenues where such granular crystals could form designer, dynamically responsive materials, which yield beneficial functionality in dynamic regimes. In recent years, the combination of self-assembly fabrication methods and laser ultrasonic experimental characterization have enabled the study of granular crystals at microscale. While our intuition may suggest that these microscale granular crystals are simply scaled-down versions of their macroscale counterparts, in fact, the relevant physics change drastically; for example, short-range adhesive forces between particles, which are negligible at macroscale, are several orders of magnitude stronger than gravity at microscale. In this thesis, we present recent advances in analytical and computational modeling of microscale granular crystals, in particular concerning the interplay of nonlinearity, shear interactions, and particle rotations, which have previously been either absent, or included separately at macroscale. Drawing inspiration from past works on phononic crystals and nonlinear lattices, we explore problems involving locally-resonant metamaterials, nonlinear localized modes, amplitude-dependent energy partition, and other rich dynamical phenomena. This work enhances our understanding of microscale granular media, which may find applicability in fields such as ultrasonic wave tailoring, signal processing

A Unique Method to Describe the Bonding Strength in a Bonded Solid–Solid Interface by Contact Acoustic Nonlinearity

International Nuclear Information System (INIS)

Jian-Jun, Chen; De, Zhang; Yi-Wei, Mao; Jian-Chun, Cheng

2009-01-01

We present a unique method to describe the bonding strength at a bonded solid–solid interface in a multilayered composite material by contact acoustic nonlinearity (CAN) parameter. A CAN model on the bonded solid–solid interface is depicted. It can be seen from the model that CAN parameter is very sensitive to the bonding strength at the interface. When an incident focusing acoustic longitudinal wave scans the interface in two dimensions, the transmitted wave can be used to extract CAN parameter. The contour of the bonding strength for a sample is obtained by CAN parameter. The results show that the region with weak bonding strength can be easily distinguished from the contour

Existence of solutions for the dynamic frictional contact problem of isotropic viscoelastic bodies

Czech Academy of Sciences Publication Activity Database

Eck, C.; Jarušek, Jiří

2003-01-01

Roč. 53, č. 2 (2003), s. 157-181 ISSN 0362-546X R&D Projects: GA AV ČR IAA1075005; GA AV ČR IAA1075707 Institutional research plan: CEZ:AV0Z1019905; CEZ:AV0Z1019905 Keywords : dynamic contact problem * parabolic equation * Coulomb law of friction Subject RIV: BA - General Mathematics Impact factor: 0.354, year: 2003

Nonlinear Damping Identification in Nonlinear Dynamic System Based on Stochastic Inverse Approach

Directory of Open Access Journals (Sweden)

S. L. Han

2012-01-01

Full Text Available The nonlinear model is crucial to prepare, supervise, and analyze mechanical system. In this paper, a new nonparametric and output-only identification procedure for nonlinear damping is studied. By introducing the concept of the stochastic state space, we formulate a stochastic inverse problem for a nonlinear damping. The solution of the stochastic inverse problem is designed as probabilistic expression via the hierarchical Bayesian formulation by considering various uncertainties such as the information insufficiency in parameter of interests or errors in measurement. The probability space is estimated using Markov chain Monte Carlo (MCMC. The applicability of the proposed method is demonstrated through numerical experiment and particular application to a realistic problem related to ship roll motion.

Non-intrusive reduced order modeling of nonlinear problems using neural networks

Science.gov (United States)

Hesthaven, J. S.; Ubbiali, S.

2018-06-01

We develop a non-intrusive reduced basis (RB) method for parametrized steady-state partial differential equations (PDEs). The method extracts a reduced basis from a collection of high-fidelity solutions via a proper orthogonal decomposition (POD) and employs artificial neural networks (ANNs), particularly multi-layer perceptrons (MLPs), to accurately approximate the coefficients of the reduced model. The search for the optimal number of neurons and the minimum amount of training samples to avoid overfitting is carried out in the offline phase through an automatic routine, relying upon a joint use of the Latin hypercube sampling (LHS) and the Levenberg-Marquardt (LM) training algorithm. This guarantees a complete offline-online decoupling, leading to an efficient RB method - referred to as POD-NN - suitable also for general nonlinear problems with a non-affine parametric dependence. Numerical studies are presented for the nonlinear Poisson equation and for driven cavity viscous flows, modeled through the steady incompressible Navier-Stokes equations. Both physical and geometrical parametrizations are considered. Several results confirm the accuracy of the POD-NN method and show the substantial speed-up enabled at the online stage as compared to a traditional RB strategy.

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods. Ph.D. Thesis - Paris Univ.; [mathematical models of transonic flow using nonlinear equations

Science.gov (United States)

Periaux, J.

1979-01-01

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Mathematical and numerical study of nonlinear boundary problems related to plasma physics

International Nuclear Information System (INIS)

Sermange, M.

1982-06-01

After the study of some equations based on the Hodgkin-Huxley model, the work presented here is concerned with nonlinear boundary problems in MHD. They are gathered in two subjects: equilibrium equations and stability equations. The axisymmetric MHD equilibrium equations with free boundary have been studied by different authors, particularly the existence, regularity, unicity and non-unicity. Here, bifurcation, convergence of calculation methods existence of solutions in a discontinuous frame are studied. MHD stability can be determined by the principle of Bernstein et al; the mathematical work concerned here bears on the equivalence, in the case of two-dimensional or axisymmetric stability, between this model and a scalar eigenvalue problem which is introduced. At last, modules for computing MHD equilibrium for the simulation of plasma confinement in a tokamak are described [fr

A limited memory BFGS method for a nonlinear inverse problem in digital breast tomosynthesis

Science.gov (United States)

Landi, G.; Loli Piccolomini, E.; Nagy, J. G.

2017-09-01

Digital breast tomosynthesis (DBT) is an imaging technique that allows the reconstruction of a pseudo three-dimensional image of the breast from a finite number of low-dose two-dimensional projections obtained by different x-ray tube angles. An issue that is often ignored in DBT is the fact that an x-ray beam is polyenergetic, i.e. it is composed of photons with different levels of energy. The polyenergetic model requires solving a large-scale, nonlinear inverse problem, which is more expensive than the typically used simplified, linear monoenergetic model. However, the polyenergetic model is much less susceptible to beam hardening artifacts, which show up as dark streaks and cupping (i.e. background nonuniformities) in the reconstructed image. In addition, it has been shown that the polyenergetic model can be exploited to obtain additional quantitative information about the material of the object being imaged. In this paper we consider the multimaterial polyenergetic DBT model, and solve the nonlinear inverse problem with a limited memory BFGS quasi-Newton method. Regularization is enforced at each iteration using a diagonally modified approximation of the Hessian matrix, and by truncating the iterations.

An analytic solution to the homogeneous EIT problem on the 2D disk and its application to estimation of electrode contact impedances

International Nuclear Information System (INIS)

Demidenko, Eugene

2011-01-01

An analytic solution of the potential distribution on a 2D homogeneous disk for electrical impedance tomography under the complete electrode model is expressed via an infinite system of linear equations. For the shunt electrode model with two electrodes, our solution coincides with the previously derived solution expressed via elliptic integral (Pidcock et al 1995 Physiol. Meas. 16 77–90). The Dirichlet-to-Neumann map is derived for statistical estimation via nonlinear least squares. The solution is validated in phantom experiments and applied for breast contact impedance estimation in vivo. Statistical hypothesis testing is used to test whether the contact impedances are the same across electrodes or all equal zero. Our solution can be especially useful for a rapid real-time test for bad surface contact in clinical setting

On Ni/Au Alloyed Contacts to Mg-Doped GaN

Science.gov (United States)

Sarkar, Biplab; Reddy, Pramod; Klump, Andrew; Kaess, Felix; Rounds, Robert; Kirste, Ronny; Mita, Seiji; Kohn, Erhard; Collazo, Ramon; Sitar, Zlatko

2018-01-01

Ni/Au contacts to p-GaN were studied as a function of free hole concentration in GaN using planar transmission line measurement structures. All contacts showed a nonlinear behavior, which became stronger for lower doping concentrations. Electrical and structural analysis indicated that the current conduction between the contact and the p-GaN was through localized nano-sized clusters. Thus, the non-linear contact behavior can be well explained using the alloyed contact model. Two contributions to the contact resistance were identified: the spreading resistance in the semiconductor developed by the current crowding around the electrically active clusters, and diode-type behavior at the interface of the electrically active clusters with the semiconductor. Hence, the equivalent Ni/Au contact model consists of a diode and a resistor in series for each active cluster. The reduced barrier height observed in the measurements is thought to be generated by the extraction of Ga from the crystalline surface and localized formation of the Au:Ga phase. The alloyed contact analyses presented in this work are in good agreement with some of the commonly observed behavior of similar contacts described in the literature.

Nonlinear approaches in engineering applications 2

CERN Document Server

Jazar, Reza N

2013-01-01

Provides updated principles and applications of the nonlinear approaches in solving engineering and physics problems Demonstrates how nonlinear approaches may open avenues to better, safer, cheaper systems with less energy consumption Has a strong emphasis on the application, physical meaning, and methodologies of nonlinear approaches in different engineering and science problems

Nonlinear dynamics of quadratically cubic systems

International Nuclear Information System (INIS)

Rudenko, O V

2013-01-01

We propose a modified form of the well-known nonlinear dynamic equations with quadratic relations used to model a cubic nonlinearity. We show that such quadratically cubic equations sometimes allow exact solutions and sometimes make the original problem easier to analyze qualitatively. Occasionally, exact solutions provide a useful tool for studying new phenomena. Examples considered include nonlinear ordinary differential equations and Hopf, Burgers, Korteweg–de Vries, and nonlinear Schrödinger partial differential equations. Some problems are solved exactly in the space–time and spectral representations. Unsolved problems potentially solvable by the proposed approach are listed. (methodological notes)

Non-linear time series extreme events and integer value problems

CERN Document Server

Turkman, Kamil Feridun; Zea Bermudez, Patrícia

2014-01-01

This book offers a useful combination of probabilistic and statistical tools for analyzing nonlinear time series. Key features of the book include a study of the extremal behavior of nonlinear time series and a comprehensive list of nonlinear models that address different aspects of nonlinearity. Several inferential methods, including quasi likelihood methods, sequential Markov Chain Monte Carlo Methods and particle filters, are also included so as to provide an overall view of the available tools for parameter estimation for nonlinear models. A chapter on integer time series models based on several thinning operations, which brings together all recent advances made in this area, is also included. Readers should have attended a prior course on linear time series, and a good grasp of simulation-based inferential methods is recommended. This book offers a valuable resource for second-year graduate students and researchers in statistics and other scientific areas who need a basic understanding of nonlinear time ...

Nonlinear robust hierarchical control for nonlinear uncertain systems

Directory of Open Access Journals (Sweden)

Leonessa Alexander

1999-01-01

Full Text Available A nonlinear robust control-system design framework predicated on a hierarchical switching controller architecture parameterized over a set of moving nominal system equilibria is developed. Specifically, using equilibria-dependent Lyapunov functions, a hierarchical nonlinear robust control strategy is developed that robustly stabilizes a given nonlinear system over a prescribed range of system uncertainty by robustly stabilizing a collection of nonlinear controlled uncertain subsystems. The robust switching nonlinear controller architecture is designed based on a generalized (lower semicontinuous Lyapunov function obtained by minimizing a potential function over a given switching set induced by the parameterized nominal system equilibria. The proposed framework robustly stabilizes a compact positively invariant set of a given nonlinear uncertain dynamical system with structured parametric uncertainty. Finally, the efficacy of the proposed approach is demonstrated on a jet engine propulsion control problem with uncertain pressure-flow map data.

Solving Boundary Value Problem for a Nonlinear Stationary Controllable System with Synthesizing Control

Directory of Open Access Journals (Sweden)

Alexander N. Kvitko

2017-01-01

Full Text Available An algorithm for constructing a control function that transfers a wide class of stationary nonlinear systems of ordinary differential equations from an initial state to a final state under certain control restrictions is proposed. The algorithm is designed to be convenient for numerical implementation. A constructive criterion of the desired transfer possibility is presented. The problem of an interorbital flight is considered as a test example and it is simulated numerically with the presented method.

Measurement and fitting techniques for the assessment of material nonlinearity using nonlinear Rayleigh waves

Energy Technology Data Exchange (ETDEWEB)

Torello, David [GW Woodruff School of Mechanical Engineering, Georgia Tech (United States); Kim, Jin-Yeon [School of Civil and Environmental Engineering, Georgia Tech (United States); Qu, Jianmin [Department of Civil and Environmental Engineering, Northwestern University (United States); Jacobs, Laurence J. [School of Civil and Environmental Engineering, Georgia Tech and GW Woodruff School of Mechanical Engineering, Georgia Tech (United States)

2015-03-31

This research considers the effects of diffraction, attenuation, and the nonlinearity of generating sources on measurements of nonlinear ultrasonic Rayleigh wave propagation. A new theoretical framework for correcting measurements made with air-coupled and contact piezoelectric receivers for the aforementioned effects is provided based on analytical models and experimental considerations. A method for extracting the nonlinearity parameter β{sub 11} is proposed based on a nonlinear least squares curve-fitting algorithm that is tailored for Rayleigh wave measurements. Quantitative experiments are conducted to confirm the predictions for the nonlinearity of the piezoelectric source and to demonstrate the effectiveness of the curve-fitting procedure. These experiments are conducted on aluminum 2024 and 7075 specimens and a β{sub 11}{sup 7075}/β{sub 11}{sup 2024} measure of 1.363 agrees well with previous literature and earlier work.

Contact kinematics of biomimetic scales

Energy Technology Data Exchange (ETDEWEB)

Ghosh, Ranajay; Ebrahimi, Hamid; Vaziri, Ashkan, E-mail: vaziri@coe.neu.edu [Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115 (United States)

2014-12-08

Dermal scales, prevalent across biological groups, considerably boost survival by providing multifunctional advantages. Here, we investigate the nonlinear mechanical effects of biomimetic scale like attachments on the behavior of an elastic substrate brought about by the contact interaction of scales in pure bending using qualitative experiments, analytical models, and detailed finite element (FE) analysis. Our results reveal the existence of three distinct kinematic phases of operation spanning linear, nonlinear, and rigid behavior driven by kinematic interactions of scales. The response of the modified elastic beam strongly depends on the size and spatial overlap of rigid scales. The nonlinearity is perceptible even in relatively small strain regime and without invoking material level complexities of either the scales or the substrate.

Contact stiffness considerations when simulating tyre/road noise

Science.gov (United States)

Winroth, Julia; Kropp, Wolfgang; Hoever, Carsten; Höstmad, Patrik

2017-11-01

Tyre/road simulation tools that can capture tyre vibrations, rolling resistance and noise generation are useful for understanding the complex processes that are involved and thereby promoting further development and optimisation. The most detailed tyre/road contact models use a spatial discretisation of the contact and assume an interfacial stiffness to account for the small-scale roughness within the elements. This interfacial stiffness has been found to have a significant impact on the simulated noise emissions but no thorough investigations of this sensitivity have been conducted. Three mechanisms are thought to be involved: The horn effect, the modal composition of the vibrational field of the tyre and the contact forces exciting the tyre vibrations. This study used a numerical tyre/road noise simulation tool based on physical relations to investigate these aspects. The model includes a detailed time-domain contact model with linear or non-linear contact springs that accounts for the effect of local tread deformation on smaller length scales. Results confirm that an increase in contact spring stiffness causes a significant increase of the simulated tyre/road noise. This is primarily caused by a corresponding increase in the contact forces, resulting in larger vibrational amplitudes. The horn effect and the modal composition are relatively unaffected and have minor effects on the radiated noise. A more detailed non-linear contact spring formulation with lower stiffness at small indentations results in a reduced high-frequency content in the contact forces and the simulated noise.

Frictional Compliant Haptic Contact and Deformation of Soft Objects

Directory of Open Access Journals (Sweden)

Naci Zafer

2016-05-01

Full Text Available This paper is concerned with compliant haptic contact and deformation of soft objects. A human soft fingertip model is considered to act as the haptic interface and is brought into contact with and deforms a discrete surface. A nonlinear constitutive law is developed in predicting normal forces and, for the haptic display of surface texture, motions along the surface are also resisted at various rates by accounting for dynamic Lund-Grenoble (LuGre frictional forces. For the soft fingertip to apply forces over an area larger than a point, normal and frictional forces are distributed around the soft fingertip contact location on the deforming surface. The distribution is realized based on a kernel smoothing function and by a nonlinear spring-damper net around the contact point. Experiments conducted demonstrate the accuracy and effectiveness of our approach in real-time haptic rendering of a kidney surface. The resistive (interaction forces are applied at the user fingertip bone edge. A 3-DoF parallel robotic manipulator equipped with a constraint based controller is used for the implementation. By rendering forces both in lateral and normal directions, the designed haptic interface system allows the user to realistically feel both the geometrical and mechanical (nonlinear properties of the deforming kidney.

Behavior of the maximal solution of the Cauchy problem for some nonlinear pseudoparabolic equation as $|x|oinfty$

Directory of Open Access Journals (Sweden)

Tatiana Kavitova

2012-08-01

Full Text Available We prove a comparison principle for solutions of the Cauchy problem of the nonlinear pseudoparabolic equation $u_t=Delta u_t+ Deltavarphi(u +h(t,u$ with nonnegative bounded initial data. We show stabilization of a maximal solution to a maximal solution of the Cauchy problem for the corresponding ordinary differential equation $vartheta'(t=h(t,vartheta$ as $|x|oinfty$ under certain conditions on an initial datum.

Nonlinear crack mechanics

International Nuclear Information System (INIS)

Khoroshun, L.P.

1995-01-01

The characteristic features of the deformation and failure of actual materials in the vicinity of a crack tip are due to their physical nonlinearity in the stress-concentration zone, which is a result of plasticity, microfailure, or a nonlinear dependence of the interatomic forces on the distance. Therefore, adequate models of the failure mechanics must be nonlinear, in principle, although linear failure mechanics is applicable if the zone of nonlinear deformation is small in comparison with the crack length. Models of crack mechanics are based on analytical solutions of the problem of the stress-strain state in the vicinity of the crack. On account of the complexity of the problem, nonlinear models are bason on approximate schematic solutions. In the Leonov-Panasyuk-Dugdale nonlinear model, one of the best known, the actual two-dimensional plastic zone (the nonlinearity zone) is replaced by a narrow one-dimensional zone, which is then modeled by extending the crack with a specified normal load equal to the yield point. The condition of finite stress is applied here, and hence the length of the plastic zone is determined. As a result of this approximation, the displacement in the plastic zone at the abscissa is nonzero

Existence of solutions to nonlinear parabolic unilateral problems with an obstacle depending on time

Directory of Open Access Journals (Sweden)

Nabila Bellal

2014-10-01

Full Text Available Using the penalty method, we prove the existence of solutions to nonlinear parabolic unilateral problems with an obstacle depending on time. To find a solution, the original inequality is transformed into an equality by adding a positive function on the right-hand side and a complementary condition. This result can be seen as a generalization of the results by Mokrane in [11] where the obstacle is zero.

Intuitionistic Fuzzy Goal Programming Technique for Solving Non-Linear Multi-objective Structural Problem

Directory of Open Access Journals (Sweden)

Samir Dey

2015-07-01

Full Text Available This paper proposes a new multi-objective intuitionistic fuzzy goal programming approach to solve a multi-objective nonlinear programming problem in context of a structural design. Here we describe some basic properties of intuitionistic fuzzy optimization. We have considered a multi-objective structural optimization problem with several mutually conflicting objectives. The design objective is to minimize weight of the structure and minimize the vertical deflection at loading point of a statistically loaded three-bar planar truss subjected to stress constraints on each of the truss members. This approach is used to solve the above structural optimization model based on arithmetic mean and compare with the solution by intuitionistic fuzzy goal programming approach. A numerical solution is given to illustrate our approach.

Discrete and continuum links to a nonlinear coupled transport problem of interacting populations

Science.gov (United States)

Duong, M. H.; Muntean, A.; Richardson, O. M.

2017-07-01

We are interested in exploring interacting particle systems that can be seen as microscopic models for a particular structure of coupled transport flux arising when different populations are jointly evolving. The scenarios we have in mind are inspired by the dynamics of pedestrian flows in open spaces and are intimately connected to cross-diffusion and thermo-diffusion problems holding a variational structure. The tools we use include a suitable structure of the relative entropy controlling TV-norms, the construction of Lyapunov functionals and particular closed-form solutions to nonlinear transport equations, a hydrodynamics limiting procedure due to Philipowski, as well as the construction of numerical approximates to both the continuum limit problem in 2D and to the original interacting particle systems.

Solving nonlinear nonstationary problem of heat-conductivity by finite element method

Directory of Open Access Journals (Sweden)

Антон Янович Карвацький

2016-11-01

Full Text Available Methodology and effective solving algorithm of non-linear dynamic problems of thermal and electric conductivity with significant temperature dependence of thermal and physical properties are given on the basis of finite element method (FEM and Newton linearization method. Discrete equations system FEM was obtained with the use of Galerkin method, where the main function is the finite element form function. The methodology based on successive solving problems of thermal and electrical conductivity has been examined in the work in order to minimize the requirements for calculating resources (RAM. in particular. Having used Mathcad software original programming code was developed to solve the given problem. After investigation of the received results, comparative analyses of accurate solution data and results of numerical solutions, obtained with the use of Matlab programming products, was held. The geometry of one fourth part of the finite sized cylinder was used to test the given numerical model. The discretization of the calculation part was fulfilled using the open programming software for automated Gmsh nets with tetrahedral units, while ParaView, which is an open programming code as well, was used to visualize the calculation results. It was found out that the maximum value violation of potential and temperature determination doesn`t exceed 0,2-0,83% in the given work according to the problem conditions

Parallel Solution of Robust Nonlinear Model Predictive Control Problems in Batch Crystallization

Directory of Open Access Journals (Sweden)

Yankai Cao

2016-06-01

Full Text Available Representing the uncertainties with a set of scenarios, the optimization problem resulting from a robust nonlinear model predictive control (NMPC strategy at each sampling instance can be viewed as a large-scale stochastic program. This paper solves these optimization problems using the parallel Schur complement method developed to solve stochastic programs on distributed and shared memory machines. The control strategy is illustrated with a case study of a multidimensional unseeded batch crystallization process. For this application, a robust NMPC based on min–max optimization guarantees satisfaction of all state and input constraints for a set of uncertainty realizations, and also provides better robust performance compared with open-loop optimal control, nominal NMPC, and robust NMPC minimizing the expected performance at each sampling instance. The performance of robust NMPC can be improved by generating optimization scenarios using Bayesian inference. With the efficient parallel solver, the solution time of one optimization problem is reduced from 6.7 min to 0.5 min, allowing for real-time application.

Multiple Solutions of Nonlinear Boundary Value Problems of Fractional Order: A New Analytic Iterative Technique

Directory of Open Access Journals (Sweden)

Omar Abu Arqub

2014-01-01

Full Text Available The purpose of this paper is to present a new kind of analytical method, the so-called residual power series, to predict and represent the multiplicity of solutions to nonlinear boundary value problems of fractional order. The present method is capable of calculating all branches of solutions simultaneously, even if these multiple solutions are very close and thus rather difficult to distinguish even by numerical techniques. To verify the computational efficiency of the designed proposed technique, two nonlinear models are performed, one of them arises in mixed convection flows and the other one arises in heat transfer, which both admit multiple solutions. The results reveal that the method is very effective, straightforward, and powerful for formulating these multiple solutions.

Introduction to nonlinear dispersive equations

CERN Document Server

Linares, Felipe

2015-01-01

This textbook introduces the well-posedness theory for initial-value problems of nonlinear, dispersive partial differential equations, with special focus on two key models, the Korteweg–de Vries equation and the nonlinear Schrödinger equation. A concise and self-contained treatment of background material (the Fourier transform, interpolation theory, Sobolev spaces, and the linear Schrödinger equation) prepares the reader to understand the main topics covered: the initial-value problem for the nonlinear Schrödinger equation and the generalized Korteweg–de Vries equation, properties of their solutions, and a survey of general classes of nonlinear dispersive equations of physical and mathematical significance. Each chapter ends with an expert account of recent developments and open problems, as well as exercises. The final chapter gives a detailed exposition of local well-posedness for the nonlinear Schrödinger equation, taking the reader to the forefront of recent research. The second edition of Introdu...

Analysis of contact stiffness in ultrasound atomic force microscopy: three-dimensional time-dependent ultrasound modeling

International Nuclear Information System (INIS)

Piras, Daniele; Sadeghian, Hamed

2017-01-01

Ultrasound atomic force microscopy (US-AFM) has been used for subsurface imaging of nanostructures. The contact stiffness variations have been suggested as the origin of the image contrast. Therefore, to analyze the image contrast, the local changes in the contact stiffness due to the presence of subsurface features should be calculated. So far, only static simulations have been conducted to analyze the local changes in the contact stiffness and, consequently, the contrast in US-AFM. Such a static approach does not fully represent the real US-AFM experiment, where an ultrasound wave is launched either into the sample or at the tip, which modulates the contact stiffness. This is a time-dependent nonlinear dynamic problem rather than a static and stationary one. This paper presents dynamic 3D ultrasound analysis of contact stiffness in US-AFM (in contrast to static analysis) to realistically predict the changes in contact stiffness and thus the changes in the subsurface image contrast. The modulation frequency also influences the contact stiffness variations and, thus, the image contrast. The three-dimensional time-dependent ultrasound analysis will greatly aid in the contrast optimization of subsurface nano imaging with US-AFM. (paper)

Nonlinear evolution equations

CERN Document Server

Uraltseva, N N

1995-01-01

This collection focuses on nonlinear problems in partial differential equations. Most of the papers are based on lectures presented at the seminar on partial differential equations and mathematical physics at St. Petersburg University. Among the topics explored are the existence and properties of solutions of various classes of nonlinear evolution equations, nonlinear imbedding theorems, bifurcations of solutions, and equations of mathematical physics (Navier-Stokes type equations and the nonlinear Schrödinger equation). The book will be useful to researchers and graduate students working in p

Mathematical modeling and applications in nonlinear dynamics

CERN Document Server

Merdan, Hüseyin

2016-01-01

The book covers nonlinear physical problems and mathematical modeling, including molecular biology, genetics, neurosciences, artificial intelligence with classical problems in mechanics and astronomy and physics. The chapters present nonlinear mathematical modeling in life science and physics through nonlinear differential equations, nonlinear discrete equations and hybrid equations. Such modeling can be effectively applied to the wide spectrum of nonlinear physical problems, including the KAM (Kolmogorov-Arnold-Moser (KAM)) theory, singular differential equations, impulsive dichotomous linear systems, analytical bifurcation trees of periodic motions, and almost or pseudo- almost periodic solutions in nonlinear dynamical systems. Provides methods for mathematical models with switching, thresholds, and impulses, each of particular importance for discontinuous processes Includes qualitative analysis of behaviors on Tumor-Immune Systems and methods of analysis for DNA, neural networks and epidemiology Introduces...

A Procedure for 3-D Contact Stress Analysis of Spiral Bevel Gears

Science.gov (United States)

Kumar, A.; Bibel, G.

1994-01-01

Contact stress distribution of spiral bevel gears using nonlinear finite element static analysis is presented. Procedures have been developed to solve the nonlinear equations that identify the gear and pinion surface coordinates based on the kinematics of the cutting process and orientate the pinion and the gear in space to mesh with each other. Contact is simulated by connecting GAP elements along the intersection of a line from each pinion point (parallel to the normal at the contact point) with the gear surface. A three dimensional model with four gear teeth and three pinion teeth is used to determine the contact stresses at two different contact positions in a spiral bevel gearset. A summary of the elliptical contact stress distribution is given. This information will be helpful to helicopter and aircraft transmission designers who need to minimize weight of the transmission and maximize reliability.

A uniformly valid approximation algorithm for nonlinear ordinary singular perturbation problems with boundary layer solutions.

Science.gov (United States)

Cengizci, Süleyman; Atay, Mehmet Tarık; Eryılmaz, Aytekin

2016-01-01

This paper is concerned with two-point boundary value problems for singularly perturbed nonlinear ordinary differential equations. The case when the solution only has one boundary layer is examined. An efficient method so called Successive Complementary Expansion Method (SCEM) is used to obtain uniformly valid approximations to this kind of solutions. Four test problems are considered to check the efficiency and accuracy of the proposed method. The numerical results are found in good agreement with exact and existing solutions in literature. The results confirm that SCEM has a superiority over other existing methods in terms of easy-applicability and effectiveness.

Feedback control for distributed heat transfer mechanisms in direct-contact membrane distillation system

KAUST Repository

Eleiwi, Fadi; N'Doye, Ibrahima; Laleg-Kirati, Taous-Meriem

2015-01-01

In this paper, the problem of stabilization and production rate reference tracking for a Direct-Contact Membrane Distillation (DCMD) system is addressed. Sufficient conditions for the asymptotic and exponential stabilization for DCMD system are presented using the Gronwall-Bellman lemma and Linear Matrix Inequalities (LMIs) approaches, respectively. A nonlinear observer is then proposed to estimate the temperature distribution among the DCMD domain. This contributes to propose a reference production rate control design for the DCMD process via observer-based output control approach. Finally, numerical simulations are given to show the effectiveness of the proposed methods.

Feedback control for distributed heat transfer mechanisms in direct-contact membrane distillation system

KAUST Repository

Eleiwi, Fadi

2015-09-21

In this paper, the problem of stabilization and production rate reference tracking for a Direct-Contact Membrane Distillation (DCMD) system is addressed. Sufficient conditions for the asymptotic and exponential stabilization for DCMD system are presented using the Gronwall-Bellman lemma and Linear Matrix Inequalities (LMIs) approaches, respectively. A nonlinear observer is then proposed to estimate the temperature distribution among the DCMD domain. This contributes to propose a reference production rate control design for the DCMD process via observer-based output control approach. Finally, numerical simulations are given to show the effectiveness of the proposed methods.

One-dimensional singular problems involving the p-Laplacian and nonlinearities indefinite in sign

OpenAIRE

Kaufmann, Uriel; Medri, Iván

2015-01-01

Let $\\Omega$ be a bounded open interval, let $p>1$ and $\\gamma>0$, and let $m:\\Omega\\rightarrow\\mathbb{R}$ be a function that may change sign in $\\Omega $. In this article we study the existence and nonexistence of positive solutions for one-dimensional singular problems of the form $-(\\left\\vert u^{\\prime}\\right\\vert ^{p-2}u^{\\prime})^{\\prime}=m\\left( x\\right) u^{-\\gamma}$ in $\\Omega$, $u=0$ on $\\partial\\Omega$. As a consequence we also derive existence results for other related nonlinearities.

CASKETSS-HEAT: a finite difference computer program for nonlinear heat conduction problems

International Nuclear Information System (INIS)

Ikushima, Takeshi

1988-12-01

A heat conduction program CASKETSS-HEAT has been developed. CASKETSS-HEAT is a finite difference computer program used for the solution of multi-dimensional nonlinear heat conduction problems. Main features of CASKETSS-HEAT are as follows. (1) One, two and three-dimensional geometries for heat conduction calculation are available. (2) Convection and radiation heat transfer of boundry can be specified. (3) Phase change and chemical change can be treated. (4) Finned surface heat transfer can be treated easily. (5) Data memory allocation in the program is variable according to problem size. (6) The program is a compatible heat transfer analysis program to the stress analysis program SAP4 and SAP5. (7) Pre- and post-processing for input data generation and graphic representation of calculation results are available. In the paper, brief illustration of calculation method, input data and sample calculation are presented. (author)

Improved study of electric dipoles on the Si(100)-2 × 1 surface by non-contact scanning nonlinear dielectric microscopy

International Nuclear Information System (INIS)

Suzuki, Masataka; Yamasue, Kohei; Cho, Yasuo; Abe, Masayuki; Sugimoto, Yoshiaki

2014-01-01

We studied a Si(100)-2 × 1 surface by non-contact scanning nonlinear dielectric microscopy (NC-SNDM). Simultaneously taken images of the topography and electric dipole moment distribution show that negative electric dipole moments are locally formed on individual dimers on the surface. In addition, we obtained the dc bias voltage dependence of the ε local (3) signal on a specific dimer by using an atom-tracking technique with NC-SNDM. We observed that the electric dipole induced a surface potential of around −250 mV on the dimer.

Solutions of dissimilar material singularity and contact problems

International Nuclear Information System (INIS)

Yang, Y.

2003-09-01

Due to the mismatch of the material properties of joined components, after a homogeneous temperature change or under a mechanical loading, very high stresses occur near the intersection of the interface and the outer surface, or near the intersection of two interfaces. For most material combinations and joint geometries, there exists even a stress singularity. These high stresses may cause fracture of the joint. The investigation of the stress situation near the singular point, therefore, is of great interest. Especially, the relationship between the singular stress exponent, the material data and joint geometry is important for choosing a suitable material combination and joint geometry. In this work, the singular stress field is described analytically in case of the joint having a real and a complex eigenvalue. Solutions of different singularity problems are given, which are two dissimilar materials joint with free edges; dissimilar materials joint with edge tractions; joint with interface corner; joint with a given displacement at one edge; cracks in dissimilar materials joint; contact problem in dissimilar materials and logarithmic stress singularity. For an arbitrary joint geometry and material combination, the stress singular exponent, the angular function and the regular stress term can be calculated analytically. The stress intensity factors for a finite joint can be determined applying numerical methods, e.g. the finite element method (FEM). The method to determine more than one stress intensity factor is presented. The characteristics of the eigenvalues and the stress intensity factors are shown for different joint conditions. (orig.)

Initial-Boundary Value Problem Solution of the Nonlinear Shallow-water Wave Equations

Science.gov (United States)

Kanoglu, U.; Aydin, B.

2014-12-01

The hodograph transformation solutions of the one-dimensional nonlinear shallow-water wave (NSW) equations are usually obtained through integral transform techniques such as Fourier-Bessel transforms. However, the original formulation of Carrier and Greenspan (1958 J Fluid Mech) and its variant Carrier et al. (2003 J Fluid Mech) involve evaluation integrals. Since elliptic integrals are highly singular as discussed in Carrier et al. (2003), this solution methodology requires either approximation of the associated integrands by smooth functions or selection of regular initial/boundary data. It should be noted that Kanoglu (2004 J Fluid Mech) partly resolves this issue by simplifying the resulting integrals in closed form. Here, the hodograph transform approach is coupled with the classical eigenfunction expansion method rather than integral transform techniques and a new analytical model for nonlinear long wave propagation over a plane beach is derived. This approach is based on the solution methodology used in Aydın & Kanoglu (2007 CMES-Comp Model Eng) for wind set-down relaxation problem. In contrast to classical initial- or boundary-value problem solutions, here, the NSW equations are formulated to yield an initial-boundary value problem (IBVP) solution. In general, initial wave profile with nonzero initial velocity distribution is assumed and the flow variables are given in the form of Fourier-Bessel series. The results reveal that the developed method allows accurate estimation of the spatial and temporal variation of the flow quantities, i.e., free-surface height and depth-averaged velocity, with much less computational effort compared to the integral transform techniques such as Carrier et al. (2003), Kanoglu (2004), Tinti & Tonini (2005 J Fluid Mech), and Kanoglu & Synolakis (2006 Phys Rev Lett). Acknowledgments: This work is funded by project ASTARTE- Assessment, STrategy And Risk Reduction for Tsunamis in Europe. Grant 603839, 7th FP (ENV.2013.6.4-3 ENV

Exact Solution of a Faraday's Law Problem that Includes a Nonlinear Term and Its Implication for Perturbation Theory.

Science.gov (United States)

Fulcher, Lewis P.

1979-01-01

Presents an exact solution to the nonlinear Faraday's law problem of a rod sliding on frictionless rails with resistance. Compares the results with perturbation calculations based on the methods of Poisson and Pincare and of Kryloff and Bogoliuboff. (Author/GA)

FRF decoupling of nonlinear systems

Science.gov (United States)

Kalaycıoğlu, Taner; Özgüven, H. Nevzat

2018-03-01

Structural decoupling problem, i.e. predicting dynamic behavior of a particular substructure from the knowledge of the dynamics of the coupled structure and the other substructure, has been well investigated for three decades and led to several decoupling methods. In spite of the inherent nonlinearities in a structural system in various forms such as clearances, friction and nonlinear stiffness, all decoupling studies are for linear systems. In this study, decoupling problem for nonlinear systems is addressed for the first time. A method, named as FRF Decoupling Method for Nonlinear Systems (FDM-NS), is proposed for calculating FRFs of a substructure decoupled from a coupled nonlinear structure where nonlinearity can be modeled as a single nonlinear element. Depending on where nonlinear element is, i.e., either in the known or unknown subsystem, or at the connection point, the formulation differs. The method requires relative displacement information between two end points of the nonlinear element, in addition to point and transfer FRFs at some points of the known subsystem. However, it is not necessary to excite the system from the unknown subsystem even when the nonlinear element is in that subsystem. The validation of FDM-NS is demonstrated with two different case studies using nonlinear lumped parameter systems. Finally, a nonlinear experimental test structure is used in order to show the real-life application and accuracy of FDM-NS.

Existence of Positive Solutions to a Singular Semipositone Boundary Value Problem of Nonlinear Fractional Differential Systems

Directory of Open Access Journals (Sweden)

Xiaofeng Zhang

2017-12-01

Full Text Available In this paper, we consider the existence of positive solutions to a singular semipositone boundary value problem of nonlinear fractional differential equations. By applying the fixed point index theorem, some new results for the existence of positive solutions are obtained. In addition, an example is presented to demonstrate the application of our main results.

Nonlinear Dynamic Phenomena in Mechanics

CERN Document Server

Warminski, Jerzy; Cartmell, Matthew P

2012-01-01

Nonlinear phenomena should play a crucial role in the design and control of engineering systems and structures as they can drastically change the prevailing dynamical responses. This book covers theoretical and applications-based problems of nonlinear dynamics concerned with both discrete and continuous systems of interest in civil and mechanical engineering. They include pendulum-like systems, slender footbridges, shape memory alloys, sagged elastic cables and non-smooth problems. Pendulums can be used as a dynamic absorber mounted in high buildings, bridges or chimneys. Geometrical nonlinear

Implementation of the - Constraint Method in Special Class of Multi-objective Fuzzy Bi-Level Nonlinear Problems

Directory of Open Access Journals (Sweden)

Azza Hassan Amer

2017-12-01

Full Text Available Geometric programming problem is a powerful tool for solving some special type nonlinear programming problems. In the last few years we have seen a very rapid development on solving multiobjective geometric programming problem. A few mathematical programming methods namely fuzzy programming, goal programming and weighting methods have been applied in the recent past to find the compromise solution. In this paper, -constraint method has been applied in bi-level multiobjective geometric programming problem to find the Pareto optimal solution at each level. The equivalent mathematical programming problems are formulated to find their corresponding value of the objective function based on the duality theorem at eash level. Here, we have developed a new algorithm for fuzzy programming technique to solve bi-level multiobjective geometric programming problems to find an optimal compromise solution. Finally the solution procedure of the fuzzy technique is illustrated by a numerical example

Nonlinear problems in data-assimilation : Can synchronization help?

Science.gov (United States)

Tribbia, J. J.; Duane, G. S.

2009-12-01

Over the past several years, operational weather centers have initiated ensemble prediction and assimilation techniques to estimate the error covariance of forecasts in the short and the medium range. The ensemble techniques used are based on linear methods. The theory This technique s been shown to be a useful indicator of skill in the linear range where forecast errors are small relative to climatological variance. While this advance has been impressive, there are still ad hoc aspects of its use in practice, like the need for covariance inflation which are troubling. Furthermore, to be of utility in the nonlinear range an ensemble assimilation and prediction method must be capable of giving probabilistic information for the situation where a probability density forecast becomes multi-modal. A prototypical, simplest example of such a situation is the planetary-wave regime transition where the pdf is bimodal. Our recent research show how the inconsistencies and extensions of linear methodology can be consistently treated using the paradigm of synchronization which views the problems of assimilation and forecasting as that of optimizing the forecast model state with respect to the future evolution of the atmosphere.

Identification of Nonlinear Dynamic Systems Possessing Some Non-linearities

Directory of Open Access Journals (Sweden)

Y. N. Pavlov

2015-01-01

Full Text Available The subject of this work is the problem of identification of nonlinear dynamic systems based on the experimental data obtained by applying test signals to the system. The goal is to determinate coefficients of differential equations of systems by experimental frequency hodographs and separate similar, but different, in essence, forces: dissipative forces with the square of the first derivative in the motion equations and dissipative force from the action of dry friction. There was a proposal to use the harmonic linearization method to approximate each of the nonlinearity of "quadratic friction" and "dry friction" by linear friction with the appropriate harmonic linearization coefficient.Assume that a frequency transfer function of the identified system has a known form. Assume as well that there are disturbances while obtaining frequency characteristics of the realworld system. As a result, the points of experimentally obtained hodograph move randomly. Searching for solution of the identification problem was in the hodograph class, specified by the system model, which has the form of the frequency transfer function the same as the form of the frequency transfer function of the system identified. Minimizing a proximity criterion (measure of the experimentally obtained system hodograph and the system hodograph model for all the experimental points described and previously published by one of the authors allowed searching for the unknown coefficients of the frequenc ransfer function of the system model. The paper shows the possibility to identify a nonlinear dynamic system with multiple nonlinearities, obtained on the experimental samples of the frequency system hodograph. The proposed algorithm allows to select the nonlinearity of the type "quadratic friction" and "dry friction", i.e. also in the case where the nonlinearity is dependent on the same dynamic parameter, in particular, on the derivative of the system output value. For the dynamic

A numerical study of non-linear crack tip parameters

Directory of Open Access Journals (Sweden)

F.V. Antunes

2015-07-01

Full Text Available Crack closure concept has been widely used to explain different issues of fatigue crack propagation. However, different authors have questioned the relevance of crack closure and have proposed alternative concepts. The main objective here is to check the effectiveness of crack closure concept by linking the contact of crack flanks with non-linear crack tip parameters. Accordingly, 3D-FE numerical models with and without contact were developed for a wide range of loading scenarios and the crack tip parameters usually linked to fatigue crack growth, namely range of cyclic plastic strain, crack tip opening displacement, size of reversed plastic zone and total plastic dissipation per cycle, were investigated. It was demonstrated that: i LEFM concepts are applicable to the problem under study; ii the crack closure phenomenon has a great influence on crack tip parameters decreasing their values; iii the Keff concept is able to explain the variations of crack tip parameters produced by the contact of crack flanks; iv the analysis of remote compliance is the best numerical parameter to quantify the crack opening level; v without contact there is no effect of stress ratio on crack tip parameters. Therefore it is proved that the crack closure concept is valid.

Variational problems arising in classical mechanics and nonlinear elasticity

International Nuclear Information System (INIS)

Spencer, P.

1999-01-01

In this thesis we consider two different classes of variational problems. First, one-dimensional problems arising from classical mechanics where the problem is to determine whether there is a unique function η 0 (x) which minimises the energy functional of the form I(η) = ∫ a b L(x,η(x), η'(x)) dx. We will investigate uniqueness by making a change of dependent and independent variables and showing that for a class of integrands L with a particular kind of scaling invariance the resulting integrand is completely convex. The change of variables arises by applying results from Lie group theory as applied in the study of differential equations and this work is motivated by [60] and [68]. Second, the problem of minimising energy functionals of the form E(u) = ∫ A W(∇u(x)) dx in the case of a nonlinear elastic body occupying an annular region A contains R 2 with u : A-bar → A-bar. This work is motivated by [57] (in particular the example of paragraph 4). We will consider rotationally symmetric deformations satisfying prescribed boundary conditions. We will show the existence of minimisers for stored energy functions of the form W(F) = g-tilde(vertical bar-F-vertical bar, det(F)) in a class of general rotationally symmetric deformations of a compressible annulus and for stored energy functions of the form W(F) = g-bar(vertical bar-F-vertical bar) in a class of rotationally symmetric deformations of an incompressible annulus. We will also show that in each case the minimisers are solutions of the full equilibrium equations. A model problem will be considered where the energy functional is the Dirichlet integral and it will be shown that the rotationally symmetric solution obtained is a minimiser among admissible non-rotationally symmetric deformations. In the case of an incompressible annulus, we will consider the Dirichlet integral as the energy functional and show that the rotationally symmetric equilibrium solutions in this case are weak local minimisers in

Dynamics and vibrations progress in nonlinear analysis

CERN Document Server

Kachapi, Seyed Habibollah Hashemi

2014-01-01

Dynamical and vibratory systems are basically an application of mathematics and applied sciences to the solution of real world problems. Before being able to solve real world problems, it is necessary to carefully study dynamical and vibratory systems and solve all available problems in case of linear and nonlinear equations using analytical and numerical methods. It is of great importance to study nonlinearity in dynamics and vibration; because almost all applied processes act nonlinearly, and on the other hand, nonlinear analysis of complex systems is one of the most important and complicated tasks, especially in engineering and applied sciences problems. There are probably a handful of books on nonlinear dynamics and vibrations analysis. Some of these books are written at a fundamental level that may not meet ambitious engineering program requirements. Others are specialized in certain fields of oscillatory systems, including modeling and simulations. In this book, we attempt to strike a balance between th...

Mathematical problems in non-linear Physics: some results

International Nuclear Information System (INIS)

1979-01-01

The basic results presented in this report are the following: 1) Characterization of the range and Kernel of the variational derivative. 2) Determination of general conservation laws in linear evolution equations, as well as bounds for the number of polynomial conserved densities in non-linear evolution equations in two independent variables of even order. 3) Construction of the most general evolution equation which has a given family of conserved densities. 4) Regularity conditions for the validity of the Lie invariance method. 5) A simple class of perturbations in non-linear wave equations. 6) Soliton solutions in generalized KdV equations. (author)

First contact diagnosis and management of contact lens-related complications.

Science.gov (United States)

Fagan, Xavier J; Jhanji, Vishal; Constantinou, Marios; Amirul Islam, F M; Taylor, Hugh R; Vajpayee, Rasik B

2012-08-01

To describe the spectrum of contact lens-related problems in cases presenting to a tertiary referral eye hospital. A retrospective case record analysis of 111 eyes of 97 consecutive patients was undertaken over a period of five months at the Royal Victorian Eye and Ear Hospital, Melbourne, Australia. Contact lens-related complications (CLRC) were classified into microbial keratitis, sterile corneal infiltrates, corneal epitheliopathy and contact lens-related red eye (CLARE). Main parameters examined were nature of the first contact, clinical diagnosis, and management pattern. Forty-two percent of the initial presentations were to health care practitioners (HCPs) other than ophthalmologists. Mean duration from the onset of symptoms to presentation was 6.3 ± 10.9 days. Forty-nine percent (n = 54) of patients had an associated risk factor, most commonly overnight use of contact lenses (n = 14, 13 %). Most common diagnosis at presentation was corneal epitheliopathy (68 %) followed by sterile infiltrates (10 %), CLARE (8 %) and microbial keratitis (6 %). No significant differences were found in the pattern of treatment modalities administered by ophthalmologists and other HCPs. HCPs other than ophthalmologists are the first contact for contact lens-related problems in a significant proportion of patients. These HCPs manage the majority of CLRC by direct treatment or immediate referral.

H∞ Balancing for Nonlinear Systems

NARCIS (Netherlands)

Scherpen, Jacquelien M.A.

1996-01-01

In previously obtained balancing methods for nonlinear systems a past and a future energy function are used to bring the nonlinear system in balanced form. By considering a different pair of past and future energy functions that are related to the H∞ control problem for nonlinear systems we define

A HIGH ORDER SOLUTION OF THREE DIMENSIONAL TIME DEPENDENT NONLINEAR CONVECTIVE-DIFFUSIVE PROBLEM USING MODIFIED VARIATIONAL ITERATION METHOD

Directory of Open Access Journals (Sweden)

Pratibha Joshi

2014-12-01

Full Text Available In this paper, we have achieved high order solution of a three dimensional nonlinear diffusive-convective problem using modified variational iteration method. The efficiency of this approach has been shown by solving two examples. All computational work has been performed in MATHEMATICA.

Limits to Nonlinear Inversion

DEFF Research Database (Denmark)

Mosegaard, Klaus

2012-01-01

For non-linear inverse problems, the mathematical structure of the mapping from model parameters to data is usually unknown or partly unknown. Absence of information about the mathematical structure of this function prevents us from presenting an analytical solution, so our solution depends on our......-heuristics are inefficient for large-scale, non-linear inverse problems, and that the 'no-free-lunch' theorem holds. We discuss typical objections to the relevance of this theorem. A consequence of the no-free-lunch theorem is that algorithms adapted to the mathematical structure of the problem perform more efficiently than...... pure meta-heuristics. We study problem-adapted inversion algorithms that exploit the knowledge of the smoothness of the misfit function of the problem. Optimal sampling strategies exist for such problems, but many of these problems remain hard. © 2012 Springer-Verlag....

Implementation of a multi-layer perception for a non-linear control problem

International Nuclear Information System (INIS)

Lister, J.B.; Schnurrenberger, H.; Marmillod, P.

1990-12-01

We present the practical application of a 1-hidden-layer multilayer perception (MLP) to provide a non-linear continuous multi-variable mapping with 42 inputs and 13 outputs. The problem resolved is one of extracting information from experimental signals with a bandwidth of many kilohertz. We have an exact model of the inverse mapping of this problem, but we have no explicit form of the required forward mapping. This is the typical situation in data interpretation. The MLP was trained to provide this mapping by learning on 500 input-output examples. The success of the off-line solution to this problem using an MLP led us to examine the robustness of the MLP to different noise sources. We found that the MLP is more robust than an approximate linear mapping of the same problem. 12 bits of resolution in the weights are necessary to avoid a significant loss of precision. The practical implementation of large analog weight matrices using DAS-multipliers and a simple segmented sigmoid is also presented. A General Adaptive Recipe (GAR) for improving the performance of conventional back-propagation was developed. The GAR uses an adaptive step length and both the bias terms and the initial weight seeding are determined by the network size. The GAR was found to be robust and much faster than conventional back-propagation. (author) 20 figs., 1 tab., 15 refs

Contacts, non-linear transport effects and failure in multi-walled carbon nanotubes

International Nuclear Information System (INIS)

Berger, C; Yi, Y; Gezo, J; Poncharal, P; Heer, W A de

2003-01-01

Pristine arc-produced multi-walled carbon nanotubes are contacted to liquid mercury in situ in a transmission electron microscope. The conductance G(V) for all tubes increases with increasing bias voltage V. This is related to the electronic density of the nanotubes. Similar G(V) behaviour is observed for HOPG-graphite contacted in air with Hg, with dG(V)/dV∼0.3G 0 . Variations observed in the conductance are related to nanotube-Hg contact effects. For tubes barely touching the Hg surface, the conductance is low (typically G(V=0)∼0.1-0.5G 0 ); G(V) may maximize around V=1.5-2 V or continue to increase linearly depending on the MWNT-Hg contact. For good contacts the maximum low-bias conductance is 1G 0 . Non-conducting tubes are observed having a low-bias conductance smaller than 10 -3 G 0 . High-voltage tube failure usually occurs at the contact with Hg for clean tubes, or at tube defects. An important phenomenon is the formation of a Hg bubble near the contact nanotube-Hg surface when the nanotube is negatively biased, under high bias current conditions, indicating the heating effect of hot electrons injected into the mercury

Nonlinear Scattering from Partially Closed Cracks and Imperfect Interfaces

International Nuclear Information System (INIS)

Pecorari, Claudio

2004-05-01

This project has investigated the potential offered by nonlinear scattering phenomena to detect stress-corrosion, surface-breaking cracks, and regions of extended interfaces which are often invisible to conventional inspection methods because of their partial closure and/or the high background noise generated by the surrounding microstructure. The investigation has looked into the basic physics of the interaction between ultrasonic waves and rough surfaces in contact, since the latter offers a prototypical example of a mechanical system which is characterized by a dynamics similar to that of a partially closed crack. To this end, three fundamental mechanisms which may be activated by an inspecting ultrasonic wave have been considered. The first mechanism is described by the Hertz force law which governs the interaction between asperities in contact that are subjected to a normal load. The second mechanism considers the dynamics of two spherical asperities subjected to an oscillating tangential load. To this end, the model developed by Mindlin and Deresiewizc has been used. The third mechanism accounts for the effect of forces of adhesion, and can be described by a model developed by Greenwood and Johnson. The validity of this model is rather general and covers the extreme cases of very soft and very rigid contacts. This model aims at describing the effect of fluid layers with thickness of atomic size, which may be present within a crack. Statistical models accounting for the topography of the two rough surfaces in contact have been developed, and the macroscopic stiffness of the interface recovered. These results have been used to formulate effective boundary conditions to be enforced at the interface, and the reflection and transmission problem has been solved in a variety of situations of experimental significance. The main conclusion of this part of the project is that the second harmonic component is the dominant feature of the nonlinear response of an

Simulation of High-current Pulse Effect on the Electrode with Nonlinear Material Characteristics and Phase Transitions Taken into Account

Directory of Open Access Journals (Sweden)

R. V. Arutjunjan

2016-01-01

Full Text Available The article investigates the thermal and electrical processes when heating the metal electrode by a high current pulse. The aim is to understand an impact nature of the nonlinearities of thermal parameters, the phase transitions of melting and evaporation, and the type of boundary conditions in the current spot. To solve the problem was formulated a mathematical model, and were also developed a finite-difference method and computer programmes which allow an effective computer simulations of thermal and electrical processes under the high current pulse impact on the metal electrodes. The Stefan problem is solved by the through "enthalpy" method. Calculation of the electric field is performed by Seidel iteration. Thermal and current balance and comparison with solution results of model problems allow computer error monitoring.The work involved a series of calculations for an informative case of iron. It enabled to find a significant influence of the nonlinearities of thermal parameters, the phase transitions of melting and evaporation, the type of boundary conditions on the values of the temperature and electric fields, especially in the vicinity of the current spot. The presence of high current density and temperature, respectively, in the vicinity of the current spot edge confirms the well-known hypothesis about the causes of contact welding on the edges of the contact area. It has been found that the impact of losses on radiation and convection cooling is negligible. The article continues and complements the well-known research in the theory of electrical contacts and welding processes based on detailed consideration of the electrode material properties, the nonlinearities, and a type of boundary conditions for temperature and electric fields.The results can be used in the practice in research and design of electrical machines and other electrical devices.The study has revealed the need to improve the enthalpy finite- difference method for

Non-conformal contact mechanical characteristic analysis on spherical components

Energy Technology Data Exchange (ETDEWEB)

Zhen-zhi, G.; Bin, H.; Zheng-ming, G.; Feng-mei, Y.; Jin, Q [The 2. Artillery Engineering Univ., Xi' an (China)

2017-03-15

Non-conformal spherical-contact mechanical problems is a three-dimensional coordination or similar to the coordination spherical contact. Due to the complexity of the problem of spherical-contact and difficulties of solving higher-order partial differential equations, problems of three-dimensional coordination or similar to the coordination spherical-contact is still no exact analytical method for solving. It is based on three-dimensional taper model is proposed a model based on the contour surface of the spherical contact and concluded of the formula of the contact pressure and constructed of finite element model by contact pressure distribution under the non-conformal spherical. The results shows spherical contact model can reflect non-conformal spherical-contacting mechanical problems more than taper-contacting model, and apply for the actual project.

An elastic-plastic contact model for line contact structures

Science.gov (United States)

Zhu, Haibin; Zhao, Yingtao; He, Zhifeng; Zhang, Ruinan; Ma, Shaopeng

2018-06-01

Although numerical simulation tools are now very powerful, the development of analytical models is very important for the prediction of the mechanical behaviour of line contact structures for deeply understanding contact problems and engineering applications. For the line contact structures widely used in the engineering field, few analytical models are available for predicting the mechanical behaviour when the structures deform plastically, as the classic Hertz's theory would be invalid. Thus, the present study proposed an elastic-plastic model for line contact structures based on the understanding of the yield mechanism. A mathematical expression describing the global relationship between load history and contact width evolution of line contact structures was obtained. The proposed model was verified through an actual line contact test and a corresponding numerical simulation. The results confirmed that this model can be used to accurately predict the elastic-plastic mechanical behaviour of a line contact structure.

Global solutions to the initial-boundary value problem for the quasilinear viscoelastic equation with a derivative nonlinearity

Directory of Open Access Journals (Sweden)

Mitsuhiro Nakao

2014-01-01

Full Text Available We prove the existence and uniqueness of a global decaying solution to the initial boundary value problem for the quasilinear wave equation with Kelvin-Voigt dissipation and a derivative nonlinearity. To derive the required estimates of the solutions we employ a 'loan' method and use a difference inequality on the energy.

New results on the mathematical problems in nonlinear physics

International Nuclear Information System (INIS)

1980-01-01

The main topics treated in this report are: I) Existence of generalized Lagrangians. II) Conserved densities for odd-order polynomial evolution equations and linear evolution systems. III ) Conservation laws for Klein-Gordon, Di rae and Maxwell equations. IV) Stability conditions for finite-energy solutions of a non-linear Klein-Gordon equation. V) Hamiltonian approach to non-linear evolution equations and Backlund transformations. VI) Anharmonic vibrations: Status of results and new possible approaches. (Author) 83 refs

Averaging of nonlinearity-managed pulses

International Nuclear Information System (INIS)

Zharnitsky, Vadim; Pelinovsky, Dmitry

2005-01-01

We consider the nonlinear Schroedinger equation with the nonlinearity management which describes Bose-Einstein condensates under Feshbach resonance. By using an averaging theory, we derive the Hamiltonian averaged equation and compare it with other averaging methods developed for this problem. The averaged equation is used for analytical approximations of nonlinearity-managed solitons

Degenerate nonlinear diffusion equations

CERN Document Server

Favini, Angelo

2012-01-01

The aim of these notes is to include in a uniform presentation style several topics related to the theory of degenerate nonlinear diffusion equations, treated in the mathematical framework of evolution equations with multivalued m-accretive operators in Hilbert spaces. The problems concern nonlinear parabolic equations involving two cases of degeneracy. More precisely, one case is due to the vanishing of the time derivative coefficient and the other is provided by the vanishing of the diffusion coefficient on subsets of positive measure of the domain. From the mathematical point of view the results presented in these notes can be considered as general results in the theory of degenerate nonlinear diffusion equations. However, this work does not seek to present an exhaustive study of degenerate diffusion equations, but rather to emphasize some rigorous and efficient techniques for approaching various problems involving degenerate nonlinear diffusion equations, such as well-posedness, periodic solutions, asympt...

MINPACK-1, Subroutine Library for Nonlinear Equation System

International Nuclear Information System (INIS)

Garbow, Burton S.

1984-01-01

1 - Description of problem or function: MINPACK1 is a package of FORTRAN subprograms for the numerical solution of systems of non- linear equations and nonlinear least-squares problems. The individual programs are: Identification/Description: - CHKDER: Check gradients for consistency with functions, - DOGLEG: Determine combination of Gauss-Newton and gradient directions, - DPMPAR: Provide double precision machine parameters, - ENORM: Calculate Euclidean norm of vector, - FDJAC1: Calculate difference approximation to Jacobian (nonlinear equations), - FDJAC2: Calculate difference approximation to Jacobian (least squares), - HYBRD: Solve system of nonlinear equations (approximate Jacobian), - HYBRD1: Easy-to-use driver for HYBRD, - HYBRJ: Solve system of nonlinear equations (analytic Jacobian), - HYBRJ1: Easy-to-use driver for HYBRJ, - LMDER: Solve nonlinear least squares problem (analytic Jacobian), - LMDER1: Easy-to-use driver for LMDER, - LMDIF: Solve nonlinear least squares problem (approximate Jacobian), - LMDIF1: Easy-to-use driver for LMDIF, - LMPAR: Determine Levenberg-Marquardt parameter - LMSTR: Solve nonlinear least squares problem (analytic Jacobian, storage conserving), - LMSTR1: Easy-to-use driver for LMSTR, - QFORM: Accumulate orthogonal matrix from QR factorization QRFAC Compute QR factorization of rectangular matrix, - QRSOLV: Complete solution of least squares problem, - RWUPDT: Update QR factorization after row addition, - R1MPYQ: Apply orthogonal transformations from QR factorization, - R1UPDT: Update QR factorization after rank-1 addition, - SPMPAR: Provide single precision machine parameters. 4. Method of solution - MINPACK1 uses the modified Powell hybrid method and the Levenberg-Marquardt algorithm

Application of the Green's function method to some nonlinear problems of an electron storage ring

International Nuclear Information System (INIS)

Kheifets, S.

1984-01-01

One of the most important characteristics of an electron storage ring is the size of the beam. However analytical calculations of beam size are beset with problems and the computational methods and programs which are used to overcome these are inadequate for all problems in which stochastic noise is an essential part. Two examples are, for an electron storage ring, beam-size evaluation including beam-beam interactions, and finding the beam size for a nonlinear machine. The method described should overcome some of the problems. It uses the Green's function method applied to the Fokker-Planck equation governing the distribution function in the phase space of particle motion. The new step is to consider the particle motion in two degrees of freedom rather than in one dimension. The technique is described fully and is then applied to a strong-focusing machine. (U.K.)

Convergence rates and finite-dimensional approximations for nonlinear ill-posed problems involving monotone operators in Banach spaces

International Nuclear Information System (INIS)

Nguyen Buong.

1992-11-01

The purpose of this paper is to investigate convergence rates for an operator version of Tikhonov regularization constructed by dual mapping for nonlinear ill-posed problems involving monotone operators in real reflective Banach spaces. The obtained results are considered in combination with finite-dimensional approximations for the space. An example is considered for illustration. (author). 15 refs

A nonlinear plate control without linearization

Directory of Open Access Journals (Sweden)

Yildirim Kenan

2017-03-01

Full Text Available In this paper, an optimal vibration control problem for a nonlinear plate is considered. In order to obtain the optimal control function, wellposedness and controllability of the nonlinear system is investigated. The performance index functional of the system, to be minimized by minimum level of control, is chosen as the sum of the quadratic 10 functional of the displacement. The velocity of the plate and quadratic functional of the control function is added to the performance index functional as a penalty term. By using a maximum principle, the nonlinear control problem is transformed to solving a system of partial differential equations including state and adjoint variables linked by initial-boundary-terminal conditions. Hence, it is shown that optimal control of the nonlinear systems can be obtained without linearization of the nonlinear term and optimal control function can be obtained analytically for nonlinear systems without linearization.

Real-time simulation of contact and cutting of heterogeneous soft-tissues.

Science.gov (United States)

Courtecuisse, Hadrien; Allard, Jérémie; Kerfriden, Pierre; Bordas, Stéphane P A; Cotin, Stéphane; Duriez, Christian

2014-02-01

This paper presents a numerical method for interactive (real-time) simulations, which considerably improves the accuracy of the response of heterogeneous soft-tissue models undergoing contact, cutting and other topological changes. We provide an integrated methodology able to deal both with the ill-conditioning issues associated with material heterogeneities, contact boundary conditions which are one of the main sources of inaccuracies, and cutting which is one of the most challenging issues in interactive simulations. Our approach is based on an implicit time integration of a non-linear finite element model. To enable real-time computations, we propose a new preconditioning technique, based on an asynchronous update at low frequency. The preconditioner is not only used to improve the computation of the deformation of the tissues, but also to simulate the contact response of homogeneous and heterogeneous bodies with the same accuracy. We also address the problem of cutting the heterogeneous structures and propose a method to update the preconditioner according to the topological modifications. Finally, we apply our approach to three challenging demonstrators: (i) a simulation of cataract surgery (ii) a simulation of laparoscopic hepatectomy (iii) a brain tumor surgery. Copyright © 2013 Elsevier B.V. All rights reserved.

Recent developments of some asymptotic methods and their applications for nonlinear vibration equations in engineering problems: a review

Directory of Open Access Journals (Sweden)

Mahmoud Bayat

Full Text Available This review features a survey of some recent developments in asymptotic techniques and new developments, which are valid not only for weakly nonlinear equations, but also for strongly ones. Further, the achieved approximate analytical solutions are valid for the whole solution domain. The limitations of traditional perturbation methods are illustrated, various modified perturbation techniques are proposed, and some mathematical tools such as variational theory, homotopy technology, and iteration technique are introduced to over-come the shortcomings.In this review we have applied different powerful analytical methods to solve high nonlinear problems in engineering vibrations. Some patterns are given to illustrate the effectiveness and convenience of the methodologies.

Explicit Nonlinear Model Predictive Control Theory and Applications

CERN Document Server

Grancharova, Alexandra

2012-01-01

Nonlinear Model Predictive Control (NMPC) has become the accepted methodology to solve complex control problems related to process industries. The main motivation behind explicit NMPC is that an explicit state feedback law avoids the need for executing a numerical optimization algorithm in real time. The benefits of an explicit solution, in addition to the efficient on-line computations, include also verifiability of the implementation and the possibility to design embedded control systems with low software and hardware complexity. This book considers the multi-parametric Nonlinear Programming (mp-NLP) approaches to explicit approximate NMPC of constrained nonlinear systems, developed by the authors, as well as their applications to various NMPC problem formulations and several case studies. The following types of nonlinear systems are considered, resulting in different NMPC problem formulations: Ø  Nonlinear systems described by first-principles models and nonlinear systems described by black-box models; �...

A chaos-based evolutionary algorithm for general nonlinear programming problems

International Nuclear Information System (INIS)

El-Shorbagy, M.A.; Mousa, A.A.; Nasr, S.M.

2016-01-01

In this paper we present a chaos-based evolutionary algorithm (EA) for solving nonlinear programming problems named chaotic genetic algorithm (CGA). CGA integrates genetic algorithm (GA) and chaotic local search (CLS) strategy to accelerate the optimum seeking operation and to speed the convergence to the global solution. The integration of global search represented in genetic algorithm and CLS procedures should offer the advantages of both optimization methods while offsetting their disadvantages. By this way, it is intended to enhance the global convergence and to prevent to stick on a local solution. The inherent characteristics of chaos can enhance optimization algorithms by enabling it to escape from local solutions and increase the convergence to reach to the global solution. Twelve chaotic maps have been analyzed in the proposed approach. The simulation results using the set of CEC’2005 show that the application of chaotic mapping may be an effective strategy to improve the performances of EAs.

An overview of adaptive model theory: solving the problems of redundancy, resources, and nonlinear interactions in human movement control.

Science.gov (United States)

Neilson, Peter D; Neilson, Megan D

2005-09-01

Adaptive model theory (AMT) is a computational theory that addresses the difficult control problem posed by the musculoskeletal system in interaction with the environment. It proposes that the nervous system creates motor maps and task-dependent synergies to solve the problems of redundancy and limited central resources. These lead to the adaptive formation of task-dependent feedback/feedforward controllers able to generate stable, noninteractive control and render nonlinear interactions unobservable in sensory-motor relationships. AMT offers a unified account of how the nervous system might achieve these solutions by forming internal models. This is presented as the design of a simulator consisting of neural adaptive filters based on cerebellar circuitry. It incorporates a new network module that adaptively models (in real time) nonlinear relationships between inputs with changing and uncertain spectral and amplitude probability density functions as is the case for sensory and motor signals.

Numerical nonlinear complex geometrical optics algorithm for the 3D Calderón problem

DEFF Research Database (Denmark)

Delbary, Fabrice; Knudsen, Kim

2014-01-01

to the generalized Laplace equation. The 3D problem was solved in theory in late 1980s using complex geometrical optics solutions and a scattering transform. Several approximations to the reconstruction method have been suggested and implemented numerically in the literature, but here, for the first time, a complete...... computer implementation of the full nonlinear algorithm is given. First a boundary integral equation is solved by a Nystrom method for the traces of the complex geometrical optics solutions, second the scattering transform is computed and inverted using fast Fourier transform, and finally a boundary value...

Contact Problem for an Elastic Layer on an Elastic Half Plane Loaded by Means of Three Rigid Flat Punches

Directory of Open Access Journals (Sweden)

T. S. Ozsahin

2013-01-01

Full Text Available The frictionless contact problem for an elastic layer resting on an elastic half plane is considered. The problem is solved by using the theory of elasticity and integral transformation technique. The compressive loads P and Q (per unit thickness in direction are applied to the layer through three rigid flat punches. The elastic layer is also subjected to uniform vertical body force due to effect of gravity. The contact along the interface between elastic layer and half plane is continuous, if the value of the load factor, λ, is less than a critical value, . In this case, initial separation loads, and initial separation points, are determined. Also the required distance between the punches to avoid any separation between the punches and the elastic layer is studied and the limit distance between punches that ends interaction of punches is investigated for various dimensionless quantities. However, if tensile tractions are not allowed on the interface, for the layer separates from the interface along a certain finite region. Numerical results for distance determining the separation area, vertical displacement in the separation zone, contact stress distribution along the interface between elastic layer and half plane are given for this discontinuous contact case.

The forced nonlinear Schroedinger equation

International Nuclear Information System (INIS)

Kaup, D.J.; Hansen, P.J.

1985-01-01

The nonlinear Schroedinger equation describes the behaviour of a radio frequency wave in the ionosphere near the reflexion point where nonlinear processes are important. A simple model of this phenomenon leads to the forced nonlinear Schroedinger equation in terms of a nonlinear boundary value problem. A WKB analysis of the time evolution equations for the nonlinear Schroedinger equation in the inverse scattering transform formalism gives a crude order of magnitude estimation of the qualitative behaviour of the solutions. This estimation is compared with the numerical solutions. (D.Gy.)

Relation of deformed nonlinear algebras with linear ones

International Nuclear Information System (INIS)

Nowicki, A; Tkachuk, V M

2014-01-01

The relation between nonlinear algebras and linear ones is established. For a one-dimensional nonlinear deformed Heisenberg algebra with two operators we find the function of deformation for which this nonlinear algebra can be transformed to a linear one with three operators. We also establish the relation between the Lie algebra of total angular momentum and corresponding nonlinear one. This relation gives a possibility to simplify and to solve the eigenvalue problem for the Hamiltonian in a nonlinear case using the reduction of this problem to the case of linear algebra. It is demonstrated in an example of a harmonic oscillator. (paper)

Inverse atmospheric radiative transfer problems - A nonlinear minimization search method of solution. [aerosol pollution monitoring

Science.gov (United States)

Fymat, A. L.

1976-01-01

The paper studies the inversion of the radiative transfer equation describing the interaction of electromagnetic radiation with atmospheric aerosols. The interaction can be considered as the propagation in the aerosol medium of two light beams: the direct beam in the line-of-sight attenuated by absorption and scattering, and the diffuse beam arising from scattering into the viewing direction, which propagates more or less in random fashion. The latter beam has single scattering and multiple scattering contributions. In the former case and for single scattering, the problem is reducible to first-kind Fredholm equations, while for multiple scattering it is necessary to invert partial integrodifferential equations. A nonlinear minimization search method, applicable to the solution of both types of problems has been developed, and is applied here to the problem of monitoring aerosol pollution, namely the complex refractive index and size distribution of aerosol particles.

Contact Allergy in Danish Healthcare Workers

DEFF Research Database (Denmark)

Schwensen, Jakob F; Menné, Torkil; Sommerlund, Mette

2016-01-01

Contact dermatitis in healthcare workers is a pan-European problem. We conducted a retrospective observational study of the patch-test results of 1402 healthcare workers and 1402 matched controls with contact dermatitis who were treated at 3 hospitals departments in Denmark between 2007 and 2014....... The primary objective was to determine whether healthcare work was associated with contact allergy to thiuram mix. Unadjusted univariate analyses revealed that healthcare work was significantly associated with occupational contact dermatitis and hand dermatitis. Contact allergy to thiuram mix was more common...... in healthcare workers was significantly associated with having occupational contact dermatitis, hand dermatitis and older age. In conclusion, we report here a potential problem of contact allergy to thiurams in healthcare workers with contact dermatitis. Legislative authorities may in the future focus...

On a Solvability of Contact Problems with Visco-Plastic Friction in the Thermo-Visco-Plastic Bingham Rheology

Czech Academy of Sciences Publication Activity Database

Nedoma, Jiří

2006-01-01

Roč. 22, č. 4 (2006), s. 484-499 ISSN 0167-739X Institutional research plan: CEZ:AV0Z10300504 Keywords : unilateral contact problem * local visco- plastic friction * thermo-visco- plastic Bingham rheology * FEM Subject RIV: BA - General Mathematics Impact factor: 0.722, year: 2006

Modelling and solution of contact problem for infinite plate and cross-shaped embedment

Directory of Open Access Journals (Sweden)

O.B. Kozin

2016-09-01

Full Text Available Development of efficient methods of determination of an intense-strained state of thin-walled constructional designs with inclusions, reinforcements and other stress raisers is an important problem both with theoretical, and from the practical point of view, considering their wide practical application. Aim: The aim of this research is to develop the analytical mathematical method of studying of an intense-strained state of infinite plate with cross-shaped embedment at a bend. Materials and Methods: The method of boundary elements is an efficient way of the boundary value problems solution for systems of differential equations. The methods based on boundary integral equations get wide application in many branches of science and technique, calculation of plates and shells. One of methods of solution of a numerous class of the integral equations and systems arising on the basis of a method of boundary integral equations is the analytical method of construction of these equations and systems to Riemann problems with their forthcoming decision. Results: The integral equation for the analysis of deflections and the analysis of an intense-strained state of a thin rigid plate with rigid cross-shaped embedment is received. The precise solution of this boundary value problem is received by reduction to a Riemann problem and its forthcoming solution. An asymptotical behavior of contact efforts at the ends of embedment is investigated.

Nonlinear and stochastic dynamics of coherent structures

DEFF Research Database (Denmark)

Rasmussen, Kim

1997-01-01

This Thesis deals with nonlinear and stochastic dynamics in systems which can be described by nonlinear Schrödinger models. Basically three different models are investigated. The first is the continuum nonlinear Schröndinger model in one and two dimensions generalized by a tunable degree of nonli......This Thesis deals with nonlinear and stochastic dynamics in systems which can be described by nonlinear Schrödinger models. Basically three different models are investigated. The first is the continuum nonlinear Schröndinger model in one and two dimensions generalized by a tunable degree...... introduces the nonlinear Schrödinger model in one and two dimensions, discussing the soliton solutions in one dimension and the collapse phenomenon in two dimensions. Also various analytical methods are described. Then a derivation of the nonlinear Schrödinger equation is given, based on a Davydov like...... system described by a tight-binding Hamiltonian and a harmonic lattice coupled b y a deformation-type potential. This derivation results in a two-dimensional nonline ar Schrödinger model, and considering the harmonic lattice to be in thermal contact with a heat bath w e show that the nonlinear...

Inexact Newton–Landweber iteration for solving nonlinear inverse problems in Banach spaces

International Nuclear Information System (INIS)

Jin, Qinian

2012-01-01

By making use of duality mappings, we formulate an inexact Newton–Landweber iteration method for solving nonlinear inverse problems in Banach spaces. The method consists of two components: an outer Newton iteration and an inner scheme providing the increments by applying the Landweber iteration in Banach spaces to the local linearized equations. It has the advantage of reducing computational work by computing more cheap steps in each inner scheme. We first prove a convergence result for the exact data case. When the data are given approximately, we terminate the method by a discrepancy principle and obtain a weak convergence result. Finally, we test the method by reporting some numerical simulations concerning the sparsity recovery and the noisy data containing outliers. (paper)

Eigenvalue problem and nonlinear evolution of kink modes in a toroidal plasma

International Nuclear Information System (INIS)

Ogino, T.; Takeda, S.; Sanuki, H.; Kamimura, T.

1979-04-01

The internal kink modes of a cylindrical plasma are investigated by a linear eigen value problem and their nonlinear evolution is studied by 3-dimensional MHD simulation based on the rectangular column model under the fixed boundary condition. The growth rates in two cases, namely uniform and diffused current profiles are analyzed in detail, which agree with the analytical estimation by Shafranov. The time evolution of the m = 1 mode showed that q > 1 is satisfied at the relaxation time (q safety factor), a stable configuration like a horse shoe (a new equilibrium) was formed. Also, the time evolution of the pressure p for the m = 2 mode showed that a stable configuration like a pair of anchors was formed. (author)

Nonlinear analysis

CERN Document Server

Gasinski, Leszek

2005-01-01

Hausdorff Measures and Capacity. Lebesgue-Bochner and Sobolev Spaces. Nonlinear Operators and Young Measures. Smooth and Nonsmooth Analysis and Variational Principles. Critical Point Theory. Eigenvalue Problems and Maximum Principles. Fixed Point Theory.

Applications and algorithms for mixed integer nonlinear programming

International Nuclear Information System (INIS)

Leyffer, Sven; Munson, Todd; Linderoth, Jeff; Luedtke, James; Miller, Andrew

2009-01-01

The mathematical modeling of systems often requires the use of both nonlinear and discrete components. Discrete decision variables model dichotomies, discontinuities, and general logical relationships. Nonlinear functions are required to accurately represent physical properties such as pressure, stress, temperature, and equilibrium. Problems involving both discrete variables and nonlinear constraint functions are known as mixed-integer nonlinear programs (MINLPs) and are among the most challenging computational optimization problems faced by researchers and practitioners. In this paper, we describe relevant scientific applications that are naturally modeled as MINLPs, we provide an overview of available algorithms and software, and we describe ongoing methodological advances for solving MINLPs. These algorithmic advances are making increasingly larger instances of this important family of problems tractable.

Analysis of the components of a pipeline closure as a unilateral confact problem

International Nuclear Information System (INIS)

Barbosa, H.J.C.; Feijoo, R.A.

1987-03-01

In many situations in the stress analysis of piping systems and pressure vessels the need to deal with (non-linear) unilateral contact problems arises. A typical case is the analysis of the components of pipe closure. In this work, numerical algorithms developed in previous publications are apllied in an example where the system composed by the closure neck, the blanking plug and the locking ring of a pipe closure is analysed. (Author) [pt

Contact resistance problems applying ERT on low bulk density forested stony soils. Is there a solution?

Science.gov (United States)

Deraedt, Deborah; Touzé, Camille; Robert, Tanguy; Colinet, Gilles; Degré, Aurore; Garré, Sarah

2015-04-01

Electrical resistivity tomography (ERT) has often been put forward as a promising tool to quantify soil water and solute fluxes in a non-invasive way. In our experiment, we wanted to determine preferential flow processes along a forested hillslope using a saline tracer with ERT. The experiment was conducted in the Houille watershed, subcatchment of the Meuse located in the North of Belgian Ardennes (50° 1'52.6'N, 4° 53'22.5'E). The climate is continental but the soil under spruce (Picea abies (L.) Karst.) and Douglas fire stand (Pseudotsuga menziesii (Mirb.) Franco) remains quite dry (19% WVC in average) during the whole year. The soil is Cambisol and the parent rock is Devonian schist covered with variable thickness of silty loam soil. The soil density ranges from 1.13 to 1.87 g/cm3 on average. The stone content varies from 20 to 89% and the soil depth fluctuates between 70 and 130 cm. The ERT tests took place on June 1st 2012, April 1st, 2nd and 3rd 2014 and May 12th 2014. We used the Terrameter LS 12 channels (ABEM, Sweden) in 2012 test and the DAS-1 (Multi-Phase Technologies, United States) in 2014. Different electrode configurations and arrays were adopted for different dates (transect and grid arrays and Wenner - Schlumberger, Wenner alpha and dipole-dipole configurations). During all tests, we systematically faced technical problems, mainly related to bad electrode contact. The recorded data show values of contact resistance above 14873 Ω (our target value would be below 3000 Ω). Subsequently, we tried to improve the contact by predrilling the soil and pouring water in the electrode holes. The contact resistance improved to 14040 Ω as minimum. The same procedure with liquid mud was then tested to prevent quick percolation of the water from the electrode location. As a result, the lower contact resistance dropped to 11745 Ω. Finally, we applied about 25 litre of saline solution (CaCl2, 0.75g/L) homogeneously on the electrode grid. The minimum value of

On Inverse Coefficient Heat-Conduction Problems on Reconstruction of Nonlinear Components of the Thermal-Conductivity Tensor of Anisotropic Bodies

Science.gov (United States)

Formalev, V. F.; Kolesnik, S. A.

2017-11-01

The authors are the first to present a closed procedure for numerical solution of inverse coefficient problems of heat conduction in anisotropic materials used as heat-shielding ones in rocket and space equipment. The reconstructed components of the thermal-conductivity tensor depend on temperature (are nonlinear). The procedure includes the formation of experimental data, the implicit gradient-descent method, the economical absolutely stable method of numerical solution of parabolic problems containing mixed derivatives, the parametric identification, construction, and numerical solution of the problem for elements of sensitivity matrices, the development of a quadratic residual functional and regularizing functionals, and also the development of algorithms and software systems. The implicit gradient-descent method permits expanding the quadratic functional in a Taylor series with retention of the linear terms for the increments of the sought functions. This substantially improves the exactness and stability of solution of the inverse problems. Software systems are developed with account taken of the errors in experimental data and disregarding them. On the basis of a priori assumptions of the qualitative behavior of the functional dependences of the components of the thermal-conductivity tensor on temperature, regularizing functionals are constructed by means of which one can reconstruct the components of the thermal-conductivity tensor with an error no higher than the error of the experimental data. Results of the numerical solution of the inverse coefficient problems on reconstruction of nonlinear components of the thermal-conductivity tensor have been obtained and are discussed.

Nonlinear Approaches in Engineering Applications

CERN Document Server

Jazar, Reza

2012-01-01

Nonlinear Approaches in Engineering Applications focuses on nonlinear phenomena that are common in the engineering field. The nonlinear approaches described in this book provide a sound theoretical base and practical tools to design and analyze engineering systems with high efficiency and accuracy and with less energy and downtime. Presented here are nonlinear approaches in areas such as dynamic systems, optimal control and approaches in nonlinear dynamics and acoustics. Coverage encompasses a wide range of applications and fields including mathematical modeling and nonlinear behavior as applied to microresonators, nanotechnologies, nonlinear behavior in soil erosion,nonlinear population dynamics, and optimization in reducing vibration and noise as well as vibration in triple-walled carbon nanotubes. This book also: Provides a complete introduction to nonlinear behavior of systems and the advantages of nonlinearity as a tool for solving engineering problems Includes applications and examples drawn from the el...

Nonlinear dynamics of intense EM pulses in plasma

International Nuclear Information System (INIS)

Mahajan, Ranju; Gill, Tarsem Singh; Kaur, Ravinder

2010-01-01

The evolution of laser beam in underdense/overdense plasma medium which is key to understanding of several nonlinear processes and underlying physics is governed by nonlinear parabolic equation. The nonlinearity considered here is of relativistic as well as of ponderomotive type. We have set Lagrangian for the problem and reduced Lagrangian problem is solved using appropriate trial function. Equation for the beam width and phase are derived. Further, these equations are used to solve eigenvalue problem for the stability of laser beam evolution and Hurwitz condition is satisfied.

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

Different methods of solution of linear and nonlinear algebraic systems are applied to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems, methods in general use of alternating directions type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method on nonlinear conjugate gradient is studied as also Newton's method and some of its variants. It should be noted, however that Newton's method is found to be more efficient when coupled with a good method for solution of the linear system. To conclude, such methods are used to solve a nonlinear diffusion problem and the numerical results obtained are to be compared [fr

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

The object of this study is to compare different methods of solving linear and nonlinear algebraic systems and to apply them to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems the conventional methods of alternating direction type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method of nonlinear conjugate gradient is studied together with Newton's method and some of its variants. It should be noted, however, that Newton's method is found to be more efficient when coupled with a good method for solving the linear system. As a conclusion, these methods are used to solve a nonlinear diffusion problem and the numerical results obtained are compared [fr

Geometric nonlinear effects on the planar dynamics of a pivoted flexible beam encountering a point-surface impact

International Nuclear Information System (INIS)

Li Qing; Wang Tianshu; Ma Xingrui

2009-01-01

Flexible-body modeling with geometric nonlinearities remains a hot topic of research by applications in multibody system dynamics undergoing large overall motions. However, the geometric nonlinear effects on the impact dynamics of flexible multibody systems have attracted significantly less attention. In this paper, a point-surface impact problem between a rigid ball and a pivoted flexible beam is investigated. The Hertzian contact law is used to describe the impact process, and the dynamic equations are formulated in the floating frame of reference using the assumed mode method. The two important geometric nonlinear effects of the flexible beam are taken into account, i.e., the longitudinal foreshortening effect due to the transverse deformation, and the stress stiffness effect due to the axial force. The simulation results show that good consistency can be obtained with the nonlinear finite element program ABAQUS/Explicit if proper geometric nonlinearities are included in the floating frame formulation. Specifically, only the foreshortening effect should be considered in a pure transverse impact for efficiency, while the stress stiffness effect should be further considered in an oblique case with much more computational effort. It also implies that the geometric nonlinear effects should be considered properly in the impact dynamic analysis of more general flexible multibody systems

Flutter analysis of an airfoil with multiple nonlinearities and uncertainties

Directory of Open Access Journals (Sweden)

Haitao Liao

2013-09-01

Full Text Available An original method for calculating the limit cycle oscillations of nonlinear aero-elastic system is presented. The problem of determining the maximum vibration amplitude of limit cycle is transformed into a nonlinear optimization problem. The harmonic balance method and the Floquet theory are selected to construct the general nonlinear equality and inequality constraints. The resulting constrained maximization problem is then solved by using the MultiStart algorithm. Finally, the proposed approach is validated and used to analyse the limit cycle oscillations of an airfoil with multiple nonlinearities and uncertainties. Numerical examples show that the coexistence of multiple nonlinearities may lead to low amplitude limit cycle oscillation.

Colored Contact Lens Dangers

Science.gov (United States)

... Hazard: The Hidden Dangers of Buying Decorative Contact Lenses Without a Prescription Leer en Español: Peligros asociados ... truth." Real People, Real Problems with Colored Contact Lenses Julian: Teenager Blinded In One Eye By Non- ...

Extreme control of impulse transmission by cylinder-based nonlinear phononic crystals

Science.gov (United States)

Chaunsali, Rajesh; Toles, Matthew; Yang, Jinkyu; Kim, Eunho

2017-10-01

We present a novel device that can offer two extremes of elastic wave propagation - nearly complete transmission and strong attenuation under impulse excitation. The mechanism of this highly tunable device relies on intermixing effects of dispersion and nonlinearity. The device consists of identical cylinders arranged in a chain, which interact with each other as per nonlinear Hertz contact law. For a 'dimer' configuration, i.e., two different contact angles alternating in the chain, we analytically, numerically, and experimentally show that impulse excitation can either propagate as a localized wave, or it can travel as a highly dispersive wave. Remarkably, these extremes can be achieved in this periodic arrangement simply by in-situ control of contact angles between cylinders. We close the discussion by highlighting the key characteristics of the mechanisms that facilitate strong attenuation of incident impulse. These include low-to-high frequency scattering, and turbulence-like cascading in a periodic system. We thus envision that these adaptive, cylinder-based nonlinear phononic crystals, in conjunction with conventional impact mitigation mechanisms, could be used to design highly tunable and efficient impact manipulation devices.

Some nonlinear space decomposition algorithms

Energy Technology Data Exchange (ETDEWEB)

Tai, Xue-Cheng; Espedal, M. [Univ. of Bergen (Norway)

1996-12-31

Convergence of a space decomposition method is proved for a general convex programming problem. The space decomposition refers to methods that decompose a space into sums of subspaces, which could be a domain decomposition or a multigrid method for partial differential equations. Two algorithms are proposed. Both can be used for linear as well as nonlinear elliptic problems and they reduce to the standard additive and multiplicative Schwarz methods for linear elliptic problems. Two {open_quotes}hybrid{close_quotes} algorithms are also presented. They converge faster than the additive one and have better parallelism than the multiplicative method. Numerical tests with a two level domain decomposition for linear, nonlinear and interface elliptic problems are presented for the proposed algorithms.

Nature Inspired Computational Technique for the Numerical Solution of Nonlinear Singular Boundary Value Problems Arising in Physiology

Directory of Open Access Journals (Sweden)

Suheel Abdullah Malik

2014-01-01

Full Text Available We present a hybrid heuristic computing method for the numerical solution of nonlinear singular boundary value problems arising in physiology. The approximate solution is deduced as a linear combination of some log sigmoid basis functions. A fitness function representing the sum of the mean square error of the given nonlinear ordinary differential equation (ODE and its boundary conditions is formulated. The optimization of the unknown adjustable parameters contained in the fitness function is performed by the hybrid heuristic computation algorithm based on genetic algorithm (GA, interior point algorithm (IPA, and active set algorithm (ASA. The efficiency and the viability of the proposed method are confirmed by solving three examples from physiology. The obtained approximate solutions are found in excellent agreement with the exact solutions as well as some conventional numerical solutions.

Existence of Positive Solutions to a Boundary Value Problem for a Delayed Nonlinear Fractional Differential System

Directory of Open Access Journals (Sweden)

Chen Yuming

2011-01-01

Full Text Available Though boundary value problems for fractional differential equations have been extensively studied, most of the studies focus on scalar equations and the fractional order between 1 and 2. On the other hand, delay is natural in practical systems. However, not much has been done for fractional differential equations with delays. Therefore, in this paper, we consider a boundary value problem of a general delayed nonlinear fractional system. With the help of some fixed point theorems and the properties of the Green function, we establish several sets of sufficient conditions on the existence of positive solutions. The obtained results extend and include some existing ones and are illustrated with some examples for their feasibility.

Colored Contact Lens Dangers

Medline Plus

Full Text Available ... Cleveland. "This is far from the truth." Real People, Real Problems with Colored Contact Lenses Julian: Teenager ... about the members of the eye-care team . Consumer warning about the improper use of colored contact ...

Colored Contact Lens Dangers

Medline Plus

Full Text Available ... Cleveland. "This is far from the truth." Real People, Real Problems with Colored Contact Lenses Julian: Teenager ... the lenses. Never share contact lenses with another person. Get follow up exams with your eye care ...

A New Nonlinear Unit Root Test with Fourier Function

OpenAIRE

Güriş, Burak

2017-01-01

Traditional unit root tests display a tendency to be nonstationary in the case of structural breaks and nonlinearity. To eliminate this problem this paper proposes a new flexible Fourier form nonlinear unit root test. This test eliminates this problem to add structural breaks and nonlinearity together to the test procedure. In this test procedure, structural breaks are modeled by means of a Fourier function and nonlinear adjustment is modeled by means of an Exponential Smooth Threshold Autore...

A dynamic elastic-visco-plastic unilateral contact problem with normal damped response and Coulomb friction

Czech Academy of Sciences Publication Activity Database

Eck, Ch.; Jarušek, Jiří; Sofonea, M.

2010-01-01

Roč. 21, č. 3 (2010), s. 229-251 ISSN 0956-7925 R&D Projects: GA AV ČR IAA100750802 Institutional research plan: CEZ:AV0Z10190503 Keywords : elastic-vosco plastic material * dynamic contact problem * normal damped response * unilateral constraint * Coulomb friction * weak solution * penalitazion * smoothing Subject RIV: BA - General Mathematics Impact factor: 1.480, year: 2010 http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=7675484&fileId=S0956792510000045

Nonlinear surface Alfven waves

International Nuclear Information System (INIS)

Cramer, N.F.

1991-01-01

The problem of nonlinear surface Alfven waves propagating on an interface between a plasma and a vacuum is discussed, with dispersion provided by the finite-frequency effect, i.e. the finite ratio of the frequency to the ion-cyclotron frequency. A set of simplified nonlinear wave equations is derived using the method of stretched co-ordinates, and another approach uses the generation of a second-harmonic wave and its interaction with the first harmonic to obtain a nonlinear dispersion relation. A nonlinear Schroedinger equation is then derived, and soliton solutions found that propagate as solitary pulses in directions close to parallel and antiparallel to the background magnetic field. (author)

Nonlinear Physics of Plasmas

CERN Document Server

Kono, Mitsuo

2010-01-01

A nonlinearity is one of the most important notions in modern physics. A plasma is rich in nonlinearities and provides a variety of behaviors inherent to instabilities, coherent wave structures and turbulence. The book covers the basic concepts and mathematical methods, necessary to comprehend nonlinear problems widely encountered in contemporary plasmas, but also in other fields of physics and current research on self-organized structures and magnetized plasma turbulence. The analyses make use of strongly nonlinear models solved by analytical techniques backed by extensive simulations and available experiments. The text is written for senior undergraduates, graduate students, lecturers and researchers in laboratory, space and fusion plasmas.

An algorithm for nonlinear thermal analysis of fuel bearing pads

International Nuclear Information System (INIS)

Attia, M.H.; D'Silva, N.

1983-01-01

An algorithm has been developed for accurate prediction of the temperature field in a CANDU fuel bearing pad and the extent of the nucleate boiling in the crevice region. The methodology recognizes the nonlinear nature of the problem due to the fact that local boiling is both controlling and being controlled by the conditions of heat transfer at the boundaries. The finite difference model accounts for the volumetric effect of the thermal contact resistance at the bearing pad/pressure tube interface. It also allows the evaluation of the thermal barrier effect caused by applying an oxide film on the radiused surface of the bearing pad. Information pertaining to the distribution of the coefficient of heat transfer over water-cooled surfaces has been generated. Analysis of the results indicated the significance of considering the nonlinear behaviour of the system in determining its state of equilibrium. It also indicated that, depending on the thickness of the oxide layer and the position of the bearing pad along the core of the reactor, the nucleate boiling process can be prevented

Minimax terminal approach problem in two-level hierarchical nonlinear discrete-time dynamical system

Energy Technology Data Exchange (ETDEWEB)

Shorikov, A. F., E-mail: afshorikov@mail.ru [Ural Federal University, 19 S. Mira, Ekaterinburg, 620002, Russia Institute of Mathematics and Mechanics, Ural Branch of Russian Academy of Sciences, 16 S. Kovalevskaya, Ekaterinburg, 620990 (Russian Federation)

2015-11-30

We consider a discrete–time dynamical system consisting of three controllable objects. The motions of all objects are given by the corresponding vector nonlinear or linear discrete–time recurrent vector relations, and control system for its has two levels: basic (first or I level) that is dominating and subordinate level (second or II level) and both have different criterions of functioning and united a priori by determined informational and control connections defined in advance. For the dynamical system in question, we propose a mathematical formalization in the form of solving a multistep problem of two-level hierarchical minimax program control over the terminal approach process with incomplete information and give a general scheme for its solving.

Recent topics in nonlinear PDE

International Nuclear Information System (INIS)

Mimura, Masayasu; Nishida, Takaaki

1984-01-01

The meeting on the subject of nonlinear partial differential equations was held at Hiroshima University in February, 1983. Leading and active mathematicians were invited to talk on their current research interests in nonlinear pdes occuring in the areas of fluid dynamics, free boundary problems, population dynamics and mathematical physics. This volume contains the theory of nonlinear pdes and the related topics which have been recently developed in Japan. (Auth.)

Colored Contact Lens Dangers

Medline Plus

Full Text Available ... Hazard: The Hidden Dangers of Buying Decorative Contact Lenses Without a Prescription Leer en Español: Peligros asociados ... truth." Real People, Real Problems with Colored Contact Lenses Julian: Teenager Blinded In One Eye By Non- ...

Regularity of the solutions to a nonlinear boundary problem with indefinite weight

Directory of Open Access Journals (Sweden)

Aomar Anane

2011-01-01

Full Text Available In this paper we study the regularity of the solutions to the problemDelta_p u = |u|^{p−2}u in the bounded smooth domainOmega ⊂ R^N,with|∇u|^{p−2} partial_{nu} u = lambda V (x|u|^{p−2}u + h as a nonlinear boundary condition, where partial Omega is C^{2,beta}, with beta ∈]0, 1[, and V is a weight in L^s(partial Omega and h ∈ L^s(partial Omega for some s ≥ 1. We prove that all solutions are in L^{infty}(Omega cap L^{infty}(Omega, and using the D.Debenedetto’s theorem of regularity in [1] we conclude that those solutions are in C^{1,alpha} overline{Omega} for some alpha ∈ ]0, 1[.

On the Uniqueness of Solutions of a Nonlinear Elliptic Problem Arising in the Confinement of a Plasma in a Stellarator Device

International Nuclear Information System (INIS)

Diaz, J. I.; Galiano, G.; Padial, J. F.

1999-01-01

We study the uniqueness of solutions of a semilinear elliptic problem obtained from an inverse formulation when the nonlinear terms of the equation are prescribed in a general class of real functions. The inverse problem arises in the modeling of the magnetic confinement of a plasma in a Stellarator device. The uniqueness proof relies on an L ∞ -estimate on the solution of an auxiliary nonlocal problem formulated in terms of the relative rearrangement of a datum with respect to the solution

Dynamic Model of Contact Interface between Stator and Rotor

Directory of Open Access Journals (Sweden)

ZengHui Zhao

2013-01-01

Full Text Available Based on the equivalent principle, a linear spring contact model was established for the friction layer between stator and rotor. Different contact conditions were described by a distance index δ. Detailed analysis of the nonlinear contact behavior especially the static and dynamic slipping was carried on using a space-time equation. A contact deflection angle was proposed to quantitatively express the influence of friction force on the output performance. A more precision simulation model was established based on the theoretical analysis, and influences of different preload pressures and elastic modulus Em of friction layer on output performance were analyzed. The results showed the simulation results had very good consistency with experimental results, and the model could well reflect the output characteristics of contact interface.

Numerical simulation of a nonlinear coupled fluid-structure problem. Application to the design of naval nuclear propulsion structures; Modelisation et simulation numerique d'un probleme couple fluide/structure non lineaire: application au dimensionnement de structures nucleaires de propulsion navale

Energy Technology Data Exchange (ETDEWEB)

Sigrist, J.F

2004-11-15

The present work deals with the numerical simulation of a coupled fluid/structure problem with fluid free surface. A generic coupled fluid/structure system is defined, on which a linear problem (modal analysis) and a non-linear problem (temporal analysis) are stated. In the linear case, a strong coupled method is used. It is based on a finite element approach of the structure problem and a finite or a boundary element approach of the fluid problem. The coupled problem is formulated in terms of pressure and displacement, leading to a non-symmetric problem which is solved with an appropriate algorithm. In the non-linear case, the structure problem is described with non-linear equations of motion, whereas the fluid problem is modeled with the Stokes equations. The numerical resolution of the coupled problem is based on a weak coupling procedure. The fluid problem is solved with a finite volume technique, using a moving mesh technique to adjust the structure motion, a VOF method for the description of the free surface and the PISO algorithm for the time integration. The structure problem is solved with a finite element technique, using an explicit/implicit time integration algorithm. A procedure is developed in order to handle the coupling in space (fluid forces and structure displacement exchanges between fluid and structure mesh, fluid re-meshing) and in time (staggered explicit algorithm, dynamic filtering of numerical oscillations). The non linear coupled problem is solved using a CFD code, whose use for FSI problem is validated with a benchmark presented in this work. A comparison is proposed between numerical results and analytical solution for two elementary fluid problems. The validation process can be applied for any CFD numerical code. A numerical study is then proposed on the generic coupled case in order to describe the fluid/structure interaction phenomenon (added mass, displaced mass, mode coupling, influence of structural non-linearity). An industrial

SPARSE ELECTROMAGNETIC IMAGING USING NONLINEAR LANDWEBER ITERATIONS

KAUST Repository

Desmal, Abdulla

2015-07-29

A scheme for efficiently solving the nonlinear electromagnetic inverse scattering problem on sparse investigation domains is described. The proposed scheme reconstructs the (complex) dielectric permittivity of an investigation domain from fields measured away from the domain itself. Least-squares data misfit between the computed scattered fields, which are expressed as a nonlinear function of the permittivity, and the measured fields is constrained by the L0/L1-norm of the solution. The resulting minimization problem is solved using nonlinear Landweber iterations, where at each iteration a thresholding function is applied to enforce the sparseness-promoting L0/L1-norm constraint. The thresholded nonlinear Landweber iterations are applied to several two-dimensional problems, where the ``measured\\'\\' fields are synthetically generated or obtained from actual experiments. These numerical experiments demonstrate the accuracy, efficiency, and applicability of the proposed scheme in reconstructing sparse profiles with high permittivity values.

Shape optimization in 2D contact problems with given friction and a solution-dependent coefficient of friction

Czech Academy of Sciences Publication Activity Database

Haslinger, J.; Outrata, Jiří; Pathó, R.

2012-01-01

Roč. 20, č. 1 (2012), s. 31-59 ISSN 1877-0533 R&D Projects: GA AV ČR IAA100750802 Institutional research plan: CEZ:AV0Z10750506 Institutional support: RVO:67985556 Keywords : shape optimization * Signorini problem * model with given frinction * solution-dependent coefficient of friction * mathematical probrams with equilibrium constraints Subject RIV: BA - General Mathematics Impact factor: 1.036, year: 2012 http://library.utia.cas.cz/separaty/2012/MTR/outrata-shape optimization in 2d contact problems with given friction and a solution-dependent coefficient of friction .pdf

Cascaded nonlinearities for ultrafast nonlinear optical science and applications

DEFF Research Database (Denmark)

Bache, Morten

the cascading nonlinearity is investigated in detail, especially with focus on femtosecond energetic laser pulses being subjected to this nonlinear response. Analytical, numerical and experimental results are used to understand the cascading interaction and applications are demonstrated. The defocusing soliton...... observations with analogies in fiber optics are observed numerically and experimentally, including soliton self-compression, soliton-induced resonant radiation, supercontinuum generation, optical wavebreaking and shock-front formation. All this happens despite no waveguide being present, thanks...... is of particular interest here, since it is quite unique and provides the solution to a number of standing challenges in the ultrafast nonlinear optics community. It solves the problem of catastrophic focusing and formation of a filaments in bulk glasses, which even under controlled circumstances is limited...

The contact sport of rough surfaces

Science.gov (United States)

Carpick, Robert W.

2018-01-01

Describing the way two surfaces touch and make contact may seem simple, but it is not. Fully describing the elastic deformation of ideally smooth contacting bodies, under even low applied pressure, involves second-order partial differential equations and fourth-rank elastic constant tensors. For more realistic rough surfaces, the problem becomes a multiscale exercise in surface-height statistics, even before including complex phenomena such as adhesion, plasticity, and fracture. A recent research competition, the “Contact Mechanics Challenge” (1), was designed to test various approximate methods for solving this problem. A hypothetical rough surface was generated, and the community was invited to model contact with this surface with competing theories for the calculation of properties, including contact area and pressure. A supercomputer-generated numerical solution was kept secret until competition entries were received. The comparison of results (2) provides insights into the relative merits of competing models and even experimental approaches to the problem.

Patch testing with markers of fragrance contact allergy. Do clinical tests correspond to patients' self-reported problems?

Science.gov (United States)

Johansen, J D; Andersen, T F; Veien, N; Avnstorp, C; Andersen, K E; Menné, T

1997-03-01

The aim of the present study was to investigate the relationship between patients' own recognition of skin problems using consumer products and the results of patch testing with markers of fragrance sensitization. Eight hundred and eighty-four consecutive eczema patients, 18-69 years of age, filled in a questionnaire prior to patch testing with the European standard series. The questionnaire contained questions about skin symptoms from the use of scented and unscented products as well as skin reactions from contact with spices, flowers and citrus fruits that could indicate fragrance sensitivity. A highly significant association was found between reporting a history of visible skin symptoms from using scented products and a positive patch test to the fragrance mix, whereas no such relationship could be established to the Peru balsam in univariate or multivariate analysis. Our results suggest that the role of Peru balsam in detecting relevant fragrance contact allergy is limited, while most fragrance mix-positive patients are aware that the use of scented products may cause skin problems.

Shape Optimization in Contact Problems with Coulomb Friction and a Solution-Dependent Friction Coefficient

Czech Academy of Sciences Publication Activity Database

Beremlijski, P.; Outrata, Jiří; Haslinger, Jaroslav; Pathó, R.

2014-01-01

Roč. 52, č. 5 (2014), s. 3371-3400 ISSN 0363-0129 R&D Projects: GA ČR(CZ) GAP201/12/0671 Grant - others:GA MŠK(CZ) CZ.1.05/1.1.00/02.0070; GA MŠK(CZ) CZ.1.07/2.3.00/20.0070 Institutional support: RVO:67985556 ; RVO:68145535 Keywords : shape optimization * contact problems * Coulomb friction * solution-dependent coefficient of friction * mathematical programs with equilibrium constraints Subject RIV: BA - General Mathematics Impact factor: 1.463, year: 2014 http://library.utia.cas.cz/separaty/2014/MTR/outrata-0434234.pdf

Nonlinear filtering for LIDAR signal processing

Directory of Open Access Journals (Sweden)

D. G. Lainiotis

1996-01-01

Full Text Available LIDAR (Laser Integrated Radar is an engineering problem of great practical importance in environmental monitoring sciences. Signal processing for LIDAR applications involves highly nonlinear models and consequently nonlinear filtering. Optimal nonlinear filters, however, are practically unrealizable. In this paper, the Lainiotis's multi-model partitioning methodology and the related approximate but effective nonlinear filtering algorithms are reviewed and applied to LIDAR signal processing. Extensive simulation and performance evaluation of the multi-model partitioning approach and its application to LIDAR signal processing shows that the nonlinear partitioning methods are very effective and significantly superior to the nonlinear extended Kalman filter (EKF, which has been the standard nonlinear filter in past engineering applications.

Approximate Solution of Nonlinear Klein-Gordon Equation Using Sobolev Gradients

Directory of Open Access Journals (Sweden)

Nauman Raza

2016-01-01

Full Text Available The nonlinear Klein-Gordon equation (KGE models many nonlinear phenomena. In this paper, we propose a scheme for numerical approximation of solutions of the one-dimensional nonlinear KGE. A common approach to find a solution of a nonlinear system is to first linearize the equations by successive substitution or the Newton iteration method and then solve a linear least squares problem. Here, we show that it can be advantageous to form a sum of squared residuals of the nonlinear problem and then find a zero of the gradient. Our scheme is based on the Sobolev gradient method for solving a nonlinear least square problem directly. The numerical results are compared with Lattice Boltzmann Method (LBM. The L2, L∞, and Root-Mean-Square (RMS values indicate better accuracy of the proposed method with less computational effort.

Contact problem for a solid indenter and a viscoelastic half-space described by the spectrum of relaxation and retardation times

Science.gov (United States)

Stepanov, F. I.

2018-04-01

The mechanical properties of a material which is modeled by an exponential creep kernel characterized by a spectrum of relaxation and retardation times are studied. The research is carried out considering a contact problem for a solid indenter sliding over a viscoelastic half-space. The contact pressure, indentation depth of the indenter, and the deformation component of the friction coefficient are analyzed with respect to the case of half-space material modeled by single relaxation and retardation times.

An approximate JKR solution for a general contact, including rough contacts

Science.gov (United States)

Ciavarella, M.

2018-05-01

In the present note, we suggest a simple closed form approximate solution to the adhesive contact problem under the so-called JKR regime. The derivation is based on generalizing the original JKR energetic derivation assuming calculation of the strain energy in adhesiveless contact, and unloading at constant contact area. The underlying assumption is that the contact area distributions are the same as under adhesiveless conditions (for an appropriately increased normal load), so that in general the stress intensity factors will not be exactly equal at all contact edges. The solution is simply that the indentation is δ =δ1 -√{ 2 wA‧ /P″ } where w is surface energy, δ1 is the adhesiveless indentation, A‧ is the first derivative of contact area and P‧‧ the second derivative of the load with respect to δ1. The solution only requires macroscopic quantities, and not very elaborate local distributions, and is exact in many configurations like axisymmetric contacts, but also sinusoidal waves contact and correctly predicts some features of an ideal asperity model used as a test case and not as a real description of a rough contact problem. The solution permits therefore an estimate of the full solution for elastic rough solids with Gaussian multiple scales of roughness, which so far was lacking, using known adhesiveless simple results. The result turns out to depend only on rms amplitude and slopes of the surface, and as in the fractal limit, slopes would grow without limit, tends to the adhesiveless result - although in this limit the JKR model is inappropriate. The solution would also go to adhesiveless result for large rms amplitude of roughness hrms, irrespective of the small scale details, and in agreement with common sense, well known experiments and previous models by the author.

Convex models and probabilistic approach of nonlinear fatigue failure

International Nuclear Information System (INIS)

Qiu Zhiping; Lin Qiang; Wang Xiaojun

2008-01-01

This paper is concerned with the nonlinear fatigue failure problem with uncertainties in the structural systems. In the present study, in order to solve the nonlinear problem by convex models, the theory of ellipsoidal algebra with the help of the thought of interval analysis is applied. In terms of the inclusion monotonic property of ellipsoidal functions, the nonlinear fatigue failure problem with uncertainties can be solved. A numerical example of 25-bar truss structures is given to illustrate the efficiency of the presented method in comparison with the probabilistic approach

Point source identification in nonlinear advection–diffusion–reaction systems

International Nuclear Information System (INIS)

Mamonov, A V; Tsai, Y-H R

2013-01-01

We consider a problem of identification of point sources in time-dependent advection–diffusion systems with a nonlinear reaction term. The linear counterpart of the problem in question can be reduced to solving a system of nonlinear algebraic equations via the use of adjoint equations. We extend this approach by constructing an algorithm that solves the problem iteratively to account for the nonlinearity of the reaction term. We study the question of improving the quality of source identification by adding more measurements adaptively using the solution obtained previously with a smaller number of measurements. (paper)

Recent advances in multiparametric nonlinear programming

KAUST Repository

Domí nguez, Luis F.; Narciso, Diogo A.; Pistikopoulos, Efstratios N.

2010-01-01

In this paper, we present recent developments in multiparametric nonlinear programming. For the case of convex problems, we highlight key issues regarding the full characterization of the parametric solution space and we discuss, through an illustrative example problem, four alternative state-of-the-art multiparametric nonlinear programming algorithms. We also identify a number of main challenges for the non-convex case and highlight future research directions. © 2009 Elsevier Ltd. All rights reserved.

Recent advances in multiparametric nonlinear programming

KAUST Repository

Domínguez, Luis F.

2010-05-01

In this paper, we present recent developments in multiparametric nonlinear programming. For the case of convex problems, we highlight key issues regarding the full characterization of the parametric solution space and we discuss, through an illustrative example problem, four alternative state-of-the-art multiparametric nonlinear programming algorithms. We also identify a number of main challenges for the non-convex case and highlight future research directions. © 2009 Elsevier Ltd. All rights reserved.

Nonlinear diffusion problem arising in plasma physics

International Nuclear Information System (INIS)

Berryman, J.G.; Holland, C.J.

1978-01-01

In earlier studies of plasma diffusion with Okuda-Dawson scaling (D approx. n/sup -1/2/), perturbation theory indicated that arbitrary initial data should evolve rapidly toward the separation solution of the relevant nonlinear diffusion equation. Now a Lyapunov functional has been found which is strictly decreasing in time and bounded below. The rigorous proof that arbitrary initial data evolve toeard the separable solution is summarized. Rigorous bounds on the decay time are also presented

Use Residual Correction Method and Monotone Iterative Technique to Calculate the Upper and Lower Approximate Solutions of Singularly Perturbed Non-linear Boundary Value Problems

Directory of Open Access Journals (Sweden)

Chi-Chang Wang

2013-09-01

Full Text Available This paper seeks to use the proposed residual correction method in coordination with the monotone iterative technique to obtain upper and lower approximate solutions of singularly perturbed non-linear boundary value problems. First, the monotonicity of a non-linear differential equation is reinforced using the monotone iterative technique, then the cubic-spline method is applied to discretize and convert the differential equation into the mathematical programming problems of an inequation, and finally based on the residual correction concept, complex constraint solution problems are transformed into simpler questions of equational iteration. As verified by the four examples given in this paper, the method proposed hereof can be utilized to fast obtain the upper and lower solutions of questions of this kind, and to easily identify the error range between mean approximate solutions and exact solutions.

Multiple Positive Solutions of a Nonlinear Four-Point Singular Boundary Value Problem with a p-Laplacian Operator on Time Scales

Directory of Open Access Journals (Sweden)

Shihuang Hong

2009-01-01

Full Text Available We present sufficient conditions for the existence of at least twin or triple positive solutions of a nonlinear four-point singular boundary value problem with a p-Laplacian dynamic equation on a time scale. Our results are obtained via some new multiple fixed point theorems.

Exact multiplicity results for quasilinear boundary-value problems with cubic-like nonlinearities

Directory of Open Access Journals (Sweden)

Idris Addou

2000-01-01

Full Text Available We consider the boundary-value problem $$displaylines{ -(varphi_p (u'' =lambda f(u mbox{ in }(0,1 cr u(0 = u(1 =0,, }$$ where $p>1$, $lambda >0$ and $varphi_p (x =| x|^{p-2}x$. The nonlinearity $f$ is cubic-like with three distinct roots 0=a less than b less than c. By means of a quadrature method, we provide the exact number of solutions for all $lambda >0$. This way we extend a recent result, for $p=2$, by Korman et al. cite{KormanLiOuyang} to the general case $p>1$. We shall prove that when 1less than $pleq 2$ the structure of the solution set is exactly the same as that studied in the case $p=2$ by Korman et al. cite{KormanLiOuyang}, and strictly different in the case $p>2$.

Nonlinear optics principles and applications

CERN Document Server

Rottwitt, Karsten

2014-01-01

IntroductionReview of linear opticsInduced polarizationHarmonic oscillator modelLocal field correctionsEstimated nonlinear responseSummaryTime-domain material responseThe polarization time-response functionThe Born-Oppenheimer approximationRaman scattering response function of silicaSummaryMaterial response in the frequency domain, susceptibility tensorsThe susceptibility tensorThe induced polarization in the frequency domainSum of monochromatic fieldsThe prefactor to the induced polarizationThird-order polarization in the Born-Oppenheimer approximation in the frequency domainKramers-Kronig relationsSummarySymmetries in nonlinear opticsSpatial symmetriesSecond-order materialsThird-order nonlinear materialsCyclic coordinate-systemContracted notation for second-order susceptibility tensorsSummaryThe nonlinear wave equationMono and quasi-monochromatic beamsPlane waves - the transverse problemWaveguidesVectorial approachNonlinear birefringenceSummarySecond-order nonlinear effectsGeneral theoryCoupled wave theoryP...

Nonlinear differential equations

CERN Document Server

Struble, Raimond A

2017-01-01

Detailed treatment covers existence and uniqueness of a solution of the initial value problem, properties of solutions, properties of linear systems, stability of nonlinear systems, and two-dimensional systems. 1962 edition.

Network science, nonlinear science and infrastructure systems

CERN Document Server

2007-01-01

Network Science, Nonlinear Science and Infrastructure Systems has been written by leading scholars in these areas. Its express purpose is to develop common theoretical underpinnings to better solve modern infrastructural problems. It is felt by many who work in these fields that many modern communication problems, ranging from transportation networks to telecommunications, Internet, supply chains, etc., are fundamentally infrastructure problems. Moreover, these infrastructure problems would benefit greatly from a confluence of theoretical and methodological work done with the areas of Network Science, Dynamical Systems and Nonlinear Science. This book is dedicated to the formulation of infrastructural tools that will better solve these types of infrastructural problems. .

Finite elements of nonlinear continua

CERN Document Server

Oden, John Tinsley

1972-01-01

Geared toward undergraduate and graduate students, this text extends applications of the finite element method from linear problems in elastic structures to a broad class of practical, nonlinear problems in continuum mechanics. It treats both theory and applications from a general and unifying point of view.The text reviews the thermomechanical principles of continuous media and the properties of the finite element method, and then brings them together to produce discrete physical models of nonlinear continua. The mathematical properties of these models are analyzed, along with the numerical s

Advances in iterative methods for nonlinear equations

CERN Document Server

Busquier, Sonia

2016-01-01

This book focuses on the approximation of nonlinear equations using iterative methods. Nine contributions are presented on the construction and analysis of these methods, the coverage encompassing convergence, efficiency, robustness, dynamics, and applications. Many problems are stated in the form of nonlinear equations, using mathematical modeling. In particular, a wide range of problems in Applied Mathematics and in Engineering can be solved by finding the solutions to these equations. The book reveals the importance of studying convergence aspects in iterative methods and shows that selection of the most efficient and robust iterative method for a given problem is crucial to guaranteeing a good approximation. A number of sample criteria for selecting the optimal method are presented, including those regarding the order of convergence, the computational cost, and the stability, including the dynamics. This book will appeal to researchers whose field of interest is related to nonlinear problems and equations...

A non-linear branch and cut method for solving discrete minimum compliance problems to global optimality

DEFF Research Database (Denmark)

Stolpe, Mathias; Bendsøe, Martin P.

2007-01-01

This paper present some initial results pertaining to a search for globally optimal solutions to a challenging benchmark example proposed by Zhou and Rozvany. This means that we are dealing with global optimization of the classical single load minimum compliance topology design problem with a fixed...... finite element discretization and with discrete design variables. Global optimality is achieved by the implementation of some specially constructed convergent nonlinear branch and cut methods, based on the use of natural relaxations and by applying strengthening constraints (linear valid inequalities...

Controllability of nonlinear delay oscillating systems

Directory of Open Access Journals (Sweden)

Chengbin Liang

2017-05-01

Full Text Available In this paper, we study the controllability of a system governed by second order delay differential equations. We introduce a delay Gramian matrix involving the delayed matrix sine, which is used to establish sufficient and necessary conditions of controllability for the linear problem. In addition, we also construct a specific control function for controllability. For the nonlinear problem, we construct a control function and transfer the controllability problem to a fixed point problem for a suitable operator. We give a sufficient condition to guarantee the nonlinear delay system is controllable. Two examples are given to illustrate our theoretical results by calculating a specific control function and inverse of a delay Gramian matrix.

Thermal Aging Evaluation of Mod. 9Cr-1Mo Steel using Nonlinear Rayleigh Waves

Energy Technology Data Exchange (ETDEWEB)

Joo, Young-Sang; Kim, Hoe-Woong; Kim, Jong-Bum [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Marino, Daniel; Kim, Jin-Yeon; Jacobs, L.J [Georgia Institute of Technology, Atlanta (United States); Ruiz, Alberto [UMSNH, Morelia (Mexico)

2014-10-15

Thermal aging can pose a high risk to decreases in the mechanical properties such as strength or creep resistance. This can lead to an unexpected failure during long term operation. Nonlinear NDE techniques are preferred over conventional NDE techniques (linear ultrasonic measurements) because nonlinear ultrasonic techniques have shown their capability to detect a microstructural damage in the structures undergoing fatigue and creep. These nonlinear ultrasonic techniques make use of the fact that the dislocation density increases, which will create a nonlinear distortion of an ultrasonic wave; this damage causes the generation of measurable higher harmonic components in an initially mono-chromatic ultrasonic signal. This study investigates the recently developed non-contact nonlinear ultrasonic technique to detect the microstructural damage of mod. 9Cr-1Mo steel based on nonlinear Rayleigh wave with varying propagation distances. Nonlinear Rayleigh surface wave measurements using a non-contact, air-coupled ultrasonic transducer have been applied for the thermal aging evaluation of modified 9Cr-1Mo ferritic-martensitic steel. Thermal aging for various heat treatment times of mod.. 9Cr-1Mo steel specimens is performed to obtain the nucleation and growth of precipitated particles in specimens. The amplitudes of the first and second harmonics are measured along the propagation distance and the relative nonlinearity parameter is obtained from these amplitudes. The relative nonlinearity parameter shows a similar trend with the Rockwell C hardness.

Methodology for global nonlinear analysis of nuclear systems

International Nuclear Information System (INIS)

Cacuci, D.G.; Cacuci, G.L.

1987-01-01

This paper outlines a general method for globally computing the crucial features of nonlinear problems: bifurcations, limit points, saddle points, extrema (maxima and minima); our method also yields the local sensitivities (i.e., first order derivatives) of the system's state variables (e.g., fluxes, power, temperatures, flows) at any point in the system's phase space. We also present an application of this method to the nonlinear BWR model discussed in Refs. 8 and 11. The most significant novel feature of our method is the recasting of a general mathematical problem comprising three aspects: (1) nonlinear constrained optimization, (2) sensitivity analysis, into a fixed point problem of the form F[u(s), λ(s)] = 0 whose global zeros and singular points are related to the special features (i.e., extrema, bifurcations, etc.) of the original problem

Optimal contact definition for reconstruction of Contact Maps

Directory of Open Access Journals (Sweden)

Stehr Henning

2010-05-01

Full Text Available Abstract Background Contact maps have been extensively used as a simplified representation of protein structures. They capture most important features of a protein's fold, being preferred by a number of researchers for the description and study of protein structures. Inspired by the model's simplicity many groups have dedicated a considerable amount of effort towards contact prediction as a proxy for protein structure prediction. However a contact map's biological interest is subject to the availability of reliable methods for the 3-dimensional reconstruction of the structure. Results We use an implementation of the well-known distance geometry protocol to build realistic protein 3-dimensional models from contact maps, performing an extensive exploration of many of the parameters involved in the reconstruction process. We try to address the questions: a to what accuracy does a contact map represent its corresponding 3D structure, b what is the best contact map representation with regard to reconstructability and c what is the effect of partial or inaccurate contact information on the 3D structure recovery. Our results suggest that contact maps derived from the application of a distance cutoff of 9 to 11Å around the Cβ atoms constitute the most accurate representation of the 3D structure. The reconstruction process does not provide a single solution to the problem but rather an ensemble of conformations that are within 2Å RMSD of the crystal structure and with lower values for the pairwise average ensemble RMSD. Interestingly it is still possible to recover a structure with partial contact information, although wrong contacts can lead to dramatic loss in reconstruction fidelity. Conclusions Thus contact maps represent a valid approximation to the structures with an accuracy comparable to that of experimental methods. The optimal contact definitions constitute key guidelines for methods based on contact maps such as structure prediction through

Multidimensional nonlinear descriptive analysis

CERN Document Server

Nishisato, Shizuhiko

2006-01-01

Quantification of categorical, or non-numerical, data is a problem that scientists face across a wide range of disciplines. Exploring data analysis in various areas of research, such as the social sciences and biology, Multidimensional Nonlinear Descriptive Analysis presents methods for analyzing categorical data that are not necessarily sampled randomly from a normal population and often involve nonlinear relations. This reference not only provides an overview of multidimensional nonlinear descriptive analysis (MUNDA) of discrete data, it also offers new results in a variety of fields. The first part of the book covers conceptual and technical preliminaries needed to understand the data analysis in subsequent chapters. The next two parts contain applications of MUNDA to diverse data types, with each chapter devoted to one type of categorical data, a brief historical comment, and basic skills peculiar to the data types. The final part examines several problems and then concludes with suggestions for futu...

Sparse Reconstruction Schemes for Nonlinear Electromagnetic Imaging

KAUST Repository

Desmal, Abdulla

2016-03-01

Electromagnetic imaging is the problem of determining material properties from scattered fields measured away from the domain under investigation. Solving this inverse problem is a challenging task because (i) it is ill-posed due to the presence of (smoothing) integral operators used in the representation of scattered fields in terms of material properties, and scattered fields are obtained at a finite set of points through noisy measurements; and (ii) it is nonlinear simply due the fact that scattered fields are nonlinear functions of the material properties. The work described in this thesis tackles the ill-posedness of the electromagnetic imaging problem using sparsity-based regularization techniques, which assume that the scatterer(s) occupy only a small fraction of the investigation domain. More specifically, four novel imaging methods are formulated and implemented. (i) Sparsity-regularized Born iterative method iteratively linearizes the nonlinear inverse scattering problem and each linear problem is regularized using an improved iterative shrinkage algorithm enforcing the sparsity constraint. (ii) Sparsity-regularized nonlinear inexact Newton method calls for the solution of a linear system involving the Frechet derivative matrix of the forward scattering operator at every iteration step. For faster convergence, the solution of this matrix system is regularized under the sparsity constraint and preconditioned by leveling the matrix singular values. (iii) Sparsity-regularized nonlinear Tikhonov method directly solves the nonlinear minimization problem using Landweber iterations, where a thresholding function is applied at every iteration step to enforce the sparsity constraint. (iv) This last scheme is accelerated using a projected steepest descent method when it is applied to three-dimensional investigation domains. Projection replaces the thresholding operation and enforces the sparsity constraint. Numerical experiments, which are carried out using

Applications of nonlinear fiber optics

CERN Document Server

Agrawal, Govind

2008-01-01

* The only book describing applications of nonlinear fiber optics * Two new chapters on the latest developments: highly nonlinear fibers and quantum applications* Coverage of biomedical applications* Problems provided at the end of each chapterThe development of new highly nonlinear fibers - referred to as microstructured fibers, holey fibers and photonic crystal fibers - is the next generation technology for all-optical signal processing and biomedical applications. This new edition has been thoroughly updated to incorporate these key technology developments.The bo

Nonlinear optimization

CERN Document Server

Ruszczynski, Andrzej

2011-01-01

Optimization is one of the most important areas of modern applied mathematics, with applications in fields from engineering and economics to finance, statistics, management science, and medicine. While many books have addressed its various aspects, Nonlinear Optimization is the first comprehensive treatment that will allow graduate students and researchers to understand its modern ideas, principles, and methods within a reasonable time, but without sacrificing mathematical precision. Andrzej Ruszczynski, a leading expert in the optimization of nonlinear stochastic systems, integrates the theory and the methods of nonlinear optimization in a unified, clear, and mathematically rigorous fashion, with detailed and easy-to-follow proofs illustrated by numerous examples and figures. The book covers convex analysis, the theory of optimality conditions, duality theory, and numerical methods for solving unconstrained and constrained optimization problems. It addresses not only classical material but also modern top...

Advances in dynamic relaxation techniques for nonlinear finite element analysis

International Nuclear Information System (INIS)

Sauve, R.G.; Metzger, D.R.

1995-01-01

Traditionally, the finite element technique has been applied to static and steady-state problems using implicit methods. When nonlinearities exist, equilibrium iterations must be performed using Newton-Raphson or quasi-Newton techniques at each load level. In the presence of complex geometry, nonlinear material behavior, and large relative sliding of material interfaces, solutions using implicit methods often become intractable. A dynamic relaxation algorithm is developed for inclusion in finite element codes. The explicit nature of the method avoids large computer memory requirements and makes possible the solution of large-scale problems. The method described approaches the steady-state solution with no overshoot, a problem which has plagued researchers in the past. The method is included in a general nonlinear finite element code. A description of the method along with a number of new applications involving geometric and material nonlinearities are presented. They include: (1) nonlinear geometric cantilever plate; (2) moment-loaded nonlinear beam; and (3) creep of nuclear fuel channel assemblies

Automatic contact in DYNA3D for vehicle crashworthiness

International Nuclear Information System (INIS)

Whirley, R.G.; Engelmann, B.E.

1994-01-01

This paper presents a new formulation for the automatic definition and treatment of mechanical contact in explicit, nonlinear, finite element analysis. Automatic contact offers the benefits of significantly reduced model construction time and fewer opportunities for user error, but faces significant challenges in reliability and computational costs. The authors have used a new four-step automatic contact algorithm. Key aspects of the proposed method include (1) automatic identification of adjacent and opposite surfaces in the global search phase, and (2) the use of a smoothly varying surface normal that allows a consistent treatment of shell intersection and corner contact conditions without ad hoc rules. Three examples are given to illustrate the performance of the newly proposed algorithm in the public DYNA3D code

Spurious Solutions Of Nonlinear Differential Equations

Science.gov (United States)

Yee, H. C.; Sweby, P. K.; Griffiths, D. F.

1992-01-01

Report utilizes nonlinear-dynamics approach to investigate possible sources of errors and slow convergence and non-convergence of steady-state numerical solutions when using time-dependent approach for problems containing nonlinear source terms. Emphasizes implications for development of algorithms in CFD and computational sciences in general. Main fundamental conclusion of study is that qualitative features of nonlinear differential equations cannot be adequately represented by finite-difference method and vice versa.

Direct approach for solving nonlinear evolution and two-point

Indian Academy of Sciences (India)

Time-delayed nonlinear evolution equations and boundary value problems have a wide range of applications in science and engineering. In this paper, we implement the differential transform method to solve the nonlinear delay differential equation and boundary value problems. Also, we present some numerical examples ...

A new numerical scheme for non uniform homogenized problems: Application to the non linear Reynolds compressible equation

Directory of Open Access Journals (Sweden)

Buscaglia Gustavo C.

2001-01-01

Full Text Available A new numerical approach is proposed to alleviate the computational cost of solving non-linear non-uniform homogenized problems. The article details the application of the proposed approach to lubrication problems with roughness effects. The method is based on a two-parameter Taylor expansion of the implicit dependence of the homogenized coefficients on the average pressure and on the local value of the air gap thickness. A fourth-order Taylor expansion provides an approximation that is accurate enough to be used in the global problem solution instead of the exact dependence, without introducing significant errors. In this way, when solving the global problem, the solution of local problems is simply replaced by the evaluation of a polynomial. Moreover, the method leads naturally to Newton-Raphson nonlinear iterations, that further reduce the cost. The overall efficiency of the numerical methodology makes it feasible to apply rigorous homogenization techniques in the analysis of compressible fluid contact considering roughness effects. Previous work makes use of an heuristic averaging technique. Numerical comparison proves that homogenization-based methods are superior when the roughness is strongly anisotropic and not aligned with the flow direction.