Numerical solution of non-linear dual-phase-lag bioheat transfer equation within skin tissues.

Science.gov (United States)

Kumar, Dinesh; Kumar, P; Rai, K N

2017-11-01

This paper deals with numerical modeling and simulation of heat transfer in skin tissues using non-linear dual-phase-lag (DPL) bioheat transfer model under periodic heat flux boundary condition. The blood perfusion is assumed temperature-dependent which results in non-linear DPL bioheat transfer model in order to predict more accurate results. A numerical method of line which is based on finite difference and Runge-Kutta (4,5) schemes, is used to solve the present non-linear problem. Under specific case, the exact solution has been obtained and compared with the present numerical scheme, and we found that those are in good agreement. A comparison based on model selection criterion (AIC) has been made among non-linear DPL models when the variation of blood perfusion rate with temperature is of constant, linear and exponential type with the experimental data and it has been found that non-linear DPL model with exponential variation of blood perfusion rate is closest to the experimental data. In addition, it is found that due to absence of phase-lag phenomena in Pennes bioheat transfer model, it achieves steady state more quickly and always predict higher temperature than thermal and DPL non-linear models. The effect of coefficient of blood perfusion rate, dimensionless heating frequency and Kirchoff number on dimensionless temperature distribution has also been analyzed. The whole analysis is presented in dimensionless form. Copyright © 2017 Elsevier Inc. All rights reserved.

Nonlinear ordinary differential equations analytical approximation and numerical methods

CERN Document Server

Hermann, Martin

2016-01-01

The book discusses the solutions to nonlinear ordinary differential equations (ODEs) using analytical and numerical approximation methods. Recently, analytical approximation methods have been largely used in solving linear and nonlinear lower-order ODEs. It also discusses using these methods to solve some strong nonlinear ODEs. There are two chapters devoted to solving nonlinear ODEs using numerical methods, as in practice high-dimensional systems of nonlinear ODEs that cannot be solved by analytical approximate methods are common. Moreover, it studies analytical and numerical techniques for the treatment of parameter-depending ODEs. The book explains various methods for solving nonlinear-oscillator and structural-system problems, including the energy balance method, harmonic balance method, amplitude frequency formulation, variational iteration method, homotopy perturbation method, iteration perturbation method, homotopy analysis method, simple and multiple shooting method, and the nonlinear stabilized march...

Modeling nonlinearities in MEMS oscillators.

Science.gov (United States)

Agrawal, Deepak K; Woodhouse, Jim; Seshia, Ashwin A

2013-08-01

We present a mathematical model of a microelectromechanical system (MEMS) oscillator that integrates the nonlinearities of the MEMS resonator and the oscillator circuitry in a single numerical modeling environment. This is achieved by transforming the conventional nonlinear mechanical model into the electrical domain while simultaneously considering the prominent nonlinearities of the resonator. The proposed nonlinear electrical model is validated by comparing the simulated amplitude-frequency response with measurements on an open-loop electrically addressed flexural silicon MEMS resonator driven to large motional amplitudes. Next, the essential nonlinearities in the oscillator circuit are investigated and a mathematical model of a MEMS oscillator is proposed that integrates the nonlinearities of the resonator. The concept is illustrated for MEMS transimpedance-amplifier- based square-wave and sine-wave oscillators. Closed-form expressions of steady-state output power and output frequency are derived for both oscillator models and compared with experimental and simulation results, with a good match in the predicted trends in all three cases.

Evaluation of nonlinearity and validity of nonlinear modeling for complex time series.

Science.gov (United States)

Suzuki, Tomoya; Ikeguchi, Tohru; Suzuki, Masuo

2007-10-01

Even if an original time series exhibits nonlinearity, it is not always effective to approximate the time series by a nonlinear model because such nonlinear models have high complexity from the viewpoint of information criteria. Therefore, we propose two measures to evaluate both the nonlinearity of a time series and validity of nonlinear modeling applied to it by nonlinear predictability and information criteria. Through numerical simulations, we confirm that the proposed measures effectively detect the nonlinearity of an observed time series and evaluate the validity of the nonlinear model. The measures are also robust against observational noises. We also analyze some real time series: the difference of the number of chickenpox and measles patients, the number of sunspots, five Japanese vowels, and the chaotic laser. We can confirm that the nonlinear model is effective for the Japanese vowel /a/, the difference of the number of measles patients, and the chaotic laser.

Computation of Nonlinear Backscattering Using a High-Order Numerical Method

Science.gov (United States)

Fibich, G.; Ilan, B.; Tsynkov, S.

2001-01-01

The nonlinear Schrodinger equation (NLS) is the standard model for propagation of intense laser beams in Kerr media. The NLS is derived from the nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. In this study we use a fourth-order finite-difference method supplemented by special two-way artificial boundary conditions (ABCs) to solve the NLH as a boundary value problem. Our numerical methodology allows for a direct comparison of the NLH and NLS models and for an accurate quantitative assessment of the backscattered signal.

Numerical Investigation on the Directionality of Nonlinear Indicial Responses

International Nuclear Information System (INIS)

Yee, Kwan Jung; Hong, Sang Won; Lee, Dong Ho

2007-01-01

An unsteady Euler solver is modified to investigate the directionality of nonlinear indicial response to a step change in the angle of attack. An impulsive change in the angle of attack is incorporated by using the field velocity approach, which is known to decouple the step change in the angle of attack from a pitch rate. Numerical results are thoroughly compared against analytical results for two-dimensional indicial responses. The same method is applied to investigate the directionality of nonlinear indicial responses. It is found that directionality is mainly due to the asymmetry of initial shock locations. Since the directionality of the pitching moment responses is significant in the critical Mach number region, it is also shown that consideration of the directionality is crucial for accurate modeling of the nonlinear indicial functions

Numerical study of fractional nonlinear Schrodinger equations

KAUST Repository

Klein, Christian

2014-10-08

Using a Fourier spectral method, we provide a detailed numerical investigation of dispersive Schrödinger-type equations involving a fractional Laplacian in an one-dimensional case. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be identified. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states and the long-time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions of the fractional nonlinear Schrödinger equation.

Numeral eddy current sensor modelling based on genetic neural network

International Nuclear Information System (INIS)

Yu Along

2008-01-01

This paper presents a method used to the numeral eddy current sensor modelling based on the genetic neural network to settle its nonlinear problem. The principle and algorithms of genetic neural network are introduced. In this method, the nonlinear model parameters of the numeral eddy current sensor are optimized by genetic neural network (GNN) according to measurement data. So the method remains both the global searching ability of genetic algorithm and the good local searching ability of neural network. The nonlinear model has the advantages of strong robustness, on-line modelling and high precision. The maximum nonlinearity error can be reduced to 0.037% by using GNN. However, the maximum nonlinearity error is 0.075% using the least square method

Non-linear belt transient analysis. A hybrid model for numerical belt conveyor simulation

Energy Technology Data Exchange (ETDEWEB)

Harrison, A. [Scientific Solutions, Inc., Aurora, CO (United States)

2008-07-01

Frictional and rolling losses along a running conveyor are discussed due to their important influence on wave propagation during starting and stopping. Hybrid friction models allow belt rubber losses and material flexing to be included in the initial tension calculations prior to any dynamic analysis. Once running tensions are defined, a numerical integration method using non-linear stiffness gradients is used to generate transient forces during starting and stopping. A modified Euler integration technique is used to simulate the entire starting and stopping cycle in less than 0.1 seconds. The procedure enables a faster scrutiny of unforeseen conveyor design issues such as low belt tension zones and high forces at drives. (orig.)

Effect of P T symmetry on nonlinear waves for three-wave interaction models in the quadratic nonlinear media

Science.gov (United States)

Shen, Yujia; Wen, Zichao; Yan, Zhenya; Hang, Chao

2018-04-01

We study the three-wave interaction that couples an electromagnetic pump wave to two frequency down-converted daughter waves in a quadratic optical crystal and P T -symmetric potentials. P T symmetric potentials are shown to modulate stably nonlinear modes in two kinds of three-wave interaction models. The first one is a spatially extended three-wave interaction system with odd gain-and-loss distribution in the channel. Modulated by the P T -symmetric single-well or multi-well Scarf-II potentials, the system is numerically shown to possess stable soliton solutions. Via adiabatical change of system parameters, numerical simulations for the excitation and evolution of nonlinear modes are also performed. The second one is a combination of P T -symmetric models which are coupled via three-wave interactions. Families of nonlinear modes are found with some particular choices of parameters. Stable and unstable nonlinear modes are shown in distinct families by means of numerical simulations. These results will be useful to further investigate nonlinear modes in three-wave interaction models.

Feedback options in nonlinear numerical finance

DEFF Research Database (Denmark)

Hugger, Jens; Mashayekhi, Sima

2012-01-01

on an infinite slab is presented and boundary values on a bounded domain are derived. This bounded, nonlinear, 2 dimensional initial-boundary value problem is solved numerically using a number of standard finite difference schemes and the methods incorporated in the symbolic software Maple™....

Numerical method for the nonlinear Fokker-Planck equation

International Nuclear Information System (INIS)

Zhang, D.S.; Wei, G.W.; Kouri, D.J.; Hoffman, D.K.

1997-01-01

A practical method based on distributed approximating functionals (DAFs) is proposed for numerically solving a general class of nonlinear time-dependent Fokker-Planck equations. The method relies on a numerical scheme that couples the usual path-integral concept to the DAF idea. The high accuracy and reliability of the method are illustrated by applying it to an exactly solvable nonlinear Fokker-Planck equation, and the method is compared with the accurate K-point Stirling interpolation formula finite-difference method. The approach is also used successfully to solve a nonlinear self-consistent dynamic mean-field problem for which both the cumulant expansion and scaling theory have been found by Drozdov and Morillo [Phys. Rev. E 54, 931 (1996)] to be inadequate to describe the occurrence of a long-lived transient bimodality. The standard interpretation of the transient bimodality in terms of the flat region in the kinetic potential fails for the present case. An alternative analysis based on the effective potential of the Schroedinger-like Fokker-Planck equation is suggested. Our analysis of the transient bimodality is strongly supported by two examples that are numerically much more challenging than other examples that have been previously reported for this problem. copyright 1997 The American Physical Society

Analytical-numerical solution of a nonlinear integrodifferential equation in econometrics

Science.gov (United States)

Kakhktsyan, V. M.; Khachatryan, A. Kh.

2013-07-01

A mixed problem for a nonlinear integrodifferential equation arising in econometrics is considered. An analytical-numerical method is proposed for solving the problem. Some numerical results are presented.

Nonlinear behaviour of cantilevered carbon nanotube resonators based on a new nonlinear electrostatic load model

Science.gov (United States)

Farokhi, Hamed; Païdoussis, Michael P.; Misra, Arun K.

2018-04-01

The present study examines the nonlinear behaviour of a cantilevered carbon nanotube (CNT) resonator and its mass detection sensitivity, employing a new nonlinear electrostatic load model. More specifically, a 3D finite element model is developed in order to obtain the electrostatic load distribution on cantilevered CNT resonators. A new nonlinear electrostatic load model is then proposed accounting for the end effects due to finite length. Additionally, a new nonlinear size-dependent continuum model is developed for the cantilevered CNT resonator, employing the modified couple stress theory (to account for size-effects) together with the Kelvin-Voigt model (to account for nonlinear damping); the size-dependent model takes into account all sources of nonlinearity, i.e. geometrical and inertial nonlinearities as well as nonlinearities associated with damping, small-scale, and electrostatic load. The nonlinear equation of motion of the cantilevered CNT resonator is obtained based on the new models developed for the CNT resonator and the electrostatic load. The Galerkin method is then applied to the nonlinear equation of motion, resulting in a set of nonlinear ordinary differential equations, consisting of geometrical, inertial, electrical, damping, and size-dependent nonlinear terms. This high-dimensional nonlinear discretized model is solved numerically utilizing the pseudo-arclength continuation technique. The nonlinear static and dynamic responses of the system are examined for various cases, investigating the effect of DC and AC voltages, length-scale parameter, nonlinear damping, and electrostatic load. Moreover, the mass detection sensitivity of the system is examined for possible application of the CNT resonator as a nanosensor.

Aeroelastic oscillations of a cantilever with structural nonlinearities: theory and numerical simulation.

Energy Technology Data Exchange (ETDEWEB)

Robinson, Brandon [Carleton Univ., Ottawa, ON (Canada). Dept. of Civil and Environmental Engineering; Rocha da Costa, Leandro Jose [Carleton Univ., Ottawa, ON (Canada). Dept. of Civil and Environmental Engineering; Poirel, Dominique [Royal Military College of Canada, Kingston (Canada). Dept. of Mechanical and Aerospace Engineering; Pettit, Chris [US Naval Academy, Annapolis, MD (United States). Dept. of Mechanical and Aerospace Engineering; Khalil, Mohammad [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Sarkar, Abhijit [Carleton Univ., Ottawa, ON (Canada). Dept. of Civil and Environmental Engineering

2017-09-01

Our study details the derivation of the nonlinear equations of motion for the axial, biaxial bending and torsional vibrations of an aeroelastic cantilever undergoing rigid body (pitch) rotation at the base. The primary attenstion is focussed on the geometric nonlinearities of the system, whereby the aeroelastic load is modeled by the theory of linear quasisteady aerodynamics. This modelling effort is intended to mimic the wind-tunnel experimental setup at the Royal Military College of Canada. While the derivation closely follows the work of Hodges and Dowell [1] for rotor blades, this aeroelastic system contains new inertial terms which stem from the fundamentally different kinematics than those exhibited by helicopter or wind turbine blades. Using the Hamilton’s principle, a set of coupled nonlinear partial differential equations (PDEs) and an ordinary differential equation (ODE) are derived which describes the coupled axial-bending-bending-torsion-pitch motion of the aeroelastic cantilever with the pitch rotation. The finite dimensional approximation of the coupled system of PDEs are obtained using the Galerkin projection, leading to a coupled system of ODEs. Subsequently, these nonlinear ODEs are solved numerically using the built-in MATLAB implicit ODE solver and the associated numerical results are compared with those obtained using Houbolt’s method. It is demonstrated that the system undergoes coalescence flutter, leading to a limit cycle oscillation (LCO) due to coupling between the rigid body pitching mode and teh flexible mode arising from the flapwise bending motion.

The factors influencing nonlinear characteristics of the short-circuit current in dye-sensitized solar cells investigated by a numerical model.

Science.gov (United States)

Shi, Yushuai; Dong, Xiandui

2013-06-24

A numerical model for interpretation of the light-intensity-dependent nonlinear characteristics of the short-circuit current in dye-sensitized solar cells is suggested. The model is based on the continuity equation and includes the influences of the nongeminate recombination between electrons and electron acceptors in the electrolyte and the geminate recombination between electrons and oxidized dye molecules. The influences of the order and rate constant of the nongeminate recombination reaction, the light-absorption coefficient of the dye, the film thickness, the rate constant of geminate recombination, and the regeneration rate constant on the nonlinear characteristics of the short-circuit current are simulated and analyzed. It is proposed that superlinear and sublinear characteristics of the short-circuit current should be attributed to low electron-collection efficiency and low dye-regeneration efficiency, respectively. These results allow a deep understanding of the origin of the nonlinear characteristics of the short-circuit current in solar cells. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Analytical and numerical modeling of sandbanks dynamics

NARCIS (Netherlands)

Idier, Deborah; Astruc, Dominique

2003-01-01

Linear and nonlinear behavior of large-scale underwater bedform patterns like sandbanks are studied using linear stability analysis and numerical modeling. The model is based on depth-integrated hydrodynamics equations with a quadratic bottom friction law and a bed load sediment transport model

Modelization of highly nonlinear waves in coastal regions

Science.gov (United States)

Gouin, Maïté; Ducrozet, Guillaume; Ferrant, Pierre

2015-04-01

The proposed work deals with the development of a highly non-linear model for water wave propagation in coastal regions. The accurate modelization of surface gravity waves is of major interest in ocean engineering, especially in the field of marine renewable energy. These marine structures are intended to be settled in coastal regions where the effect of variable bathymetry may be significant on local wave conditions. This study presents a numerical model for the wave propagation with complex bathymetry. It is based on High-Order Spectral (HOS) method, initially limited to the propagation of non-linear wave fields over flat bottom. Such a model has been developed and validated at the LHEEA Lab. (Ecole Centrale Nantes) over the past few years and the current developments will enlarge its application range. This new numerical model will keep the interesting numerical properties of the original pseudo-spectral approach (convergence, efficiency with the use of FFTs, …) and enable the possibility to propagate highly non-linear wave fields over long time and large distance. Different validations will be provided in addition to the presentation of the method. At first, Bragg reflection will be studied with the proposed approach. If the Bragg condition is satisfied, the reflected wave generated by a sinusoidal bottom patch should be amplified as a result of resonant quadratic interactions between incident wave and bottom. Comparisons will be provided with experiments and reference solutions. Then, the method will be used to consider the transformation of a non-linear monochromatic wave as it propagates up and over a submerged bar. As the waves travel up the front slope of the bar, it steepens and high harmonics are generated due to non-linear interactions. Comparisons with experimental data will be provided. The different test cases will assess the accuracy and efficiency of the method proposed.

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; �...

Phenomenological modeling of nonlinear holograms based on metallic geometric metasurfaces.

Science.gov (United States)

Ye, Weimin; Li, Xin; Liu, Juan; Zhang, Shuang

2016-10-31

Benefiting from efficient local phase and amplitude control at the subwavelength scale, metasurfaces offer a new platform for computer generated holography with high spatial resolution. Three-dimensional and high efficient holograms have been realized by metasurfaces constituted by subwavelength meta-atoms with spatially varying geometries or orientations. Metasurfaces have been recently extended to the nonlinear optical regime to generate holographic images in harmonic generation waves. Thus far, there has been no vector field simulation of nonlinear metasurface holograms because of the tremendous computational challenge in numerically calculating the collective nonlinear responses of the large number of different subwavelength meta-atoms in a hologram. Here, we propose a general phenomenological method to model nonlinear metasurface holograms based on the assumption that every meta-atom could be described by a localized nonlinear polarizability tensor. Applied to geometric nonlinear metasurfaces, we numerically model the holographic images formed by the second-harmonic waves of different spins. We show that, in contrast to the metasurface holograms operating in the linear optical regime, the wavelength of incident fundamental light should be slightly detuned from the fundamental resonant wavelength to optimize the efficiency and quality of nonlinear holographic images. The proposed modeling provides a general method to simulate nonlinear optical devices based on metallic metasurfaces.

Numerical combination for nonlinear analysis of structures coupled to layered soils

Directory of Open Access Journals (Sweden)

Wagner Queiroz Silva

Full Text Available This paper presents an alternative coupling strategy between the Boundary Element Method (BEM and the Finite Element Method (FEM in order to create a computational code for the analysis of geometrical nonlinear 2D frames coupled to layered soils. The soil is modeled via BEM, considering multiple inclusions and internal load lines, through an alternative formulation to eliminate traction variables on subregions interfaces. A total Lagrangean formulation based on positions is adopted for the consideration of the geometric nonlinear behavior of frame structures with exact kinematics. The numerical coupling is performed by an algebraic strategy that extracts and condenses the equivalent soil's stiffness matrix and contact forces to be introduced into the frame structures hessian matrix and internal force vector, respectively. The formulation covers the analysis of shallow foundation structures and piles in any direction. Furthermore, the piles can pass through different layers. Numerical examples are shown in order to illustrate and confirm the accuracy and applicability of the proposed technique.

The numerical dynamic for highly nonlinear partial differential equations

Science.gov (United States)

Lafon, A.; Yee, H. C.

1992-01-01

Problems associated with the numerical computation of highly nonlinear equations in computational fluid dynamics are set forth and analyzed in terms of the potential ranges of spurious behaviors. A reaction-convection equation with a nonlinear source term is employed to evaluate the effects related to spatial and temporal discretizations. The discretization of the source term is described according to several methods, and the various techniques are shown to have a significant effect on the stability of the spurious solutions. Traditional linearized stability analyses cannot provide the level of confidence required for accurate fluid dynamics computations, and the incorporation of nonlinear analysis is proposed. Nonlinear analysis based on nonlinear dynamical systems complements the conventional linear approach and is valuable in the analysis of hypersonic aerodynamics and combustion phenomena.

Applicability of linear and non-linear potential flow models on a Wavestar float

DEFF Research Database (Denmark)

Bozonnet, Pauline; Dupin, Victor; Tona, Paolino

2017-01-01

as a model based on non-linear potential flow theory and weakscatterer hypothesis are successively considered. Simple tests, such as dip tests, decay tests and captive tests enable to highlight the improvements obtained with the introduction of nonlinearities. Float motion under wave actions and without...... control action, limited to small amplitude motion with a single float, is well predicted by the numerical models, including the linear one. Still, float velocity is better predicted by accounting for non-linear hydrostatic and Froude-Krylov forces.......Numerical models based on potential flow theory, including different types of nonlinearities are compared and validated against experimental data for the Wavestar wave energy converter technology. Exact resolution of the rotational motion, non-linear hydrostatic and Froude-Krylov forces as well...

From spiking neuron models to linear-nonlinear models.

Science.gov (United States)

Ostojic, Srdjan; Brunel, Nicolas

2011-01-20

Neurons transform time-varying inputs into action potentials emitted stochastically at a time dependent rate. The mapping from current input to output firing rate is often represented with the help of phenomenological models such as the linear-nonlinear (LN) cascade, in which the output firing rate is estimated by applying to the input successively a linear temporal filter and a static non-linear transformation. These simplified models leave out the biophysical details of action potential generation. It is not a priori clear to which extent the input-output mapping of biophysically more realistic, spiking neuron models can be reduced to a simple linear-nonlinear cascade. Here we investigate this question for the leaky integrate-and-fire (LIF), exponential integrate-and-fire (EIF) and conductance-based Wang-Buzsáki models in presence of background synaptic activity. We exploit available analytic results for these models to determine the corresponding linear filter and static non-linearity in a parameter-free form. We show that the obtained functions are identical to the linear filter and static non-linearity determined using standard reverse correlation analysis. We then quantitatively compare the output of the corresponding linear-nonlinear cascade with numerical simulations of spiking neurons, systematically varying the parameters of input signal and background noise. We find that the LN cascade provides accurate estimates of the firing rates of spiking neurons in most of parameter space. For the EIF and Wang-Buzsáki models, we show that the LN cascade can be reduced to a firing rate model, the timescale of which we determine analytically. Finally we introduce an adaptive timescale rate model in which the timescale of the linear filter depends on the instantaneous firing rate. This model leads to highly accurate estimates of instantaneous firing rates.

Numerical investigation of beam-driven PWFA in quasi-nonlinear regime

International Nuclear Information System (INIS)

Londrillo, P.; Gatti, C.; Ferrario, M.

2014-01-01

In beam-driven Plasma Based Wakefield Acceleration (PWFA), the quasi-nonlinear model has been designed to combine high efficient ‘blowout’ regimes, where cold and overdense driving electron beams form a totally rarefied plasma channel, with low charge beam distribution assuring the excited wakefield preserves relevant linear properties. This scheme can have applications in experimental facilities, like SPARC 150 MeV linac at LNF-INFN laboratories, where low-emittance, low-charge narrow electron beams can be produced to be injected on a preformed plasma channel. Here we present a preliminary numerical investigation of this configuration, using the fully 3D ALaDyn PIC code, as a preparatory work to design optimal conditions for the COMB experimental set-up. Specific numerical tools, having computational and diagnostic advantages in PWFA conditions and checks of the numerical outcomes with analytical results, are also presented and discussed

Numerical treatments for solving nonlinear mixed integral equation

Directory of Open Access Journals (Sweden)

M.A. Abdou

2016-12-01

Full Text Available We consider a mixed type of nonlinear integral equation (MNLIE of the second kind in the space C[0,T]×L2(Ω,T<1. The Volterra integral terms (VITs are considered in time with continuous kernels, while the Fredholm integral term (FIT is considered in position with singular general kernel. Using the quadratic method and separation of variables method, we obtain a nonlinear system of Fredholm integral equations (NLSFIEs with singular kernel. A Toeplitz matrix method, in each case, is then used to obtain a nonlinear algebraic system. Numerical results are calculated when the kernels take a logarithmic form or Carleman function. Moreover, the error estimates, in each case, are then computed.

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

A. C. Gilchrist; J. C. Suhling; T. J. Urbanik

1999-01-01

In this research, an investigation on the mechanical behavior of corrugated board has been performed using finite element analysis. Numerical finite element models for corrugated board geometries have been created and executed. Both geometric (large deformation) and material nonlinearities were included in the models. The analyses were performed using the commercial...

On the nonlinear modeling of ring oscillators

KAUST Repository

Elwakil, Ahmed S.

2009-06-01

We develop higher-order nonlinear models of three-stage and five-stage ring oscillators based on a novel inverter model. The oscillation condition and oscillation frequency are derived and compared to classical linear model analysis. Two important special cases for five-stage ring oscillators are also studied. Numerical simulations are shown. © 2009 World Scientific Publishing Company.

On the nonlinear modeling of ring oscillators

KAUST Repository

Elwakil, Ahmed S.; Salama, Khaled N.

2009-01-01

We develop higher-order nonlinear models of three-stage and five-stage ring oscillators based on a novel inverter model. The oscillation condition and oscillation frequency are derived and compared to classical linear model analysis. Two important special cases for five-stage ring oscillators are also studied. Numerical simulations are shown. © 2009 World Scientific Publishing Company.

A numerical study on the impact of nonlinear interactions on the amplitude of the migrating semidiurnal tide

Directory of Open Access Journals (Sweden)

C. M. Huang

2006-12-01

Full Text Available To quantitatively study the effects of nonlinear interactions on tide structure, a nonlinear numerical tidal model is developed, and the reliability and convergence of the adopted algorithm and coding are checked by numerical experiments. Under the same conditions as those employed by the GSWM-00 (Global Scale Wave Model 2000, our model provides the nonlinear quasi-steady solution of the migrating semidiurnal tide, which differs from the GSWM-00 result (the linear steady solution in the MLT region, especially above 100 km. Additionally, their amplitude difference displays a remarkable month-to-month variation, and its significant magnitudes occur during the month with strong semidiurnal tide. A quantitative analysis suggests that the main cause for the amplitude difference is that the initial migrating 12-h tide will interact with the mean flow as well as the nonlinearity-excited 6-h tide, and subsequently yield a new 12-h tidal part. Furthermore, our simulations also show that the mean flow/tidal interaction will significantly alter the background wind and temperature fields. The large magnitudes of the tidal amplitude difference and the background alteration indicate that the nonlinear processes involved in tidal propagations should be comprehensively considered in the description of global atmospheric dynamics in the MLT region. The comparisons among our simulations, the GSWMs and some observations of tides suggest that the nonlinearity-induced tidal structure variation could be a possible mechanism to account for some discrepancies between the GSWMs and the observations.

Time Domain Modeling and Simulation of Nonlinear Slender Viscoelastic Beams Associating Cosserat Theory and a Fractional Derivative Model

Directory of Open Access Journals (Sweden)

Adailton S. Borges

Full Text Available Abstract A broad class of engineering systems can be satisfactory modeled under the assumptions of small deformations and linear material properties. However, many mechanical systems used in modern applications, like structural elements typical of aerospace and petroleum industries, have been characterized by increased slenderness and high static and dynamic loads. In such situations, it becomes indispensable to consider the nonlinear geometric effects and/or material nonlinear behavior. At the same time, in many cases involving dynamic loads, there comes the need for attenuation of vibration levels. In this context, this paper describes the development and validation of numerical models of viscoelastic slender beam-like structures undergoing large displacements. The numerical approach is based on the combination of the nonlinear Cosserat beam theory and a viscoelastic model based on Fractional Derivatives. Such combination enables to derive nonlinear equations of motion that, upon finite element discretization, can be used for predicting the dynamic behavior of the structure in the time domain, accounting for geometric nonlinearity and viscoelastic damping. The modeling methodology is illustrated and validated by numerical simulations, the results of which are compared to others available in the literature.

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

Numerical studies on the electromagnetic properties of the nonlinear Lorentz Computational model for the dielectric media

International Nuclear Information System (INIS)

Abe, H.; Okuda, H.

1994-06-01

We study linear and nonlinear properties of a new computer simulation model developed to study the propagation of electromagnetic waves in a dielectric medium in the linear and nonlinear regimes. The model is constructed by combining a microscopic model used in the semi-classical approximation for the dielectric media and the particle model developed for the plasma simulations. It is shown that the model may be useful for studying linear and nonlinear wave propagation in the dielectric media

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

Science.gov (United States)

Shih, Tsan-Hsing; Liu, nan-Suey

2010-01-01

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

European Workshop on High Order Nonlinear Numerical Schemes for Evolutionary PDEs

CERN Document Server

Beaugendre, Héloïse; Congedo, Pietro; Dobrzynski, Cécile; Perrier, Vincent; Ricchiuto, Mario

2014-01-01

This book collects papers presented during the European Workshop on High Order Nonlinear Numerical Methods for Evolutionary PDEs (HONOM 2013) that was held at INRIA Bordeaux Sud-Ouest, Talence, France in March, 2013. The central topic is high order methods for compressible fluid dynamics. In the workshop, and in this proceedings, greater emphasis is placed on the numerical than the theoretical aspects of this scientific field. The range of topics is broad, extending through algorithm design, accuracy, large scale computing, complex geometries, discontinuous Galerkin, finite element methods, Lagrangian hydrodynamics, finite difference methods and applications and uncertainty quantification. These techniques find practical applications in such fields as fluid mechanics, magnetohydrodynamics, nonlinear solid mechanics, and others for which genuinely nonlinear methods are needed.

An efficient nonlinear finite-difference approach in the computational modeling of the dynamics of a nonlinear diffusion-reaction equation in microbial ecology.

Science.gov (United States)

Macías-Díaz, J E; Macías, Siegfried; Medina-Ramírez, I E

2013-12-01

In this manuscript, we present a computational model to approximate the solutions of a partial differential equation which describes the growth dynamics of microbial films. The numerical technique reported in this work is an explicit, nonlinear finite-difference methodology which is computationally implemented using Newton's method. Our scheme is compared numerically against an implicit, linear finite-difference discretization of the same partial differential equation, whose computer coding requires an implementation of the stabilized bi-conjugate gradient method. Our numerical results evince that the nonlinear approach results in a more efficient approximation to the solutions of the biofilm model considered, and demands less computer memory. Moreover, the positivity of initial profiles is preserved in the practice by the nonlinear scheme proposed. Copyright © 2013 Elsevier Ltd. All rights reserved.

Quasilinear Extreme Learning Machine Model Based Internal Model Control for Nonlinear Process

Directory of Open Access Journals (Sweden)

Dazi Li

2015-01-01

Full Text Available A new strategy for internal model control (IMC is proposed using a regression algorithm of quasilinear model with extreme learning machine (QL-ELM. Aimed at the chemical process with nonlinearity, the learning process of the internal model and inverse model is derived. The proposed QL-ELM is constructed as a linear ARX model with a complicated nonlinear coefficient. It shows some good approximation ability and fast convergence. The complicated coefficients are separated into two parts. The linear part is determined by recursive least square (RLS, while the nonlinear part is identified through extreme learning machine. The parameters of linear part and the output weights of ELM are estimated iteratively. The proposed internal model control is applied to CSTR process. The effectiveness and accuracy of the proposed method are extensively verified through numerical results.

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

Numerical analysis of nonlinear behavior of steel-concrete composite structures

Directory of Open Access Journals (Sweden)

Í.J.M. LEMES

Full Text Available Abstract This paper presents the development of an effective numerical formulation for the analysis of steel-concrete composite structures considering geometric and materials nonlinear effects. Thus, a methodology based on Refined Plastic Hinge Method (RPHM was developed and the stiffness parameters were obtained by homogenization of cross-section. The evaluation of structural elements strength is done through the Strain Compatibility Method (SCM. The Newton-Raphson Method with path-following strategies is adopted to solve nonlinear global and local (in cross-section level equations. The results are compared with experimental and numerical database presents in literature and a good accuracy is observed in composite cross sections, composite columns, and composite portal frames.

Numerical investigation into the highly nonlinear heat transfer equation with bremsstrahlung emission in the inertial confinement fusion plasmas

Energy Technology Data Exchange (ETDEWEB)

Habibi, M.; Oloumi, M.; Hosseinkhani, H.; Magidi, S. [Plasma and Fusion Research School, Nuclear Science and Technology Research Institute, Tehran (Iran, Islamic Republic of)

2015-10-15

A highly nonlinear parabolic partial differential equation that models the electron heat transfer process in laser inertial fusion has been solved numerically. The strong temperature dependence of the electron thermal conductivity and heat loss term (Bremsstrahlung emission) makes this a highly nonlinear process. In this case, an efficient numerical method is developed for the energy transport mechanism from the region of energy deposition into the ablation surface by a combination of the Crank-Nicolson scheme and the Newton-Raphson method. The quantitative behavior of the electron temperature and the comparison between analytic and numerical solutions are also investigated. For more clarification, the accuracy and conservation of energy in the computations are tested. The numerical results can be used to evaluate the nonlinear electron heat conduction, considering the released energy of the laser pulse at the Deuterium-Tritium (DT) targets and preheating by heat conduction ahead of a compression shock in the inertial confinement fusion (ICF) approach. (copyright 2015 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim)

Symmetries and discretizations of the O(3) nonlinear sigma model

Energy Technology Data Exchange (ETDEWEB)

Flore, Raphael [TPI, Universitaet Jena (Germany)

2011-07-01

Nonlinear sigma models possess many interesting properties like asymptotic freedom, confinement or dynamical mass generation, and hence serve as toy models for QCD and other theories. We derive a formulation of the N=2 supersymmetric extension of the O(3) nonlinear sigma model in terms of constrained field variables. Starting from this formulation, it is discussed how the model can be discretized in a way that maintains as many symmetries of the theory as possible. Finally, recent numerical results related to these discretizations are presented.

Numerical simulation of the nonlinear dynamics of packets of spiral density waves

International Nuclear Information System (INIS)

Korchagin, V.I.

1987-01-01

In a numerical experiment, the behavior of nonlinear packets of spiral density waves in a gas disk has been investigated for different initial wave amplitudes. If the amplitude of the density perturbations is small (<5%), the wave packet is drawn toward the center or toward the periphery of the disk in accordance with the linear theory. The behavior of linear packets of waves with wavelength comparable to the disk radius (R/sub d//lambda = 4) exhibits good agreement with the conclusions of the linear theory of tightly wound spiral waves. The dynamics of wave packets with initial density amplitudes 16, 30, 50% demonstrates the nonlinear nature of the behavior. THe behavior is governed by whether or not the nonlinear effects of higher than third order in the wave amplitude play a part. If the wave packet dynamics is determined by the cubic nonlinearity, the results of the numerical experiment are in qualitative and quantitative agreement with the nonlinear theory of short waves, although the characteristic scale of the packet and the wavelength are of the order of the disk radius. In the cases when the nonlinear effects of higher orders in the amplitude play an important part, the behavior of a packet does not differ qualitatively from the behavior predicted by the theory of cubic nonlinearity, but the nonlinear spreading of the packet takes place more rapidly

NON-LINEAR FINITE ELEMENT MODELING OF DEEP DRAWING PROCESS

Directory of Open Access Journals (Sweden)

Hasan YILDIZ

2004-03-01

Full Text Available Deep drawing process is one of the main procedures used in different branches of industry. Finding numerical solutions for determination of the mechanical behaviour of this process will save time and money. In die surfaces, which have complex geometries, it is hard to determine the effects of parameters of sheet metal forming. Some of these parameters are wrinkling, tearing, and determination of the flow of the thin sheet metal in the die and thickness change. However, the most difficult one is determination of material properties during plastic deformation. In this study, the effects of all these parameters are analyzed before producing the dies. The explicit non-linear finite element method is chosen to be used in the analysis. The numerical results obtained for non-linear material and contact models are also compared with the experiments. A good agreement between the numerical and the experimental results is obtained. The results obtained for the models are given in detail.

Polynomial model inversion control: numerical tests and applications

OpenAIRE

Novara, Carlo

2015-01-01

A novel control design approach for general nonlinear systems is described in this paper. The approach is based on the identification of a polynomial model of the system to control and on the on-line inversion of this model. Extensive simulations are carried out to test the numerical efficiency of the approach. Numerical examples of applicative interest are presented, concerned with control of the Duffing oscillator, control of a robot manipulator and insulin regulation in a type 1 diabetic p...

Numerical analysis of non-linear vibrations of a fractionally damped cylindrical shell under the conditions of combinational internal resonance

Directory of Open Access Journals (Sweden)

Rossikhin Yury A.

2018-01-01

Full Text Available Non-linear damped vibrations of a cylindrical shell embedded into a fractional derivative medium are investigated for the case of the combinational internal resonance, resulting in modal interaction, using two different numerical methods with further comparison of the results obtained. The damping properties of the surrounding medium are described by the fractional derivative Kelvin-Voigt model utilizing the Riemann-Liouville fractional derivatives. Within the first method, the generalized displacements of a coupled set of nonlinear ordinary differential equations of the second order are estimated using numerical solution of nonlinear multi-term fractional differential equations by the procedure based on the reduction of the problem to a system of fractional differential equations. According to the second method, the amplitudes and phases of nonlinear vibrations are estimated from the governing nonlinear differential equations describing amplitude-and-phase modulations for the case of the combinational internal resonance. A good agreement in results is declared.

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

A Nonlinear Dynamic Model and Free Vibration Analysis of Deployable Mesh Reflectors

Science.gov (United States)

Shi, H.; Yang, B.; Thomson, M.; Fang, H.

2011-01-01

This paper presents a dynamic model of deployable mesh reflectors, in which geometric and material nonlinearities of such a space structure are fully described. Then, by linearization around an equilibrium configuration of the reflector structure, a linearized model is obtained. With this linearized model, the natural frequencies and mode shapes of a reflector can be computed. The nonlinear dynamic model of deployable mesh reflectors is verified by using commercial finite element software in numerical simulation. As shall be seen, the proposed nonlinear model is useful for shape (surface) control of deployable mesh reflectors under thermal loads.

Research on nonlinear stochastic dynamical price model

International Nuclear Information System (INIS)

Li Jiaorui; Xu Wei; Xie Wenxian; Ren Zhengzheng

2008-01-01

In consideration of many uncertain factors existing in economic system, nonlinear stochastic dynamical price model which is subjected to Gaussian white noise excitation is proposed based on deterministic model. One-dimensional averaged Ito stochastic differential equation for the model is derived by using the stochastic averaging method, and applied to investigate the stability of the trivial solution and the first-passage failure of the stochastic price model. The stochastic price model and the methods presented in this paper are verified by numerical studies

Numerical simulation of nonlinear wave force on a quasi-ellipse caisson

Science.gov (United States)

Wang, Yongxue; Ren, Xiaozhong; Wang, Guoyu

2011-09-01

A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio B/D was shown to have significant influence on the non-dimensional transverse wave force.

Sinusoidal velaroidal shell – numerical modelling of the nonlinear ...

African Journals Online (AJOL)

The nonlinearity, applied to a sinusoidal velaroidal shell with the inner radius r0, the outer variables radii from 10m to 20m and the number of waves n=8, will give rise to the investigation of its nonlinear buckling resistance. The building material is a high-performant concrete. The investigation emphasizes more on the ...

Nonlinear flow model for well production in an underground formation

Directory of Open Access Journals (Sweden)

J. C. Guo

2013-05-01

Full Text Available Fluid flow in underground formations is a nonlinear process. In this article we modelled the nonlinear transient flow behaviour of well production in an underground formation. Based on Darcy's law and material balance equations, we used quadratic pressure gradients to deduce diffusion equations and discuss the origins of nonlinear flow issues. By introducing an effective-well-radius approach that considers skin factor, we established a nonlinear flow model for both gas and liquid (oil or water. The liquid flow model was solved using a semi-analytical method, while the gas flow model was solved using numerical simulations because the diffusion equation of gas flow is a stealth function of pressure. For liquid flow, a series of standard log-log type curves of pressure transients were plotted and nonlinear transient flow characteristics were analyzed. Qualitative and quantitative analyses were used to compare the solutions of the linear and nonlinear models. The effect of nonlinearity upon pressure transients should not be ignored. For gas flow, pressure transients were simulated and compared with oil flow under the same formation and well conditions, resulting in the conclusion that, under the same volume rate production, oil wells demand larger pressure drops than gas wells. Comparisons between theoretical data and field data show that nonlinear models will describe fluid flow in underground formations realistically and accurately.

Reduced-order computational model in nonlinear structural dynamics for structures having numerous local elastic modes in the low-frequency range. Application to fuel assemblies

International Nuclear Information System (INIS)

Batou, A.; Soize, C.; Brie, N.

2013-01-01

Highlights: • A ROM of a nonlinear dynamical structure is built with a global displacements basis. • The reduced order model of fuel assemblies is accurate and of very small size. • The shocks between grids of a row of seven fuel assemblies are computed. -- Abstract: We are interested in the construction of a reduced-order computational model for nonlinear complex dynamical structures which are characterized by the presence of numerous local elastic modes in the low-frequency band. This high modal density makes the use of the classical modal analysis method not suitable. Therefore the reduced-order computational model is constructed using a basis of a space of global displacements, which is constructed a priori and which allows the nonlinear dynamical response of the structure observed on the stiff part to be predicted with a good accuracy. The methodology is applied to a complex industrial structure which is made up of a row of seven fuel assemblies with possibility of collisions between grids and which is submitted to a seismic loading

Reduced-order computational model in nonlinear structural dynamics for structures having numerous local elastic modes in the low-frequency range. Application to fuel assemblies

Energy Technology Data Exchange (ETDEWEB)

Batou, A., E-mail: anas.batou@univ-paris-est.fr [Université Paris-Est, Laboratoire Modélisation et Simulation Multi Echelle, MSME UMR 8208 CNRS, 5 bd Descartes, 77454 Marne-la-Vallee (France); Soize, C., E-mail: christian.soize@univ-paris-est.fr [Université Paris-Est, Laboratoire Modélisation et Simulation Multi Echelle, MSME UMR 8208 CNRS, 5 bd Descartes, 77454 Marne-la-Vallee (France); Brie, N., E-mail: nicolas.brie@edf.fr [EDF R and D, Département AMA, 1 avenue du général De Gaulle, 92140 Clamart (France)

2013-09-15

Highlights: • A ROM of a nonlinear dynamical structure is built with a global displacements basis. • The reduced order model of fuel assemblies is accurate and of very small size. • The shocks between grids of a row of seven fuel assemblies are computed. -- Abstract: We are interested in the construction of a reduced-order computational model for nonlinear complex dynamical structures which are characterized by the presence of numerous local elastic modes in the low-frequency band. This high modal density makes the use of the classical modal analysis method not suitable. Therefore the reduced-order computational model is constructed using a basis of a space of global displacements, which is constructed a priori and which allows the nonlinear dynamical response of the structure observed on the stiff part to be predicted with a good accuracy. The methodology is applied to a complex industrial structure which is made up of a row of seven fuel assemblies with possibility of collisions between grids and which is submitted to a seismic loading.

Simplified Model of Nonlinear Landau Damping

International Nuclear Information System (INIS)

Yampolsky, N.A.; Fisch, N.J.

2009-01-01

The nonlinear interaction of a plasma wave with resonant electrons results in a plateau in the electron distribution function close to the phase velocity of the plasma wave. As a result, Landau damping of the plasma wave vanishes and the resonant frequency of the plasma wave downshifts. However, this simple picture is invalid when the external driving force changes the plasma wave fast enough so that the plateau cannot be fully developed. A new model to describe amplification of the plasma wave including the saturation of Landau damping and the nonlinear frequency shift is proposed. The proposed model takes into account the change of the plasma wave amplitude and describes saturation of the Landau damping rate in terms of a single fluid equation, which simplifies the description of the inherently kinetic nature of Landau damping. A proposed fluid model, incorporating these simplifications, is verified numerically using a kinetic Vlasov code.

Unstructured Spectral Element Model for Dispersive and Nonlinear Wave Propagation

DEFF Research Database (Denmark)

Engsig-Karup, Allan Peter; Eskilsson, Claes; Bigoni, Daniele

2016-01-01

We introduce a new stabilized high-order and unstructured numerical model for modeling fully nonlinear and dispersive water waves. The model is based on a nodal spectral element method of arbitrary order in space and a -transformed formulation due to Cai, Langtangen, Nielsen and Tveito (1998). In...

Controlled Nonlinear Stochastic Delay Equations: Part I: Modeling and Approximations

International Nuclear Information System (INIS)

Kushner, Harold J.

2012-01-01

This two-part paper deals with “foundational” issues that have not been previously considered in the modeling and numerical optimization of nonlinear stochastic delay systems. There are new classes of models, such as those with nonlinear functions of several controls (such as products), each with is own delay, controlled random Poisson measure driving terms, admissions control with delayed retrials, and others. There are two basic and interconnected themes for these models. The first, dealt with in this part, concerns the definition of admissible control. The classical definition of an admissible control as a nonanticipative relaxed control is inadequate for these models and needs to be extended. This is needed for the convergence proofs of numerical approximations for optimal controls as well as to have a well-defined model. It is shown that the new classes of admissible controls do not enlarge the range of the value functions, is closed (together with the associated paths) under weak convergence, and is approximatable by ordinary controls. The second theme, dealt with in Part II, concerns transportation equation representations, and their role in the development of numerical algorithms with much reduced memory and computational requirements.

On a method of numerical calculation of nonlinear radial pulsations of stars

International Nuclear Information System (INIS)

Kosovichev, A.G.

1984-01-01

Some features of using the finite difference method for numerical investigation of nonradial pulsations of stars were considered. The mathematical model of these pulsations is described by time-dependent gasdynaMic equations with gravity. A one-dimentional (spherically-symmetric) case is considered. It was obtained a two-parametric family of ultimate conservative difference schemes where the diffepence analogy of the main conservative laws as well as the additional relations for the balance to individual kinds of energy are performed. Such difference schemes provide more exact calculation of nonlinear flows with shocks as compared with the other difference schemes with the same order of approximation. The methods of numerical solution of implicit (absolute stable) difference schemes for a given family were considered. The coupled equations are solved through iterative Newton method Using martrix and separate successive eliminations. Numerical method can be used for calculation of large amplitude radial pulsations of stars

Numerical studies on soliton propagation in the dielectric media by the nonlinear Lorentz computational model

International Nuclear Information System (INIS)

Abe, H.; Okuda, H.

1994-06-01

Soliton propagation in the dielectric media has been simulated by using the nonlinear Lorentz computational model, which was recently developed to study the propagation of electromagnetic waves in a linear and a nonlinear dielectric. The model is constructed by combining a microscopic model used in the semi-classical approximation for dielectric media and the particle model developed for the plasma simulations. The carrier wave frequency is retained in the simulation so that not only the envelope of the soliton but also its phase can be followed in time. It is shown that the model may be useful for studying pulse propagation in the dielectric media

A stochastic delay model for pricing debt and equity: Numerical techniques and applications

Science.gov (United States)

Tambue, Antoine; Kemajou Brown, Elisabeth; Mohammed, Salah

2015-01-01

Delayed nonlinear models for pricing corporate liabilities and European options were recently developed. Using self-financed strategy and duplication we were able to derive a Random Partial Differential Equation (RPDE) whose solutions describe the evolution of debt and equity values of a corporate in the last delay period interval in the accompanied paper (Kemajou et al., 2012) [14]. In this paper, we provide robust numerical techniques to solve the delayed nonlinear model for the corporate value, along with the corresponding RPDEs modeling the debt and equity values of the corporate. Using financial data from some firms, we forecast and compare numerical solutions from both the nonlinear delayed model and classical Merton model with the real corporate data. From this comparison, it comes up that in corporate finance the past dependence of the firm value process may be an important feature and therefore should not be ignored.

Generalized Roe's numerical scheme for a two-fluid model

International Nuclear Information System (INIS)

Toumi, I.; Raymond, P.

1993-01-01

This paper is devoted to a mathematical and numerical study of a six equation two-fluid model. We will prove that the model is strictly hyperbolic due to the inclusion of the virtual mass force term in the phasic momentum equations. The two-fluid model is naturally written under a nonconservative form. To solve the nonlinear Riemann problem for this nonconservative hyperbolic system, a generalized Roe's approximate Riemann solver, is used, based on a linearization of the nonconservative terms. A Godunov type numerical scheme is built, using this approximate Riemann solver. 10 refs., 5 figs,

A Weakly Nonlinear Model for Kelvin–Helmholtz Instability in Incompressible Fluids

International Nuclear Information System (INIS)

Li-Feng, Wang; Wen-Hua, Ye; Zheng-Feng, Fan; Chuang, Xue; Ying-Jun, Li

2009-01-01

A weakly nonlinear model is proposed for the Kelvin–Helmholtz instability in two-dimensional incompressible fluids by expanding the perturbation velocity potential to third order. The third-order harmonic generation effects of single-mode perturbation are analyzed, as well as the nonlinear correction to the exponential growth of the fundamental modulation. The weakly nonlinear results are supported by numerical simulations. Density and resonance effects exist in the development of mode coupling. (fundamental areas of phenomenology (including applications))

Nonlinear model predictive control of a wave energy converter based on differential flatness parameterisation

Science.gov (United States)

Li, Guang

2017-01-01

This paper presents a fast constrained optimization approach, which is tailored for nonlinear model predictive control of wave energy converters (WEC). The advantage of this approach relies on its exploitation of the differential flatness of the WEC model. This can reduce the dimension of the resulting nonlinear programming problem (NLP) derived from the continuous constrained optimal control of WEC using pseudospectral method. The alleviation of computational burden using this approach helps to promote an economic implementation of nonlinear model predictive control strategy for WEC control problems. The method is applicable to nonlinear WEC models, nonconvex objective functions and nonlinear constraints, which are commonly encountered in WEC control problems. Numerical simulations demonstrate the efficacy of this approach.

Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms

KAUST Repository

Carrillo, José A.

2016-09-22

In this paper we present a numerical scheme for nonlinear continuity equations, which is based on the gradient flow formulation of an energy functional with respect to the quadratic transportation distance. It can be applied to a large class of nonlinear continuity equations, whose dynamics are driven by internal energies, given external potentials and/or interaction energies. The solver is based on its variational formulation as a gradient flow with respect to the Wasserstein distance. Positivity of solutions as well as energy decrease of the semi-discrete scheme are guaranteed by its construction. We illustrate this property with various examples in spatial dimension one and two.

Landau fluid model for weakly nonlinear dispersive magnetohydrodynamics

International Nuclear Information System (INIS)

Passot, T.; Sulem, P. L.

2005-01-01

In may astrophysical plasmas such as the solar wind, the terrestrial magnetosphere, or in the interstellar medium at small enough scales, collisions are negligible. When interested in the large-scale dynamics, a hydrodynamic approach is advantageous not only because its numerical simulations is easier than of the full Vlasov-Maxwell equations, but also because it provides a deep understanding of cross-scale nonlinear couplings. It is thus of great interest to construct fluid models that extended the classical magnetohydrodynamic (MHD) equations to collisionless situations. Two ingredients need to be included in such a model to capture the main kinetic effects: finite Larmor radius (FLR) corrections and Landau damping, the only fluid-particle resonance that can affect large scales and can be modeled in a relatively simple way. The Modelization of Landau damping in a fluid formalism is hardly possible in the framework of a systematic asymptotic expansion and was addressed mainly by means of parameter fitting in a linearized setting. We introduced a similar Landau fluid model but, that has the advantage of taking dispersive effects into account. This model properly describes dispersive MHD waves in quasi-parallel propagation. Since, by construction, the system correctly reproduces their linear dynamics, appropriate tests should address the nonlinear regime. In a first case, we show analytically that the weakly nonlinear modulational dynamics of quasi-parallel propagating Alfven waves is well captured. As a second test we consider the parametric decay instability of parallel Alfven waves and show that numerical simulations of the dispersive Landau fluid model lead to results that closely match the outcome of hybrid simulations. (Author)

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.

Force-controlled absorption in a fully-nonlinear numerical wave tank

International Nuclear Information System (INIS)

Spinneken, Johannes; Christou, Marios; Swan, Chris

2014-01-01

An active control methodology for the absorption of water waves in a numerical wave tank is introduced. This methodology is based upon a force-feedback technique which has previously been shown to be very effective in physical wave tanks. Unlike other methods, an a-priori knowledge of the wave conditions in the tank is not required; the absorption controller being designed to automatically respond to a wide range of wave conditions. In comparison to numerical sponge layers, effective wave absorption is achieved on the boundary, thereby minimising the spatial extent of the numerical wave tank. In contrast to the imposition of radiation conditions, the scheme is inherently capable of absorbing irregular waves. Most importantly, simultaneous generation and absorption can be achieved. This is an important advance when considering inclusion of reflective bodies within the numerical wave tank. In designing the absorption controller, an infinite impulse response filter is adopted, thereby eliminating the problem of non-causality in the controller optimisation. Two alternative controllers are considered, both implemented in a fully-nonlinear wave tank based on a multiple-flux boundary element scheme. To simplify the problem under consideration, the present analysis is limited to water waves propagating in a two-dimensional domain. The paper presents an extensive numerical validation which demonstrates the success of the method for a wide range of wave conditions including regular, focused and random waves. The numerical investigation also highlights some of the limitations of the method, particularly in simultaneously generating and absorbing large amplitude or highly-nonlinear waves. The findings of the present numerical study are directly applicable to related fields where optimum absorption is sought; these include physical wavemaking, wave power absorption and a wide range of numerical wave tank schemes

Numerical bifurcation analysis of a class of nonlinear renewal equations

NARCIS (Netherlands)

Breda, Dimitri; Diekmann, Odo; Liessi, Davide; Scarabel, Francesca

2016-01-01

We show, by way of an example, that numerical bifurcation tools for ODE yield reliable bifurcation diagrams when applied to the pseudospectral approximation of a one-parameter family of nonlinear renewal equations. The example resembles logistic-and Ricker-type population equations and exhibits

NUMERICAL WITHOUT ITERATION METHOD OF MODELING OF ELECTROMECHANICAL PROCESSES IN ASYNCHRONOUS ENGINES

Directory of Open Access Journals (Sweden)

D. G. Patalakh

2018-02-01

Full Text Available Purpose. Development of calculation of electromagnetic and electromechanic transients is in asynchronous engines without iterations. Methodology. Numeral methods of integration of usual differential equations, programming. Findings. As the system of equations, describing the dynamics of asynchronous engine, contents the products of rotor and stator currents and product of rotation frequency of rotor and currents, so this system is nonlinear one. The numeral solution of nonlinear differential equations supposes an iteration process on every step of integration. Time-continuing and badly converging iteration process may be the reason of calculation slowing. The improvement of numeral method by the way of an iteration process removing is offered. As result the modeling time is reduced. The improved numeral method is applied for integration of differential equations, describing the dynamics of asynchronous engine. Originality. The improvement of numeral method allowing to execute numeral integrations of differential equations containing product of functions is offered, that allows to avoid an iteration process on every step of integration and shorten modeling time. Practical value. On the basis of the offered methodology the universal program of modeling of electromechanics processes in asynchronous engines could be developed as taking advantage on fast-acting.

A non-linear model of economic production processes

Science.gov (United States)

Ponzi, A.; Yasutomi, A.; Kaneko, K.

2003-06-01

We present a new two phase model of economic production processes which is a non-linear dynamical version of von Neumann's neoclassical model of production, including a market price-setting phase as well as a production phase. The rate of an economic production process is observed, for the first time, to depend on the minimum of its input supplies. This creates highly non-linear supply and demand dynamics. By numerical simulation, production networks are shown to become unstable when the ratio of different products to total processes increases. This provides some insight into observed stability of competitive capitalist economies in comparison to monopolistic economies. Capitalist economies are also shown to have low unemployment.

Modeling of nonlinear responses for reciprocal transducers involving polarization switching

DEFF Research Database (Denmark)

Willatzen, Morten; Wang, Linxiang

2007-01-01

Nonlinearities and hysteresis effects in a reciprocal PZT transducer are examined by use of a dynamical mathematical model on the basis of phase-transition theory. In particular, we consider the perovskite piezoelectric ceramic in which the polarization process in the material can be modeled...... by Landau theory for the first-order phase transformation, in which each polarization state is associated with a minimum of the Landau free-energy function. Nonlinear constitutive laws are obtained by using thermodynamical equilibrium conditions, and hysteretic behavior of the material can be modeled...... intrinsically. The time-dependent Ginzburg-Landau theory is used in the parameter identification involving hysteresis effects. We use the Chebyshev collocation method in the numerical simulations. The elastic field is assumed to be coupled linearly with other fields, and the nonlinearity is in the E-D coupling...

Numerical solution of a non-linear conservation law applicable to the interior dynamics of partially molten planets

Science.gov (United States)

Bower, Dan J.; Sanan, Patrick; Wolf, Aaron S.

2018-01-01

The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. Crucially, in this formulation the effective or eddy diffusivity depends on the entropy gradient, ∂S / ∂r , as well as entropy itself. First we present a simplified model with semi-analytical solutions that highlights the large dynamic range of ∂S / ∂r -around 12 orders of magnitude-for physically-relevant parameters. It also elucidates the thermal structure of a magma ocean during the earliest stage of crystal formation. This motivates the development of a simple yet stable numerical scheme able to capture the large dynamic range of ∂S / ∂r and hence provide a flexible and robust method for time-integrating the energy equation. Using insight gained from the simplified model, we consider a full model, which includes energy fluxes associated with convection, mixing, gravitational separation, and conduction that all depend on the thermophysical properties of the melt and solid phases. This model is discretised and evolved by applying the finite volume method (FVM), allowing for extended precision calculations and using ∂S / ∂r as the solution variable. The FVM is well-suited to this problem since it is naturally energy conserving, flexible, and intuitive to incorporate arbitrary non-linear fluxes that rely on lookup data. Special attention is given to the numerically challenging scenario in which crystals first form in the centre of a magma ocean. The computational framework we devise is immediately applicable to modelling high melt fraction phenomena in Earth and planetary science research. Furthermore, it provides a template for solving similar non-linear diffusion equations that arise in other science and engineering disciplines, particularly for non-linear functional forms of the diffusion coefficient.

Anti-symmetrically fused model and non-linear integral equations in the three-state Uimin-Sutherland model

International Nuclear Information System (INIS)

Fujii, Akira; Kluemper, Andreas

1999-01-01

We derive the non-linear integral equations determining the free energy of the three-state pure bosonic Uimin-Sutherland model. In order to find a complete set of auxiliary functions, the anti-symmetric fusion procedure is utilized. We solve the non-linear integral equations numerically and see that the low-temperature behavior coincides with that predicted by conformal field theory. The magnetization and magnetic susceptibility are also calculated by means of the non-linear integral equation

International Benchmark on Numerical Simulations for 1D, Nonlinear Site Response (PRENOLIN) : Verification Phase Based on Canonical Cases

NARCIS (Netherlands)

Régnier, Julie; Bonilla, Luis-Fabian; Bard, Pierre-Yves; Bertrand, Etienne; Hollender, Fabrice; Kawase, Hiroshi; Sicilia, Deborah; Arduino, Pedro; Amorosi, Angelo; Asimaki, Dominiki; Pisano, F.

2016-01-01

PREdiction of NOn‐LINear soil behavior (PRENOLIN) is an international benchmark aiming to test multiple numerical simulation codes that are capable of predicting nonlinear seismic site response with various constitutive models. One of the objectives of this project is the assessment of the

Geometrically Nonlinear Transient Response of Laminated Plates with Nonlinear Elastic Restraints

Directory of Open Access Journals (Sweden)

Shaochong Yang

2017-01-01

Full Text Available To investigate the dynamic behavior of laminated plates with nonlinear elastic restraints, a varied constraint force model and a systematic numerical procedure are presented in this work. Several kinds of typical relationships of force-displacement for spring are established to simulate the nonlinear elastic restraints. In addition, considering the restraining moments of flexible pads, the pads are modeled by translational and rotational springs. The displacement- dependent constraint forces are added to the right-hand side of equations of motion and treated as additional applied loads. These loads can be explicitly defined, via an independent set of nonlinear load functions. The time histories of transverse displacements at typical points of the laminated plate are obtained through the transient analysis. Numerical examples show that the present method can effectively treat the geometrically nonlinear transient response of plates with nonlinear elastic restraints.

Mathematical and numerical study of nonlinear hyperbolic equations: model coupling and nonclassical shocks

International Nuclear Information System (INIS)

Boutin, B.

2009-11-01

This thesis concerns the mathematical and numerical study of nonlinear hyperbolic partial differential equations. A first part deals with an emergent problematic: the coupling of hyperbolic equations. The pursued applications are linked with the mathematical coupling of computing platforms, dedicated to an adaptative simulation of multi-scale phenomena. We propose and analyze a new coupling formalism based on extended PDE systems avoiding the geometric treatment of the interfaces. In addition, it allows to formulate the problem in a multidimensional setting, with possible covering of the coupled models. This formalism allows in particular to equip the coupling procedure with viscous regularization mechanisms, useful in the selection of natural discontinuous solutions. We analyze existence and uniqueness in the framework of a parabolic regularization a la Dafermos. Existence of a solution holds true under very general conditions but failure of uniqueness may naturally arise as soon as resonance occurs at the interfaces. Next, we highlight that our extended PDE framework gives rise to another regularization strategy based on thick interfaces. In this setting, we prove existence and uniqueness of the solutions of the Cauchy problem for initial data in L ∞ . The main tool consists in the derivation of a flexible and robust finite volume method for general triangulation which is analyzed in the setting of entropy measure-valued solutions by DiPerna. The second part is devoted to the definition of a finite volume scheme for the computing of nonclassical solutions of a scalar conservation law based on a kinetic relation. This scheme offers the feature to be stricto sensu conservative, in opposition to a Glimm approach that is only statistically conservative. The validity of our approach is illustrated through numerical examples. (author)

Numerical Modelling Of Pumpkin Balloon Instability

Science.gov (United States)

Wakefield, D.

Tensys have been involved in the numerical formfinding and load analysis of architectural stressed membrane structures for 15 years. They have recently broadened this range of activities into the `lighter than air' field with significant involvement in aerostat and heavy-lift hybrid airship design. Since early 2004 they have been investigating pumpkin balloon instability on behalf of the NASA ULDB programme. These studies are undertaken using inTENS, an in-house finite element program suite based upon the Dynamic Relaxation solution method and developed especially for the non-linear analysis and patterning of membrane structures. The paper describes the current state of an investigation that started with a numerical simulation of the lobed cylinder problem first studied by Calladine. The influence of material properties and local geometric deformation on stability is demonstrated. A number of models of complete pumpkin balloons have then been established, including a 64-gore balloon with geometry based upon Julian Nott's Endeavour. This latter clefted dramatically upon initial inflation, a phenomenon that has been reproduced in the numerical model. Ongoing investigations include the introduction of membrane contact modelling into inTENS and correlation studies with the series of large-scale ULDB models currently in preparation.

Design and numerical implementation of a 3-D non-linear viscoelastic constitutive model for brain tissue during impact

NARCIS (Netherlands)

Brands, D.W.A.; Peters, G.W.M.; Bovendeerd, P.H.M.

2004-01-01

Finite Element (FE) head models are often used to understand mechanical response of the head and its contents during impact loading in the head. CurrentFE models do not account for non-linear viscoelastic material behavior of brain tissue. We developed a new non-linear viscoelastic material model

An Efficient Numerical Approach for Solving Nonlinear Coupled Hyperbolic Partial Differential Equations with Nonlocal Conditions

Directory of Open Access Journals (Sweden)

A. H. Bhrawy

2014-01-01

Full Text Available One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initial-boundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initial-boundary coupled hyperbolic PDEs with variable coefficients as well as initial-nonlocal conditions. Using triangular, soliton, and exponential-triangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate.

Modeling and Analysis of a Nonlinear Age-Structured Model for Tumor Cell Populations with Quiescence

Science.gov (United States)

Liu, Zijian; Chen, Jing; Pang, Jianhua; Bi, Ping; Ruan, Shigui

2018-05-01

We present a nonlinear first-order hyperbolic partial differential equation model to describe age-structured tumor cell populations with proliferating and quiescent phases at the avascular stage in vitro. The division rate of the proliferating cells is assumed to be nonlinear due to the limitation of the nutrient and space. The model includes a proportion of newborn cells that enter directly the quiescent phase with age zero. This proportion can reflect the effect of treatment by drugs such as erlotinib. The existence and uniqueness of solutions are established. The local and global stabilities of the trivial steady state are investigated. The existence and local stability of the positive steady state are also analyzed. Numerical simulations are performed to verify the results and to examine the impacts of parameters on the nonlinear dynamics of the model.

An analog model for quantum lightcone fluctuations in nonlinear optics

International Nuclear Information System (INIS)

Ford, L.H.; De Lorenci, V.A.; Menezes, G.; Svaiter, N.F.

2013-01-01

We propose an analog model for quantum gravity effects using nonlinear dielectrics. Fluctuations of the spacetime lightcone are expected in quantum gravity, leading to variations in the flight times of pulses. This effect can also arise in a nonlinear material. We propose a model in which fluctuations of a background electric field, such as that produced by a squeezed photon state, can cause fluctuations in the effective lightcone for probe pulses. This leads to a variation in flight times analogous to that in quantum gravity. We make some numerical estimates which suggest that the effect might be large enough to be observable. - Highlights: ► Lightcone fluctuations, quantum fluctuations of the effective speed of light, are a feature of quantum gravity. ► Nonlinear dielectrics have a variable speed of light, analogous to the effects of gravity. ► Fluctuating electric fields create the effect of lightcone fluctuations in a nonlinear material. ► We propose to use squeezed light in a nonlinear material as an analog model of lightcone fluctuations. ► Variation in the speed of propagation of pulses is the observational signature of lightcone fluctuations.

Isotropic damage model and serial/parallel mix theory applied to nonlinear analysis of ferrocement thin walls. Experimental and numerical analysis

Directory of Open Access Journals (Sweden)

Jairo A. Paredes

2016-01-01

Full Text Available Ferrocement thin walls are the structural elements that comprise the earthquake resistant system of dwellings built with this material. This article presents the results drawn from an experimental campaign carried out over full-scale precast ferrocement thin walls that were assessed under lateral static loading conditions. The tests allowed the identification of structural parameters and the evaluation of the performance of the walls under static loading conditions. Additionally, an isotropic damage model for modelling the mortar was applied, as well as the classic elasto-plastic theory for modelling the meshes and reinforcing bars. The ferrocement is considered as a composite material, thus the serial/parallel mix theory is used for modelling its mechanical behavior. In this work a methodology for the numerical analysis that allows modeling the nonlinear behavior exhibited by ferrocement walls under static loading conditions, as well as their potential use in earthquake resistant design, is proposed.

Accuracy of three-dimensional seismic ground response analysis in time domain using nonlinear numerical simulations

Science.gov (United States)

Liang, Fayun; Chen, Haibing; Huang, Maosong

2017-07-01

To provide appropriate uses of nonlinear ground response analysis for engineering practice, a three-dimensional soil column with a distributed mass system and a time domain numerical analysis were implemented on the OpenSees simulation platform. The standard mesh of a three-dimensional soil column was suggested to be satisfied with the specified maximum frequency. The layered soil column was divided into multiple sub-soils with a different viscous damping matrix according to the shear velocities as the soil properties were significantly different. It was necessary to use a combination of other one-dimensional or three-dimensional nonlinear seismic ground analysis programs to confirm the applicability of nonlinear seismic ground motion response analysis procedures in soft soil or for strong earthquakes. The accuracy of the three-dimensional soil column finite element method was verified by dynamic centrifuge model testing under different peak accelerations of the earthquake. As a result, nonlinear seismic ground motion response analysis procedures were improved in this study. The accuracy and efficiency of the three-dimensional seismic ground response analysis can be adapted to the requirements of engineering practice.

NUMERICAL MODELLING OF THE SOIL BEHAVIOUR BY USING NEWLY DEVELOPED ADVANCED MATERIAL MODEL

Directory of Open Access Journals (Sweden)

Jan Veselý

2017-02-01

Full Text Available This paper describes a theoretical background, implementation and validation of the newly developed Jardine plastic hardening-softening model (JPHS model, which can be used for numerical modelling of the soils behaviour. Although the JPHS model is based on the elasto-plastic theory, like the Mohr-Coulomb model that is widely used in geotechnics, it contains some improvements, which removes the main disadvantages of the MC model. The presented model is coupled with an isotopically hardening and softening law, non-linear elastic stress-strain law, non-associated elasto-plastic material description and a cap yield surface. The validation of the model is done by comparing the numerical results with real measured data from the laboratory tests and by testing of the model on the real project of the tunnel excavation. The 3D numerical analysis is performed and the comparison between the JPHS, Mohr-Coulomb, Modified Cam-Clay, Hardening small strain model and monitoring in-situ data is done.

Reproducing the nonlinear dynamic behavior of a structured beam with a generalized continuum model

Science.gov (United States)

Vila, J.; Fernández-Sáez, J.; Zaera, R.

2018-04-01

In this paper we study the coupled axial-transverse nonlinear vibrations of a kind of one dimensional structured solids by application of the so called Inertia Gradient Nonlinear continuum model. To show the accuracy of this axiomatic model, previously proposed by the authors, its predictions are compared with numeric results from a previously defined finite discrete chain of lumped masses and springs, for several number of particles. A continualization of the discrete model equations based on Taylor series allowed us to set equivalent values of the mechanical properties in both discrete and axiomatic continuum models. Contrary to the classical continuum model, the inertia gradient nonlinear continuum model used herein is able to capture scale effects, which arise for modes in which the wavelength is comparable to the characteristic distance of the structured solid. The main conclusion of the work is that the proposed generalized continuum model captures the scale effects in both linear and nonlinear regimes, reproducing the behavior of the 1D nonlinear discrete model adequately.

Study of dispersive and nonlinear effects of coastal wave dynamics with a fully nonlinear potential flow model

Science.gov (United States)

Benoit, Michel; Yates, Marissa L.; Raoult, Cécile

2017-04-01

Efficient and accurate numerical models simulating wave propagation are required for a variety of engineering projects including the evaluation of coastal risks, the design of protective coastal structures, and the estimation of the potential for marine renewable energy devices. Nonlinear and dispersive effects are particularly significant in the coastal zone where waves interact with the bottom, the shoreline, and coastal structures. The main challenge in developing a numerical models is finding a compromise between computational efficiency and the required accuracy of the simulated wave field. Here, a potential approach is selected and the (fully nonlinear) water wave problem is formulated using the Euler-Zakharov equations (Zakharov, 1968) describing the temporal evolution of the free surface elevation and velocity potential. The proposed model (Yates and Benoit, 2015) uses a spectral approach in the vertical (i.e. the vertical variation of the potential is approximated by a linear combination of the first NT+1 Chebyshev polynomials, following the work of Tian and Sato (2008)). The Zakharov equations are integrated in time using a fourth-order Runge-Kutta scheme with a constant time step. At each sub-timestep, the Laplace Boundary Value Problem (BVP) is solved to estimate the free surface vertical velocity using the spectral approach, with typical values of NT between 5 to 8 for practical applications. The 1DH version of the code is validated with comparisons to the experimental data set of Becq-Girard et al. (1999), which studied the propagation of irregular waves over a beach profile with a submerged bar. The nonlinear and dispersive capacities of the model are verified with the correct representation of wave-wave interactions, in particular the transfer of energy between different harmonic components during wave propagation (analysis of the transformation of the variance spectrum along the channel). Evolution of wave skewness, asymmetry and kurtosis along the

Numerical Modeling and Mechanical Analysis of Flexible Risers

Directory of Open Access Journals (Sweden)

J. Y. Li

2015-01-01

Full Text Available ABAQUS is used to create a detailed finite element model for a 10-layer unbonded flexible riser to simulate the riser’s mechanical behavior under three load conditions: tension force and internal and external pressure. It presents a technique to create detailed finite element model and to analyze flexible risers. In FEM model, all layers are modeled separately with contact interfaces; interaction between steel trips in certain layers has been considered as well. FEM model considering contact interaction, geometric nonlinearity, and friction has been employed to accurately simulate the structural behavior of riser. The model includes the main features of the riser geometry with very little simplifying assumptions. The model was solved using a fully explicit time-integration scheme implemented in a parallel environment on an eight-processor cluster and 24 G memory computer. There is a very good agreement obtained from numerical and analytical comparisons, which validates the use of numerical model here. The results from the numerical simulation show that the numerical model takes into account various details of the riser. It has been shown that the detailed finite element model can be used to predict riser’s mechanics behavior under various load cases and bound conditions.

A Nonlinear Ship Manoeuvering Model: Identification and adaptive control with experiments for a model ship

Directory of Open Access Journals (Sweden)

Roger Skjetne

2004-01-01

Full Text Available Complete nonlinear dynamic manoeuvering models of ships, with numerical values, are hard to find in the literature. This paper presents a modeling, identification, and control design where the objective is to manoeuver a ship along desired paths at different velocities. Material from a variety of references have been used to describe the ship model, its difficulties, limitations, and possible simplifications for the purpose of automatic control design. The numerical values of the parameters in the model is identified in towing tests and adaptive manoeuvering experiments for a small ship in a marine control laboratory.

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.

Numerical solutions of a general coupled nonlinear system of parabolic and hyperbolic equations of thermoelasticity

Science.gov (United States)

Sweilam, N. H.; Abou Hasan, M. M.

2017-05-01

In this paper, the weighted-average non-standard finite-difference (WANSFD) method is used to study numerically the general time-fractional nonlinear, one-dimensional problem of thermoelasticity. This model contains the standard system arising in thermoelasticity as a special case. The stability of the proposed method is analyzed by a procedure akin to the standard John von Neumann technique. Moreover, the accuracy of the proposed scheme is proved. Numerical results are presented graphically, which reveal that the WANSFD method is easy to implement, effective and convenient for solving the proposed system. The proposed method could also be easily extended to solve other systems of fractional partial differential equations.

A deterministic combination of numerical and physical models for coastal waves

DEFF Research Database (Denmark)

Zhang, Haiwen

2006-01-01

of numerical and physical modelling hence provides an attractive alternative to the use of either tool on it's own. The goal of this project has been to develop a deterministically combined numerical/physical model where the physical wave tank is enclosed in a much larger computational domain, and the two......Numerical and physical modelling are the two main tools available for predicting the influence of water waves on coastlines and structures placed in the near-shore environment. Numerical models can cover large areas at the correct scale, but are limited in their ability to capture strong...... nonlinearities, wave breaking, splash, mixing, and other such complicated physics. Physical models naturally include the real physics (at the model scale), but are limited by the physical size of the facility and must contend with the fact that different physical effects scale differently. An integrated use...

Global-local nonlinear model reduction for flows in heterogeneous porous media

KAUST Repository

AlOtaibi, Manal; Calo, Victor M.; Efendiev, Yalchin R.; Galvis, Juan; Ghommem, Mehdi

2015-01-01

In this paper, we combine discrete empirical interpolation techniques, global mode decomposition methods, and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM), to reduce the computational complexity associated with nonlinear flows in highly-heterogeneous porous media. To solve the nonlinear governing equations, we employ the GMsFEM to represent the solution on a coarse grid with multiscale basis functions and apply proper orthogonal decomposition on a coarse grid. Computing the GMsFEM solution involves calculating the residual and the Jacobian on a fine grid. As such, we use local and global empirical interpolation concepts to circumvent performing these computations on the fine grid. The resulting reduced-order approach significantly reduces the flow problem size while accurately capturing the behavior of fully-resolved solutions. We consider several numerical examples of nonlinear multiscale partial differential equations that are numerically integrated using fully-implicit time marching schemes to demonstrate the capability of the proposed model reduction approach to speed up simulations of nonlinear flows in high-contrast porous media.

Global-local nonlinear model reduction for flows in heterogeneous porous media

KAUST Repository

AlOtaibi, Manal

2015-08-01

In this paper, we combine discrete empirical interpolation techniques, global mode decomposition methods, and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM), to reduce the computational complexity associated with nonlinear flows in highly-heterogeneous porous media. To solve the nonlinear governing equations, we employ the GMsFEM to represent the solution on a coarse grid with multiscale basis functions and apply proper orthogonal decomposition on a coarse grid. Computing the GMsFEM solution involves calculating the residual and the Jacobian on a fine grid. As such, we use local and global empirical interpolation concepts to circumvent performing these computations on the fine grid. The resulting reduced-order approach significantly reduces the flow problem size while accurately capturing the behavior of fully-resolved solutions. We consider several numerical examples of nonlinear multiscale partial differential equations that are numerically integrated using fully-implicit time marching schemes to demonstrate the capability of the proposed model reduction approach to speed up simulations of nonlinear flows in high-contrast porous media.

Nonlinear modeling and identification of a DC motor for bidirectional operation with real time experiments

International Nuclear Information System (INIS)

Kara, Tolgay; Eker, Ilyas

2004-01-01

Modeling and identification of mechanical systems constitute an essential stage in practical control design and applications. Controllers commanding systems that operate at varying conditions or require high precision operation raise the need for a nonlinear approach in modeling and identification. Most mechanical systems used in industry are composed of masses moving under the action of position and velocity dependent forces. These forces exhibit nonlinear behavior in certain regions of operation. For a multi-mass rotational system, the nonlinearities, like Coulomb friction and dead zone, significantly influence the system operation when the rotation changes direction. The paper presents nonlinear modeling and identification of a DC motor rotating in two directions together with real time experiments. Linear and nonlinear models for the system are obtained for identification purposes, and the major nonlinearities in the system, such as Coulomb friction and dead zone, are investigated and integrated in the nonlinear model. The Hammerstein nonlinear system approach is used for identification of the nonlinear system model. Online identification of the linear and nonlinear system models is performed using the recursive least squares method. Results of the real time experiments are graphically and numerically presented, and the advantages of the nonlinear identification approach are revealed

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.

Numerical Simulation of Freak Waves Based on the Four-Order Nonlinear Schr(o)dinger Equation

Institute of Scientific and Technical Information of China (English)

ZHANG Yun-qiu; ZHANG Ning-chuan; PEI Yu-guo

2007-01-01

A numerical wave model based on the modified four-order nonlinear Schrodinger (NLS) equation in deep water is developed to simulate freak waves. A standard split-step, pseudo-spectral method is used to solve NLS equation. The validation of the model is firstly verified, and then the simulation of freak waves is performed by changing sideband conditions. Results show that freak waves entirely consistent with the definition in the evolution of wave trains are obtained. The possible occurrence mechanism of freak waves is discussed and the relevant characteristics are also analyzed.

Experimental assessment and numerical modeling of the nonlinear behavior of the masonry shear walls under in-plane cyclic loading considering the brickwork-setting effect

Directory of Open Access Journals (Sweden)

Amir Hossein Karimi

2017-08-01

Full Text Available In this article, the main purpose is nonlinear analysis of the cyclic behavior of the masonry shear walls including brickwork setting using finite element method. Three different brickwork-settings including running bond style, herringbone style and Zarbi style (herreh style were investigated. To this end, the walls (in dimension of 195×1500×1720 mm were tested in the laboratory and then were simulated using macro modeling method by Abaqus software, and their hysteretic curves was drawn. The concrete damaged plasticity criteria in the Abaqus software is a model used in this research.In this method, the main failure mechanisms of fracture are cracking in tension and crushing in compression. The macro modeling method was used for numerical assessment of the masonry walls. After numerical modeling and drawing hysteretic curves and contrasting them with laboratory results, it was proven that the concrete damaged plasticity model, which is behavioral model for simulating concrete material, can be used for modeling masonry materials under seismic loading. However, this model cannot be used to simulate pinching effect in hysteretic curve drawn from seismic loading. The envelope curve resulted from the numerical analysis of all three brickwork layouts had a good agreement with the results of the laboratory tests, but in Hysteretic curve of Herringbone style and Zarbi style the pinching effect did not match experimental results

Numerical experiment for nonlinear full-wave tomography. 3; Hisenkei full wave tomography no suchi jikken

Energy Technology Data Exchange (ETDEWEB)

Tsuchiya, T [Dia Consultants Company, Tokyo (Japan)

1996-10-01

Nonlinear full-wave tomography (FWT) is under investigation to improve the estimation accuracy of Vp/Vs distributions. Full-wave tomography is one of the underground structure exploration methods mainly using Tarantola`s nonlinear local optimization method (LOM). Numerical experiment for FWT was carried out assuming relatively weak nonlinear underground structure. In the case of inversion by local optimization method, adequate preconditioning is important. Utilization of geological information is also effective in estimating low-frequency components of a model. As far as data are obtained under proper observation arrangement, even in actual field, precise estimation of Vp/Vs distributions is possible by FWT using explosion in a hole as wave source. In full-wave tomography, selection of observation arrangement is essential for both Vp and Vs. However, the proper arrangement is different between Vp and Vs. Approach to different analyses for Vp and Vs is also necessary by using only proper data for Vp and Vs among obtained data sets. 4 figs.

Parameter estimation in nonlinear models for pesticide degradation

International Nuclear Information System (INIS)

Richter, O.; Pestemer, W.; Bunte, D.; Diekkrueger, B.

1991-01-01

A wide class of environmental transfer models is formulated as ordinary or partial differential equations. With the availability of fast computers, the numerical solution of large systems became feasible. The main difficulty in performing a realistic and convincing simulation of the fate of a substance in the biosphere is not the implementation of numerical techniques but rather the incomplete data basis for parameter estimation. Parameter estimation is a synonym for statistical and numerical procedures to derive reasonable numerical values for model parameters from data. The classical method is the familiar linear regression technique which dates back to the 18th century. Because it is easy to handle, linear regression has long been established as a convenient tool for analysing relationships. However, the wide use of linear regression has led to an overemphasis of linear relationships. In nature, most relationships are nonlinear and linearization often gives a poor approximation of reality. Furthermore, pure regression models are not capable to map the dynamics of a process. Therefore, realistic models involve the evolution in time (and space). This leads in a natural way to the formulation of differential equations. To establish the link between data and dynamical models, numerical advanced parameter identification methods have been developed in recent years. This paper demonstrates the application of these techniques to estimation problems in the field of pesticide dynamics. (7 refs., 5 figs., 2 tabs.)

A new differential equations-based model for nonlinear history-dependent magnetic behaviour

International Nuclear Information System (INIS)

Aktaa, J.; Weth, A. von der

2000-01-01

The paper presents a new kind of numerical model describing nonlinear magnetic behaviour. The model is formulated as a set of differential equations taking into account history dependence phenomena like the magnetisation hysteresis as well as saturation effects. The capability of the model is demonstrated carrying out comparisons between measurements and calculations

Numerical Solution of Nonlinear Volterra Integral Equations System Using Simpson’s 3/8 Rule

Directory of Open Access Journals (Sweden)

Adem Kılıçman

2012-01-01

Full Text Available The Simpson’s 3/8 rule is used to solve the nonlinear Volterra integral equations system. Using this rule the system is converted to a nonlinear block system and then by solving this nonlinear system we find approximate solution of nonlinear Volterra integral equations system. One of the advantages of the proposed method is its simplicity in application. Further, we investigate the convergence of the proposed method and it is shown that its convergence is of order O(h4. Numerical examples are given to show abilities of the proposed method for solving linear as well as nonlinear systems. Our results show that the proposed method is simple and effective.

Geometrical nonlinear deformation model and its experimental study on bimorph giant magnetostrictive thin film

Institute of Scientific and Technical Information of China (English)

Wei LIU; Zhenyuan JIA; Fuji WANG; Yongshun ZHANG; Dongming GUO

2008-01-01

The geometrical nonlinearity of a giant magne-tostrictive thin film (GMF) can be clearly detected under the magnetostriction effect. Thus, using geometrical linear elastic theory to describe the strain, stress, and constitutive relationship of GMF is inaccurate. According to nonlinear elastic theory, a nonlinear deformation model of the bimorph GMF is established based on assumptions that the magnetostriction effect is equivalent to the effect of body force loaded on the GMF. With Taylor series method, the numerical solution is deduced. Experiments on TbDyFe/Polyimide (PI)/SmFe and TbDyFe/Cu/SmFe are then conducted to verify the proposed model, respectively. Results indicate that the nonlinear deflection curve model is in good conformity with the experimental data.

Development of a numerical 2-dimensional beach evolution model

DEFF Research Database (Denmark)

Baykal, Cüneyt

2014-01-01

This paper presents the description of a 2-dimensional numerical model constructed for the simulation of beach evolution under the action of wind waves only over the arbitrary land and sea topographies around existing coastal structures and formations. The developed beach evolution numerical model...... is composed of 4 submodels: a nearshore spectral wave transformation model based on an energy balance equation including random wave breaking and diffraction terms to compute the nearshore wave characteristics, a nearshore wave-induced circulation model based on the nonlinear shallow water equations...... to compute the nearshore depth-averaged wave-induced current velocities and mean water level changes, a sediment transport model to compute the local total sediment transport rates occurring under the action of wind waves, and a bottom evolution model to compute the bed level changes in time based...

Comparative Study of FDTD-Adopted Numerical Algorithms for Kerr Nonlinearities

DEFF Research Database (Denmark)

Maksymov, Ivan S.; Sukhorukov, Andrey A.; Lavrinenko, Andrei

2011-01-01

Accurate finite-difference time-domain (FDTD) modeling of optical pulse propagation in nonlinear media usually implies the use of auxiliary differential equation (ADE) techniques. The updating of electric field in full-vectorial 3-D ADE FDTD modeling of the optical Kerr effect and two-photon abso...... approaches. Such schemes can significantly reduce the CPU time for nonlinear computations, especially in 3-D models....

Adaptive numerical modeling of dynamic crack propagation

International Nuclear Information System (INIS)

Adouani, H.; Tie, B.; Berdin, C.; Aubry, D.

2006-01-01

We propose an adaptive numerical strategy that aims at developing reliable and efficient numerical tools to model dynamic crack propagation and crack arrest. We use the cohesive zone theory as behavior of interface-type elements to model crack. Since the crack path is generally unknown beforehand, adaptive meshing is proposed to model the dynamic crack propagation. The dynamic study requires the development of specific solvers for time integration. As both geometry and finite element mesh of the studied structure evolve in time during transient analysis, the stability behavior of dynamic solver becomes a major concern. For this purpose, we use the space-time discontinuous Galerkin finite element method, well-known to provide a natural framework to manage meshes that evolve in time. As an important result, we prove that the space-time discontinuous Galerkin solver is unconditionally stable, when the dynamic crack propagation is modeled by the cohesive zone theory, which is highly non-linear. (authors)

Numerical Oscillations Analysis for Nonlinear Delay Differential Equations in Physiological Control Systems

Directory of Open Access Journals (Sweden)

Qi Wang

2012-01-01

Full Text Available This paper deals with the oscillations of numerical solutions for the nonlinear delay differential equations in physiological control systems. The exponential θ-method is applied to p′(t=β0ωμp(t−τ/(ωμ+pμ(t−τ−γp(t and it is shown that the exponential θ-method has the same order of convergence as that of the classical θ-method. Several conditions under which the numerical solutions oscillate are derived. Moreover, it is proven that every nonoscillatory numerical solution tends to positive equilibrium of the continuous system. Finally, the main results are illustrated with numerical examples.

A new nonlinear turbulence model based on Partially-Averaged Navier-Stokes Equations

International Nuclear Information System (INIS)

Liu, J T; Wu, Y L; Cai, C; Liu, S H; Wang, L Q

2013-01-01

Partially-averaged Navier-Stokes (PANS) Model was recognized as a Reynolds-averaged Navier-Stokes (RANS) to direct numerical simulation (DNS) bridging method. PANS model was purported for any filter width-from RANS to DNS. PANS method also shared some similarities with the currently popular URANS (unsteady RANS) method. In this paper, a new PANS model was proposed, which was based on RNG k-ε turbulence model. The Standard and RNG k-ε turbulence model were both isotropic models, as well as PANS models. The sheer stress in those PANS models was solved by linear equation. The linear hypothesis was not accurate in the simulation of complex flow, such as stall phenomenon. The sheer stress here was solved by nonlinear method proposed by Ehrhard. Then, the nonlinear PANS model was set up. The pressure coefficient of the suction side of the NACA0015 hydrofoil was predicted. The result of pressure coefficient agrees well with experimental result, which proves that the nonlinear PANS model can capture the high pressure gradient flow. A low specific centrifugal pump was used to verify the capacity of the nonlinear PANS model. The comparison between the simulation results of the centrifugal pump and Particle Image Velocimetry (PIV) results proves that the nonlinear PANS model can be used in the prediction of complex flow field

Generalized Nonlinear Yule Models

OpenAIRE

Lansky, Petr; Polito, Federico; Sacerdote, Laura

2016-01-01

With the aim of considering models with persistent memory we propose a fractional nonlinear modification of the classical Yule model often studied in the context of macrovolution. Here the model is analyzed and interpreted in the framework of the development of networks such as the World Wide Web. Nonlinearity is introduced by replacing the linear birth process governing the growth of the in-links of each specific webpage with a fractional nonlinear birth process with completely general birth...

Comparisons of measurements and numerical simulations of a nonlinear beam subjected to a multi-harmonic non-ideal input signal

International Nuclear Information System (INIS)

Claeys, M.; Sinou, J.J.; Lambelin, J.P.; Alcoverro, B.

2014-01-01

This study presents a direct comparison of measured and predicted nonlinear vibrations of a clamped-clamped steel beam. A multi-harmonic comparison of simulations with measurements is performed at the vicinity of the primary resonance. First of all, a nonlinear analytical model of the beam is developed taking into account non-ideal boundary conditions. The Harmonic Balance Method is implemented to estimate the nonlinear behavior of the clamped-clamped beam. This nonlinear method enables to simulate the vibration stationary response of a nonlinear system projected on several harmonics. This study then proposes a method to compare numerical simulations with measurements on all these harmonics. A signal analysis tool is developed to extract the system harmonics' frequency responses from a temporal signal of a swept sine experiment. An evolutionary updating algorithm (Covariance Matrix Adaptation Evolution Strategy), coupled with highly selective filters is used to identify both fundamental frequency and harmonics' amplitude in the temporal signal, at every moment. This tool enables to extract the harmonic amplitudes of the output signal as well as the input signal. The input of the Harmonic Balance Method can then be either an ideal mono-harmonic signal or a multi-harmonic experimental input signal. Finally, the present work focuses on the comparison of experimental and simulated results. From experimental output harmonics and numerical simulations, it is shown that it is possible to distinguish the nonlinearities of the clamped-clamped beam and the effect of the non-ideal input signal. (authors)

Numerical model for wind-driven circulation in the Bay of Bengal

Digital Repository Service at National Institute of Oceanography (India)

Bahulayan, N.; Varadachari, V.V.R.

Wind-driven circulation in the Bay of Bengal, generated by a southwest wind of constant speed (10 m.sec -1) and direction (225 degrees TN), is presented. A non-linear hydrodynamic model is used for the simulation of circulation. Numerical...

NUMERICAL MODELING AND DYNAMIC SIMULATIONS OF NONLINEAR AEROTHERMOELASTIC OF A DOUBLE-WEDGE LIFTING SURFACE

Directory of Open Access Journals (Sweden)

ARIF A. EBRAHEEM AL-QASSAR

2008-12-01

Full Text Available The design of the re-entry space vehicles and high-speed aircrafts requires special attention to the nonlinear thermoelastic and aerodynamic instabilities of their structural components. The thermal effects are important since temperature environment influences significantly the static and dynamic behaviors of flight structures in supersonic/hypersonic regimes. To contribute to the understanding of dynamic behavior of these “hot” structures, a double-wedge lifting surface with combined freeplay and cubic stiffening structural nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes has been analyzed. A third order Piston Theory Aerodynamics is used to evaluate the applied nonlinear unsteady aerodynamic loads. The loss of torsional stiffness that may be incurred by lifting surfaces subjected to axial stresses induced by aerodynamic heating is also considered. The aerodynamic heating effect is estimated based on the adiabatic wall temperature due to high speed airstreams. Modelling issues as well as simulation results have been presented and pertinent conclusions outlined. It is highlighted that a serious loss of torsional stiffness may induce the dynamic instability of the lifting surfaces. The influence of various parameters such as flight condition, thickness ratio, freeplays and pitching stiffness nonlinearity are also discussed.

The influence of noise on nonlinear time series detection based on Volterra-Wiener-Korenberg model

Energy Technology Data Exchange (ETDEWEB)

Lei Min [State Key Laboratory of Vibration, Shock and Noise, Shanghai Jiao Tong University, Shanghai 200030 (China)], E-mail: leimin@sjtu.edu.cn; Meng Guang [State Key Laboratory of Vibration, Shock and Noise, Shanghai Jiao Tong University, Shanghai 200030 (China)

2008-04-15

This paper studies the influence of noises on Volterra-Wiener-Korenberg (VWK) nonlinear test model. Our numerical results reveal that different types of noises lead to different behavior of VWK model detection. For dynamic noise, it is difficult to distinguish chaos from nonchaotic but nonlinear determinism. For time series, measure noise has no impact on chaos determinism detection. This paper also discusses various behavior of VWK model detection with surrogate data for different noises.

Convergence Guaranteed Nonlinear Constraint Model Predictive Control via I/O Linearization

Directory of Open Access Journals (Sweden)

Xiaobing Kong

2013-01-01

Full Text Available Constituting reliable optimal solution is a key issue for the nonlinear constrained model predictive control. Input-output feedback linearization is a popular method in nonlinear control. By using an input-output feedback linearizing controller, the original linear input constraints will change to nonlinear constraints and sometimes the constraints are state dependent. This paper presents an iterative quadratic program (IQP routine on the continuous-time system. To guarantee its convergence, another iterative approach is incorporated. The proposed algorithm can reach a feasible solution over the entire prediction horizon. Simulation results on both a numerical example and the continuous stirred tank reactors (CSTR demonstrate the effectiveness of the proposed method.

A direct method for numerical solution of a class of nonlinear Volterra integro-differential equations and its application to the nonlinear fission and fusion reactor kinetics

International Nuclear Information System (INIS)

Nakahara, Yasuaki; Ise, Takeharu; Kobayashi, Kensuke; Itoh, Yasuyuki

1975-12-01

A new method has been developed for numerical solution of a class of nonlinear Volterra integro-differential equations with quadratic nonlinearity. After dividing the domain of the variable into subintervals, piecewise approximations are applied in the subintervals. The equation is first integrated over a subinterval to obtain the piecewise equation, to which six approximate treatments are applied, i.e. fully explicit, fully implicit, Crank-Nicolson, linear interpolation, quadratic and cubic spline. The numerical solution at each time step is obtained directly as a positive root of the resulting algebraic quadratic equation. The point reactor kinetics with a ramp reactivity insertion, linear temperature feedback and delayed neutrons can be described by one of this type of nonlinear Volterra integro-differential equations. The algorithm is applied to the Argonne benchmark problem and a model problem for a fast reactor without delayed neutrons. The fully implicit method has been found to be unconditionally stable in the sense that it always gives the positive real roots. The cubic spline method is divergent, and the other four methods are intermediate in between. From the estimation of the stability, convergency, accuracy and CPU time, it is concluded that the Crank-Nicolson method is best, then the linear interpolation method comes closely next to it. Discussions are also made on the possibility of applying the algorithm to the fusion reactor kinetics in the form of a nonlinear partial differential equation. (auth.)

A model reduction approach to numerical inversion for a parabolic partial differential equation

International Nuclear Information System (INIS)

Borcea, Liliana; Druskin, Vladimir; Zaslavsky, Mikhail; Mamonov, Alexander V

2014-01-01

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss–Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments. (paper)

A model reduction approach to numerical inversion for a parabolic partial differential equation

Science.gov (United States)

Borcea, Liliana; Druskin, Vladimir; Mamonov, Alexander V.; Zaslavsky, Mikhail

2014-12-01

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss-Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments.

PD/PID controller tuning based on model approximations: Model reduction of some unstable and higher order nonlinear models

Directory of Open Access Journals (Sweden)

Christer Dalen

2017-10-01

Full Text Available A model reduction technique based on optimization theory is presented, where a possible higher order system/model is approximated with an unstable DIPTD model by using only step response data. The DIPTD model is used to tune PD/PID controllers for the underlying possible higher order system. Numerous examples are used to illustrate the theory, i.e. both linear and nonlinear models. The Pareto Optimal controller is used as a reference controller.

Nonlinear System Identification Using Quasi-ARX RBFN Models with a Parameter-Classified Scheme

Directory of Open Access Journals (Sweden)

Lan Wang

2017-01-01

Full Text Available Quasi-linear autoregressive with exogenous inputs (Quasi-ARX models have received considerable attention for their usefulness in nonlinear system identification and control. In this paper, identification methods of quasi-ARX type models are reviewed and categorized in three main groups, and a two-step learning approach is proposed as an extension of the parameter-classified methods to identify the quasi-ARX radial basis function network (RBFN model. Firstly, a clustering method is utilized to provide statistical properties of the dataset for determining the parameters nonlinear to the model, which are interpreted meaningfully in the sense of interpolation parameters of a local linear model. Secondly, support vector regression is used to estimate the parameters linear to the model; meanwhile, an explicit kernel mapping is given in terms of the nonlinear parameter identification procedure, in which the model is transformed from the nonlinear-in-nature to the linear-in-parameter. Numerical and real cases are carried out finally to demonstrate the effectiveness and generalization ability of the proposed method.

Numerical Treatment of Fixed Point Applied to the Nonlinear Fredholm Integral Equation

Directory of Open Access Journals (Sweden)

Berenguer MI

2009-01-01

Full Text Available The authors present a method of numerical approximation of the fixed point of an operator, specifically the integral one associated with a nonlinear Fredholm integral equation, that uses strongly the properties of a classical Schauder basis in the Banach space .

Linear theory for filtering nonlinear multiscale systems with model error.

Science.gov (United States)

Berry, Tyrus; Harlim, John

2014-07-08

In this paper, we study filtering of multiscale dynamical systems with model error arising from limitations in resolving the smaller scale processes. In particular, the analysis assumes the availability of continuous-time noisy observations of all components of the slow variables. Mathematically, this paper presents new results on higher order asymptotic expansion of the first two moments of a conditional measure. In particular, we are interested in the application of filtering multiscale problems in which the conditional distribution is defined over the slow variables, given noisy observation of the slow variables alone. From the mathematical analysis, we learn that for a continuous time linear model with Gaussian noise, there exists a unique choice of parameters in a linear reduced model for the slow variables which gives the optimal filtering when only the slow variables are observed. Moreover, these parameters simultaneously give the optimal equilibrium statistical estimates of the underlying system, and as a consequence they can be estimated offline from the equilibrium statistics of the true signal. By examining a nonlinear test model, we show that the linear theory extends in this non-Gaussian, nonlinear configuration as long as we know the optimal stochastic parametrization and the correct observation model. However, when the stochastic parametrization model is inappropriate, parameters chosen for good filter performance may give poor equilibrium statistical estimates and vice versa; this finding is based on analytical and numerical results on our nonlinear test model and the two-layer Lorenz-96 model. Finally, even when the correct stochastic ansatz is given, it is imperative to estimate the parameters simultaneously and to account for the nonlinear feedback of the stochastic parameters into the reduced filter estimates. In numerical experiments on the two-layer Lorenz-96 model, we find that the parameters estimated online , as part of a filtering

Numerical Material Model for Composite Laminates in High-Velocity Impact Simulation

Directory of Open Access Journals (Sweden)

Tao Liu

Full Text Available Abstract A numerical material model for composite laminate, was developed and integrated into the nonlinear dynamic explicit finite element programs as a material user subroutine. This model coupling nonlinear state of equation (EOS, was a macro-mechanics model, which was used to simulate the major mechanical behaviors of composite laminate under high-velocity impact conditions. The basic theoretical framework of the developed material model was introduced. An inverse flyer plate simulation was conducted, which demonstrated the advantage of the developed model in characterizing the nonlinear shock response. The developed model and its implementation were validated through a classic ballistic impact issue, i.e. projectile impacting on Kevlar29/Phenolic laminate. The failure modes and ballistic limit velocity were analyzed, and a good agreement was achieved when comparing with the analytical and experimental results. The computational capacity of this model, for Kevlar/Epoxy laminates with different architectures, i.e. plain-woven and cross-plied laminates, was further evaluated and the residual velocity curves and damage cone were accurately predicted.

A Numerical Model for Trickle Bed Reactors

Science.gov (United States)

Propp, Richard M.; Colella, Phillip; Crutchfield, William Y.; Day, Marcus S.

2000-12-01

Trickle bed reactors are governed by equations of flow in porous media such as Darcy's law and the conservation of mass. Our numerical method for solving these equations is based on a total-velocity splitting, sequential formulation which leads to an implicit pressure equation and a semi-implicit mass conservation equation. We use high-resolution finite-difference methods to discretize these equations. Our solution scheme extends previous work in modeling porous media flows in two ways. First, we incorporate physical effects due to capillary pressure, a nonlinear inlet boundary condition, spatial porosity variations, and inertial effects on phase mobilities. In particular, capillary forces introduce a parabolic component into the recast evolution equation, and the inertial effects give rise to hyperbolic nonconvexity. Second, we introduce a modification of the slope-limiting algorithm to prevent our numerical method from producing spurious shocks. We present a numerical algorithm for accommodating these difficulties, show the algorithm is second-order accurate, and demonstrate its performance on a number of simplified problems relevant to trickle bed reactor modeling.

Numerical simulation and comparison of nonlinear self-focusing based on iteration and ray tracing

Science.gov (United States)

Li, Xiaotong; Chen, Hao; Wang, Weiwei; Ruan, Wangchao; Zhang, Luwei; Cen, Zhaofeng

2017-05-01

Self-focusing is observed in nonlinear materials owing to the interaction between laser and matter when laser beam propagates. Some of numerical simulation strategies such as the beam propagation method (BPM) based on nonlinear Schrödinger equation and ray tracing method based on Fermat's principle have applied to simulate the self-focusing process. In this paper we present an iteration nonlinear ray tracing method in that the nonlinear material is also cut into massive slices just like the existing approaches, but instead of paraxial approximation and split-step Fourier transform, a large quantity of sampled real rays are traced step by step through the system with changing refractive index and laser intensity by iteration. In this process a smooth treatment is employed to generate a laser density distribution at each slice to decrease the error caused by the under-sampling. The characteristics of this method is that the nonlinear refractive indices of the points on current slice are calculated by iteration so as to solve the problem of unknown parameters in the material caused by the causal relationship between laser intensity and nonlinear refractive index. Compared with the beam propagation method, this algorithm is more suitable for engineering application with lower time complexity, and has the calculation capacity for numerical simulation of self-focusing process in the systems including both of linear and nonlinear optical media. If the sampled rays are traced with their complex amplitudes and light paths or phases, it will be possible to simulate the superposition effects of different beam. At the end of the paper, the advantages and disadvantages of this algorithm are discussed.

Lattice Boltzmann model for high-order nonlinear partial differential equations.

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂_{t}ϕ+∑_{k=1}^{m}α_{k}∂_{x}^{k}Π_{k}(ϕ)=0 (1≤k≤m≤6), α_{k} are constant coefficients, Π_{k}(ϕ) are some known differential functions of ϕ. As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K(n,n)-Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009)1672-179910.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009)PHYADX0378-437110.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Lattice Boltzmann model for high-order nonlinear partial differential equations

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂tϕ +∑k=1mαk∂xkΠk(ϕ ) =0 (1 ≤k ≤m ≤6 ), αk are constant coefficients, Πk(ϕ ) are some known differential functions of ϕ . As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K (n ,n ) -Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009), 10.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009), 10.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Analytical and Numerical solutions of a nonlinear alcoholism model via variable-order fractional differential equations

Science.gov (United States)

Gómez-Aguilar, J. F.

2018-03-01

In this paper, we analyze an alcoholism model which involves the impact of Twitter via Liouville-Caputo and Atangana-Baleanu-Caputo fractional derivatives with constant- and variable-order. Two fractional mathematical models are considered, with and without delay. Special solutions using an iterative scheme via Laplace and Sumudu transform were obtained. We studied the uniqueness and existence of the solutions employing the fixed point postulate. The generalized model with variable-order was solved numerically via the Adams method and the Adams-Bashforth-Moulton scheme. Stability and convergence of the numerical solutions were presented in details. Numerical examples of the approximate solutions are provided to show that the numerical methods are computationally efficient. Therefore, by including both the fractional derivatives and finite time delays in the alcoholism model studied, we believe that we have established a more complete and more realistic indicator of alcoholism model and affect the spread of the drinking.

Numerical Methods for Nonlinear PDEs in Finance

DEFF Research Database (Denmark)

Mashayekhi, Sima

Nonlinear Black-Scholes equations arise from considering parameters such as feedback and illiquid markets eects or large investor preferences, volatile portfolio and nontrivial transaction costs into option pricing models to have more accurate option price. Here some nite dierence schemes have be...

Enzymatic Synthesis of Ampicillin: Nonlinear Modeling, Kinetics Estimation, and Adaptive Control

Directory of Open Access Journals (Sweden)

Monica Roman

2012-01-01

Full Text Available Nowadays, the use of advanced control strategies in biotechnology is quite low. A main reason is the lack of quality of the data, and the fact that more sophisticated control strategies must be based on a model of the dynamics of bioprocesses. The nonlinearity of the bioprocesses and the absence of cheap and reliable instrumentation require an enhanced modeling effort and identification strategies for the kinetics. The present work approaches modeling and control strategies for the enzymatic synthesis of ampicillin that is carried out inside a fed-batch bioreactor. First, a nonlinear dynamical model of this bioprocess is obtained by using a novel modeling procedure for biotechnology: the bond graph methodology. Second, a high gain observer is designed for the estimation of the imprecisely known kinetics of the synthesis process. Third, by combining an exact linearizing control law with the on-line estimation kinetics algorithm, a nonlinear adaptive control law is designed. The case study discussed shows that a nonlinear feedback control strategy applied to the ampicillin synthesis bioprocess can cope with disturbances, noisy measurements, and parametric uncertainties. Numerical simulations performed with MATLAB environment are included in order to test the behavior and the performances of the proposed estimation and control strategies.

A New Method to Solve Numeric Solution of Nonlinear Dynamic System

Directory of Open Access Journals (Sweden)

Min Hu

2016-01-01

Full Text Available It is well known that the cubic spline function has advantages of simple forms, good convergence, approximation, and second-order smoothness. A particular class of cubic spline function is constructed and an effective method to solve the numerical solution of nonlinear dynamic system is proposed based on the cubic spline function. Compared with existing methods, this method not only has high approximation precision, but also avoids the Runge phenomenon. The error analysis of several methods is given via two numeric examples, which turned out that the proposed method is a much more feasible tool applied to the engineering practice.

Nonlinear Modeling by Assembling Piecewise Linear Models

Science.gov (United States)

Yao, Weigang; Liou, Meng-Sing

2013-01-01

To preserve nonlinearity of a full order system over a parameters range of interest, we propose a simple modeling approach by assembling a set of piecewise local solutions, including the first-order Taylor series terms expanded about some sampling states. The work by Rewienski and White inspired our use of piecewise linear local solutions. The assembly of these local approximations is accomplished by assigning nonlinear weights, through radial basis functions in this study. The efficacy of the proposed procedure is validated for a two-dimensional airfoil moving at different Mach numbers and pitching motions, under which the flow exhibits prominent nonlinear behaviors. All results confirm that our nonlinear model is accurate and stable for predicting not only aerodynamic forces but also detailed flowfields. Moreover, the model is robustness-accurate for inputs considerably different from the base trajectory in form and magnitude. This modeling preserves nonlinearity of the problems considered in a rather simple and accurate manner.

Mathematical and numerical study of non-linear models used in plasma physics

International Nuclear Information System (INIS)

Ebrard, G.

2005-12-01

We study the interaction of several crossing beams with a plasma in the Laser-Megajoule context. We start from Euler-Maxwell. The formal asymptotic is the Zakharov system. For simplified systems of Klein-Gordon-wave type, we justify an approximation by a Zakharov equation for solutions of large amplitude. We construct a new system that simulates the interaction of 2 beams and present a whole hierarchy of models. We introduce a numerical scheme using the known results on Zakharov-wave equations which are valid for short pulses. We give a scheme which eliminate the backscattering wave. We give some numerical results. Finally, we do several numerical simulations of laser-plasma interaction for the initial value problem and the boundary value problem. (author)

Mathematical models of non-linear phenomena, processes and systems: from molecular scale to planetary atmosphere

CERN Document Server

2013-01-01

This book consists of twenty seven chapters, which can be divided into three large categories: articles with the focus on the mathematical treatment of non-linear problems, including the methodologies, algorithms and properties of analytical and numerical solutions to particular non-linear problems; theoretical and computational studies dedicated to the physics and chemistry of non-linear micro-and nano-scale systems, including molecular clusters, nano-particles and nano-composites; and, papers focused on non-linear processes in medico-biological systems, including mathematical models of ferments, amino acids, blood fluids and polynucleic chains.

Space-time least-squares Petrov-Galerkin projection in nonlinear model reduction.

Energy Technology Data Exchange (ETDEWEB)

Choi, Youngsoo [Sandia National Laboratories (SNL-CA), Livermore, CA (United States). Extreme-scale Data Science and Analytics Dept.; Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Carlberg, Kevin Thomas [Sandia National Laboratories (SNL-CA), Livermore, CA (United States). Extreme-scale Data Science and Analytics Dept.

2017-09-01

Our work proposes a space-time least-squares Petrov-Galerkin (ST-LSPG) projection method for model reduction of nonlinear dynamical systems. In contrast to typical nonlinear model-reduction methods that first apply Petrov-Galerkin projection in the spatial dimension and subsequently apply time integration to numerically resolve the resulting low-dimensional dynamical system, the proposed method applies projection in space and time simultaneously. To accomplish this, the method first introduces a low-dimensional space-time trial subspace, which can be obtained by computing tensor decompositions of state-snapshot data. The method then computes discrete-optimal approximations in this space-time trial subspace by minimizing the residual arising after time discretization over all space and time in a weighted ℓ²-norm. This norm can be de ned to enable complexity reduction (i.e., hyper-reduction) in time, which leads to space-time collocation and space-time GNAT variants of the ST-LSPG method. Advantages of the approach relative to typical spatial-projection-based nonlinear model reduction methods such as Galerkin projection and least-squares Petrov-Galerkin projection include: (1) a reduction of both the spatial and temporal dimensions of the dynamical system, (2) the removal of spurious temporal modes (e.g., unstable growth) from the state space, and (3) error bounds that exhibit slower growth in time. Numerical examples performed on model problems in fluid dynamics demonstrate the ability of the method to generate orders-of-magnitude computational savings relative to spatial-projection-based reduced-order models without sacrificing accuracy.

Adaptive regression for modeling nonlinear relationships

CERN Document Server

Knafl, George J

2016-01-01

This book presents methods for investigating whether relationships are linear or nonlinear and for adaptively fitting appropriate models when they are nonlinear. Data analysts will learn how to incorporate nonlinearity in one or more predictor variables into regression models for different types of outcome variables. Such nonlinear dependence is often not considered in applied research, yet nonlinear relationships are common and so need to be addressed. A standard linear analysis can produce misleading conclusions, while a nonlinear analysis can provide novel insights into data, not otherwise possible. A variety of examples of the benefits of modeling nonlinear relationships are presented throughout the book. Methods are covered using what are called fractional polynomials based on real-valued power transformations of primary predictor variables combined with model selection based on likelihood cross-validation. The book covers how to formulate and conduct such adaptive fractional polynomial modeling in the s...

The spectral cell method in nonlinear earthquake modeling

Science.gov (United States)

Giraldo, Daniel; Restrepo, Doriam

2017-12-01

This study examines the applicability of the spectral cell method (SCM) to compute the nonlinear earthquake response of complex basins. SCM combines fictitious-domain concepts with the spectral-version of the finite element method to solve the wave equations in heterogeneous geophysical domains. Nonlinear behavior is considered by implementing the Mohr-Coulomb and Drucker-Prager yielding criteria. We illustrate the performance of SCM with numerical examples of nonlinear basins exhibiting physically and computationally challenging conditions. The numerical experiments are benchmarked with results from overkill solutions, and using MIDAS GTS NX, a finite element software for geotechnical applications. Our findings show good agreement between the two sets of results. Traditional spectral elements implementations allow points per wavelength as low as PPW = 4.5 for high-order polynomials. Our findings show that in the presence of nonlinearity, high-order polynomials (p ≥ 3) require mesh resolutions above of PPW ≥ 10 to ensure displacement errors below 10%.

Finite Element Analysis of Biot’s Consolidation with a Coupled Nonlinear Flow Model

Directory of Open Access Journals (Sweden)

Yue-bao Deng

2016-01-01

Full Text Available A nonlinear flow relationship, which assumes that the fluid flow in the soil skeleton obeys the Hansbo non-Darcian flow and that the coefficient of permeability changes with void ratio, was incorporated into Biot’s general consolidation theory for a consolidation simulation of normally consolidated soft ground with or without vertical drains. The governing equations with the coupled nonlinear flow model were presented first for the force equilibrium condition and then for the continuity condition. Based on the weighted residual method, the finite element (FE formulations were then derived, and an existing FE program was modified accordingly to take the nonlinear flow model into consideration. Comparative analyses using established theoretical solutions and numerical solutions were completed, and the results were satisfactory. On this basis, we investigated the effect of the coupled nonlinear flow on consolidation development.

A 2D nonlinear multiring model for blood flow in large elastic arteries

Science.gov (United States)

Ghigo, Arthur R.; Fullana, Jose-Maria; Lagrée, Pierre-Yves

2017-12-01

In this paper, we propose a two-dimensional nonlinear ;multiring; model to compute blood flow in axisymmetric elastic arteries. This model is designed to overcome the numerical difficulties of three-dimensional fluid-structure interaction simulations of blood flow without using the over-simplifications necessary to obtain one-dimensional blood flow models. This multiring model is derived by integrating over concentric rings of fluid the simplified long-wave Navier-Stokes equations coupled to an elastic model of the arterial wall. The resulting system of balance laws provides a unified framework in which both the motion of the fluid and the displacement of the wall are dealt with simultaneously. The mathematical structure of the multiring model allows us to use a finite volume method that guarantees the conservation of mass and the positivity of the numerical solution and can deal with nonlinear flows and large deformations of the arterial wall. We show that the finite volume numerical solution of the multiring model provides at a reasonable computational cost an asymptotically valid description of blood flow velocity profiles and other averaged quantities (wall shear stress, flow rate, ...) in large elastic and quasi-rigid arteries. In particular, we validate the multiring model against well-known solutions such as the Womersley or the Poiseuille solutions as well as against steady boundary layer solutions in quasi-rigid constricted and expanded tubes.

Single-step emulation of nonlinear fiber-optic link with gaussian mixture model

DEFF Research Database (Denmark)

Borkowski, Robert; Doberstein, Andy; Haisch, Hansjörg

2015-01-01

We use a fast and low-complexity statistical signal processing method to emulate nonlinear noise in fiber links. The proposed emulation technique stands in good agreement with the numerical NLSE simulation for 32 Gbaud DP-16QAM nonlinear transmission.......We use a fast and low-complexity statistical signal processing method to emulate nonlinear noise in fiber links. The proposed emulation technique stands in good agreement with the numerical NLSE simulation for 32 Gbaud DP-16QAM nonlinear transmission....

Pressure fluctuation prediction of a model pump turbine at no load opening by a nonlinear k-ε turbulence model

International Nuclear Information System (INIS)

Liu, J T; Zuo, Z G; Liu, S H; Wu, Y L

2014-01-01

In this paper, a new nonlinear k-ε turbulence model based on RNG k-ε turbulence model and Wilcox's k-ω turbulence model was proposed to simulate the unsteady flow and to predict the pressure fluctuation through a model pump turbine for engineering application. Calculations on a curved rectangular duct proved that the nonlinear k-ε turbulence model is applicable for high pressure gradient flows and large curvature flows. The numerically predicted relative pressure amplitude (peak to peak) in time domain to the pump turbine head at no load condition is very close to the experimental data. It is indicated that the prediction of the pressure fluctuation is valid by the present nonlinear k-ε method. The high pressure fluctuation in this area is the main issue for pump turbine design, especially at high head condition

Efficient Hybrid-Spectral Model for Fully Nonlinear Numerical Wave Tank

DEFF Research Database (Denmark)

Christiansen, Torben; Bingham, Harry B.; Engsig-Karup, Allan Peter

2013-01-01

A new hybrid-spectral solution strategy is proposed for the simulation of the fully nonlinear free surface equations based on potential flow theory. A Fourier collocation method is adopted horisontally for the discretization of the free surface equations. This is combined with a modal Chebyshev Tau...... method in the vertical for the discretization of the Laplace equation in the fluid domain, which yields a sparse and spectrally accurate Dirichletto-Neumann operator. The Laplace problem is solved with an efficient Defect Correction method preconditioned with a spectral discretization of the linearised...... wave problem, ensuring fast convergence and optimal scaling with the problem size. Preliminary results for very nonlinear waves show expected convergence rates and a clear advantage of using spectral schemes....

Nonlinear modelling of polymer electrolyte membrane fuel cell stack using nonlinear cancellation technique

Energy Technology Data Exchange (ETDEWEB)

Barus, R. P. P., E-mail: rismawan.ppb@gmail.com [Engineering Physics, Faculty of Industrial Technology, Institut Teknologi Bandung, Jalan Ganesa 10 Bandung and Centre for Material and Technical Product, Jalan Sangkuriang No. 14 Bandung (Indonesia); Tjokronegoro, H. A.; Leksono, E. [Engineering Physics, Faculty of Industrial Technology, Institut Teknologi Bandung, Jalan Ganesa 10 Bandung (Indonesia); Ismunandar [Chemistry Study, Faculty of Mathematics and Science, Institut Teknologi Bandung, Jalan Ganesa 10 Bandung (Indonesia)

2014-09-25

Fuel cells are promising new energy conversion devices that are friendly to the environment. A set of control systems are required in order to operate a fuel cell based power plant system optimally. For the purpose of control system design, an accurate fuel cell stack model in describing the dynamics of the real system is needed. Currently, linear model are widely used for fuel cell stack control purposes, but it has limitations in narrow operation range. While nonlinear models lead to nonlinear control implemnetation whos more complex and hard computing. In this research, nonlinear cancellation technique will be used to transform a nonlinear model into a linear form while maintaining the nonlinear characteristics. The transformation is done by replacing the input of the original model by a certain virtual input that has nonlinear relationship with the original input. Then the equality of the two models is tested by running a series of simulation. Input variation of H2, O2 and H2O as well as disturbance input I (current load) are studied by simulation. The error of comparison between the proposed model and the original nonlinear model are less than 1 %. Thus we can conclude that nonlinear cancellation technique can be used to represent fuel cell nonlinear model in a simple linear form while maintaining the nonlinear characteristics and therefore retain the wide operation range.

Nonlinear modelling of polymer electrolyte membrane fuel cell stack using nonlinear cancellation technique

International Nuclear Information System (INIS)

Barus, R. P. P.; Tjokronegoro, H. A.; Leksono, E.; Ismunandar

2014-01-01

Fuel cells are promising new energy conversion devices that are friendly to the environment. A set of control systems are required in order to operate a fuel cell based power plant system optimally. For the purpose of control system design, an accurate fuel cell stack model in describing the dynamics of the real system is needed. Currently, linear model are widely used for fuel cell stack control purposes, but it has limitations in narrow operation range. While nonlinear models lead to nonlinear control implemnetation whos more complex and hard computing. In this research, nonlinear cancellation technique will be used to transform a nonlinear model into a linear form while maintaining the nonlinear characteristics. The transformation is done by replacing the input of the original model by a certain virtual input that has nonlinear relationship with the original input. Then the equality of the two models is tested by running a series of simulation. Input variation of H2, O2 and H2O as well as disturbance input I (current load) are studied by simulation. The error of comparison between the proposed model and the original nonlinear model are less than 1 %. Thus we can conclude that nonlinear cancellation technique can be used to represent fuel cell nonlinear model in a simple linear form while maintaining the nonlinear characteristics and therefore retain the wide operation range

Nonlinear behavior analysis of split-winding dry-type transformer using a new star model and a coupled field-circuit approach

Directory of Open Access Journals (Sweden)

Azizian Davood

2016-12-01

Full Text Available Regarding the importance of short circuit and inrush current simulations in the split-winding transformer, a novel nonlinear equivalent circuit is introduced in this paper for nonlinear simulation of this transformer. The equivalent circuit is extended using the nonlinear inductances. Employing a numerical method, leakage and magnetizing inductances in the split-winding transformer are extracted and the nonlinear model inductances are estimated using these inductances. The introduced model is validated and using this nonlinear model, inrush and short-circuit currents are calculated. It has been seen that the introduced model is valid and suitable for simulations of the split-winding transformer due to various loading conditions. Finally, the effects of nonlinearity of the model inductances are discussed in the following.

Eddy current modeling in linear and nonlinear multifilamentary composite materials

Science.gov (United States)

Menana, Hocine; Farhat, Mohamad; Hinaje, Melika; Berger, Kevin; Douine, Bruno; Lévêque, Jean

2018-04-01

In this work, a numerical model is developed for a rapid computation of eddy currents in composite materials, adaptable for both carbon fiber reinforced polymers (CFRPs) for NDT applications and multifilamentary high temperature superconductive (HTS) tapes for AC loss evaluation. The proposed model is based on an integro-differential formulation in terms of the electric vector potential in the frequency domain. The high anisotropy and the nonlinearity of the considered materials are easily handled in the frequency domain.

Nonlinear saturation controller for vibration supersession of a nonlinear composite beam

Energy Technology Data Exchange (ETDEWEB)

Hamed, Y. S. [Menofia University, Menouf (Egypt); Amer, Y. A. [Zagazig University, Zagazig (Egypt)

2014-08-15

In this paper, a study for nonlinear saturation controller (NSC) is presented that used to suppress the vibration amplitude of a structural dynamic model simulating nonlinear composite beam at simultaneous sub-harmonic and internal resonance excitation. The absorber exploits the saturation phenomenon that is known to occur in dynamical systems with quadratic non-linearities of the feedback gain and a two-to-one internal resonance. The analytical solution for the system and the nonlinear saturation controller are obtained using method of multiple time scales perturbation up to the second order approximation. All possible resonance cases were extracted at this approximation order and studied numerically. The stability of the system at the worst resonance case (Ω = 2ω{sub s} and ω{sub s} =2ω{sub C}) is investigated using both frequency response equations and phase-plane trajectories. The effects of different parameters on the system and the controller are studied numerically. The effect of some types of controller on the system is investigated numerically. The simulation results are achieved using Matlab and Maple programs.

Numerical modelling of multimode fibre-optic communication lines

Energy Technology Data Exchange (ETDEWEB)

Sidelnikov, O S; Fedoruk, M P [Novosibirsk State University, Novosibirsk (Russian Federation); Sygletos, S; Ferreira, F [Aston University, England, Birmingham, B4 7ET (United Kingdom)

2016-01-31

The results of numerical modelling of nonlinear propagation of an optical signal in multimode fibres with a small differential group delay are presented. It is found that the dependence of the error vector magnitude (EVM) on the differential group delay can be reduced by increasing the number of ADC samples per symbol in the numerical implementation of the differential group delay compensation algorithm in the receiver. The possibility of using multimode fibres with a small differential group delay for data transmission in modern digital communication systems is demonstrated. It is shown that with increasing number of modes the strong coupling regime provides a lower EVM level than the weak coupling one. (fibre-optic communication lines)

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

International Nuclear Information System (INIS)

Lin-Jie, Chen; Chang-Feng, Ma

2010-01-01

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form u t + αuu x + βu n u x + γu xx + δu xxx + ζu xxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman–Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions. (general)

Numerical solution of continuous-time DSGE models under Poisson uncertainty

DEFF Research Database (Denmark)

Posch, Olaf; Trimborn, Timo

We propose a simple and powerful method for determining the transition process in continuous-time DSGE models under Poisson uncertainty numerically. The idea is to transform the system of stochastic differential equations into a system of functional differential equations of the retarded type. We...... classes of models. We illustrate the algorithm simulating both the stochastic neoclassical growth model and the Lucas model under Poisson uncertainty which is motivated by the Barro-Rietz rare disaster hypothesis. We find that, even for non-linear policy functions, the maximum (absolute) error is very...

Spurious Numerical Solutions Of Differential Equations

Science.gov (United States)

Lafon, A.; Yee, H. C.

1995-01-01

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Model reduction and frequency residuals for a robust estimation of nonlinearities in subspace identification

Science.gov (United States)

De Filippis, G.; Noël, J. P.; Kerschen, G.; Soria, L.; Stephan, C.

2017-09-01

The introduction of the frequency-domain nonlinear subspace identification (FNSI) method in 2013 constitutes one in a series of recent attempts toward developing a realistic, first-generation framework applicable to complex structures. If this method showed promising capabilities when applied to academic structures, it is still confronted with a number of limitations which needs to be addressed. In particular, the removal of nonphysical poles in the identified nonlinear models is a distinct challenge. In the present paper, it is proposed as a first contribution to operate directly on the identified state-space matrices to carry out spurious pole removal. A modal-space decomposition of the state and output matrices is examined to discriminate genuine from numerical poles, prior to estimating the extended input and feedthrough matrices. The final state-space model thus contains physical information only and naturally leads to nonlinear coefficients free of spurious variations. Besides spurious variations due to nonphysical poles, vibration modes lying outside the frequency band of interest may also produce drifts of the nonlinear coefficients. The second contribution of the paper is to include residual terms, accounting for the existence of these modes. The proposed improved FNSI methodology is validated numerically and experimentally using a full-scale structure, the Morane-Saulnier Paris aircraft.

Finite Element Model for Nonlinear Analysis of Reinforced Concrete Beams and Plane Frames

Directory of Open Access Journals (Sweden)

R.S.B. STRAMANDINOLI

Full Text Available Abstract In this work, a two-dimensional finite element (FE model for physical and geometric nonlinear analysis of reinforced concrete beams and plane frames, developed by the authors, is presented. The FE model is based on the Euler-Bernoulli Beam Theory, in which shear deformations are neglected. The bar elements have three nodes with a total of seven degrees of freedom. Three Gauss-points are utilized for the element integration, with the element section discretized into layers at each Gauss point (Fiber Model. It is assumed that concrete and reinforcing bars are perfectly bonded, and each section layer is assumed to be under a uniaxial stress-state. Nonlinear constitutive laws are utilized for both concrete and reinforcing steel layers, and a refined tension-stiffening model, developed by the authors, is included. The Total Lagrangean Formulation is adopted for geometric nonlinear consideration and several methods can be utilized to achieve equilibrium convergence of the nonlinear equations. The developed model is implemented into a computer program named ANEST/CA, which is validated by comparison with some tests on RC beams and plane frames, showing an excellent correlation between numerical and experimental results.

A non-linear procedure for the numerical analysis of crack development in beams failing in shear

Directory of Open Access Journals (Sweden)

P. Bernardi

2016-01-01

Full Text Available In this work, a consistent formulation for the representation of concrete behavior before and after cracking has been implemented into a non-linear model for the analysis of reinforced concrete structures, named 2D-PARC. Several researches have indeed pointed out that the adoption of an effective modeling for concrete, combined with an accurate failure criterion, is crucial for the correct prediction of the structural behavior, not only in terms of failure load, but also with reference to a realistic representation of crack initiation and development. This last aspect is particularly relevant at serviceability conditions in order to verify the fulfillment of structural requirements provided by Design Codes, which limit the maximum crack width due to appearance and durability issues. In more details, a constitutive model originally proposed by Ottosen and based on non-linear elasticity has been here incorporated into 2D-PARC in order to improve the numerical efficiency of the adopted algorithm, providing at the same time an accurate prediction of the structural response. The effectiveness of this procedure has been verified against significant experimental results available in the technical literature and relative to reinforced concrete beams without stirrups failing in shear, which represent a problem of great theoretical and practical importance in the field of structural engineering. Numerical results have been compared to experimental evidences not only in terms of global structural response (i.e. applied load vs. midspan deflection, but also in terms of crack pattern evolution and maximum crack widths.

Non-perturbative aspects of nonlinear sigma models

Energy Technology Data Exchange (ETDEWEB)

Flore, Raphael

2012-12-07

The aim of this thesis was the study and further development of non-perturbative methods of quantum field theory by means of their application to nonlinear sigma models. While a large part of the physical phenomena of quantum field theory can be successfully predicted by the perturbation theory, some aspects in the region of large coupling strengths are not definitively understood and require suited non-perturbative methods for its analysis. This thesis is concentrated on two approaches, the numerical treatment of field theories on discrete space-time lattices and the functional renormalization group (FRG) as description of the renormalization flux of effective actions. Considerations of the nonlinear O(N) models have shown that for the correct analysis of the critical properties in the framework of the FRG an approach must be chosen, which contained fourth-derivation orders. For this a covariant formalism was developed, which is based on a background-field expansion and the development of a heat kernel. Apart from a destabilizing coupling the results suggest a nontrivial fixed point and by this a non-perturbative renormalizability of these models. The resulting flow diagrams were finally still compared with the results of a numerical analysis of the renormalization flow by means of the Monte-Carlo renormalization group, and hereby qualitative agreement was found. Furthermore an alternative formulation of the FRG in phase-space coordinates was studied and their consistency tested on simple examples. Beyond this an alternative expansion of the effective action in orders of the canonical momenta was applied to the nonlinear O(N) models with the result of a stable non-trivial fixed point, the critical properties of which however show not the expected N-dependence. By means of the FRG finally still the renormalization of topological operators was studied by means of the winding number of the O(3){approx_equal}CP{sup 1} model. By the generalization of the topological

Non-perturbative aspects of nonlinear sigma models

International Nuclear Information System (INIS)

Flore, Raphael

2012-01-01

The aim of this thesis was the study and further development of non-perturbative methods of quantum field theory by means of their application to nonlinear sigma models. While a large part of the physical phenomena of quantum field theory can be successfully predicted by the perturbation theory, some aspects in the region of large coupling strengths are not definitively understood and require suited non-perturbative methods for its analysis. This thesis is concentrated on two approaches, the numerical treatment of field theories on discrete space-time lattices and the functional renormalization group (FRG) as description of the renormalization flux of effective actions. Considerations of the nonlinear O(N) models have shown that for the correct analysis of the critical properties in the framework of the FRG an approach must be chosen, which contained fourth-derivation orders. For this a covariant formalism was developed, which is based on a background-field expansion and the development of a heat kernel. Apart from a destabilizing coupling the results suggest a nontrivial fixed point and by this a non-perturbative renormalizability of these models. The resulting flow diagrams were finally still compared with the results of a numerical analysis of the renormalization flow by means of the Monte-Carlo renormalization group, and hereby qualitative agreement was found. Furthermore an alternative formulation of the FRG in phase-space coordinates was studied and their consistency tested on simple examples. Beyond this an alternative expansion of the effective action in orders of the canonical momenta was applied to the nonlinear O(N) models with the result of a stable non-trivial fixed point, the critical properties of which however show not the expected N-dependence. By means of the FRG finally still the renormalization of topological operators was studied by means of the winding number of the O(3)≅CP 1 model. By the generalization of the topological operator and the

Modeling, numerical simulation, and nonlinear dynamic behavior analysis of PV microgrid-connected inverter with capacitance catastrophe

Science.gov (United States)

Li, Sichen; Liao, Zhixian; Luo, Xiaoshu; Wei, Duqu; Jiang, Pinqun; Jiang, Qinghong

2018-02-01

The value of the output capacitance (C) should be carefully considered when designing a photovoltaic (PV) inverter since it can cause distortion in the working state of the circuit, and the circuit produces nonlinear dynamic behavior. According to Kirchhoff’s laws and the characteristics of an ideal operational amplifier for a strict piecewise linear state equation, a circuit simulation model is constructed to study the system parameters (time, C) for the current passing through an inductor with an inductance of L and the voltage across the capacitor with a capacitance of C. The developed simulation model uses Runge-Kutta methods to solve the state equations. This study focuses on predicting the fault of the circuit from the two aspects of the harmonic distortion and simulation results. Moreover, the presented model is also used to research the working state of the system in the case of a load capacitance catastrophe. The nonlinear dynamic behaviors in the inverter are simulated and verified.

Stationary nonlinear Airy beams

International Nuclear Information System (INIS)

Lotti, A.; Faccio, D.; Couairon, A.; Papazoglou, D. G.; Panagiotopoulos, P.; Tzortzakis, S.; Abdollahpour, D.

2011-01-01

We demonstrate the existence of an additional class of stationary accelerating Airy wave forms that exist in the presence of third-order (Kerr) nonlinearity and nonlinear losses. Numerical simulations and experiments, in agreement with the analytical model, highlight how these stationary solutions sustain the nonlinear evolution of Airy beams. The generic nature of the Airy solution allows extension of these results to other settings, and a variety of applications are suggested.

On a New Method for Computing the Numerical Solution of Systems of Nonlinear Equations

Directory of Open Access Journals (Sweden)

H. Montazeri

2012-01-01

Full Text Available We consider a system of nonlinear equations F(x=0. A new iterative method for solving this problem numerically is suggested. The analytical discussions of the method are provided to reveal its sixth order of convergence. A discussion on the efficiency index of the contribution with comparison to the other iterative methods is also given. Finally, numerical tests illustrate the theoretical aspects using the programming package Mathematica.

Generalized Nonlinear Yule Models

Science.gov (United States)

Lansky, Petr; Polito, Federico; Sacerdote, Laura

2016-11-01

With the aim of considering models related to random graphs growth exhibiting persistent memory, we propose a fractional nonlinear modification of the classical Yule model often studied in the context of macroevolution. Here the model is analyzed and interpreted in the framework of the development of networks such as the World Wide Web. Nonlinearity is introduced by replacing the linear birth process governing the growth of the in-links of each specific webpage with a fractional nonlinear birth process with completely general birth rates. Among the main results we derive the explicit distribution of the number of in-links of a webpage chosen uniformly at random recognizing the contribution to the asymptotics and the finite time correction. The mean value of the latter distribution is also calculated explicitly in the most general case. Furthermore, in order to show the usefulness of our results, we particularize them in the case of specific birth rates giving rise to a saturating behaviour, a property that is often observed in nature. The further specialization to the non-fractional case allows us to extend the Yule model accounting for a nonlinear growth.

CROSS DIFFUSION AND NONLINEAR DIFFUSION PREVENTING BLOW UP IN THE KELLER–SEGEL MODEL

KAUST Repository

CARRILLO, JOSÉ ANTONIO

2012-12-01

A parabolic-parabolic (Patlak-)Keller-Segel model in up to three space dimensions with nonlinear cell diffusion and an additional nonlinear cross-diffusion term is analyzed. The main feature of this model is that there exists a new entropy functional, yielding gradient estimates for the cell density and chemical concentration. For arbitrarily small cross-diffusion coefficients and for suitable exponents of the nonlinear diffusion terms, the global-in-time existence of weak solutions is proved, thus preventing finite-time blow up of the cell density. The global existence result also holds for linear and fast diffusion of the cell density in a certain parameter range in three dimensions. Furthermore, we show L∞ bounds for the solutions to the parabolic-elliptic system. Sufficient conditions leading to the asymptotic stability of the constant steady state are given for a particular choice of the nonlinear diffusion exponents. Numerical experiments in two and three space dimensions illustrate the theoretical results. © 2012 World Scientific Publishing Company.

Study of intermittent bifurcations and chaos in boost PFC converters by nonlinear discrete models

International Nuclear Information System (INIS)

Zhang Hao; Ma Xikui; Xue Bianling; Liu Weizeng

2005-01-01

This paper mainly deals with nonlinear phenomena like intermittent bifurcations and chaos in boost PFC converters under peak-current control mode. Two nonlinear models in the form of discrete maps are derived to describe precisely the nonlinear dynamics of boost PFC converters from two points of view, i.e., low- and high-frequency regimes. Based on the presented discrete models, both the evolution of intermittent behavior and the periodicity of intermittency are investigated in detail from the fast and slow-scale aspects, respectively. Numerical results show that the occurrence of intermittent bifurcations and chaos with half one line period is one of the most distinguished dynamical characteristics. Finally, we make some instructive conclusions, which prove to be helpful in improving the performances of practical circuits

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models.

Science.gov (United States)

Pozo, Carlos; Marín-Sanguino, Alberto; Alves, Rui; Guillén-Gosálbez, Gonzalo; Jiménez, Laureano; Sorribas, Albert

2011-08-25

Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Directory of Open Access Journals (Sweden)

Sorribas Albert

2011-08-01

Full Text Available Abstract Background Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Relaxation and Numerical Approximation of a Two-Fluid Two-Pressure Diphasic Model

International Nuclear Information System (INIS)

Ambroso, A.; Chalons, Ch.; Galie, Th.; Chalons, Ch.; Coquel, F.; Coquel, F.

2009-01-01

This paper is concerned with the numerical approximation of the solutions of a two-fluid two-pressure model used in the modelling of two-phase flows. We present a relaxation strategy for easily dealing with both the nonlinearities associated with the pressure laws and the nonconservative terms that are inherently present in the set of convective equations and that couple the two phases. In particular, the proposed approximate Riemann solver is given by explicit formulas, preserves the natural phase space, and exactly captures the coupling waves between the two phases. Numerical evidences are given to corroborate the validity of our approach. (authors)

Linear and nonlinear stability analysis in BWRs applying a reduced order model

Energy Technology Data Exchange (ETDEWEB)

Olvera G, O. A.; Espinosa P, G.; Prieto G, A., E-mail: omar_olverag@hotmail.com [Universidad Autonoma Metropolitana, Unidad Iztapalapa, San Rafael Atlixco No. 186, Col. Vicentina, 09340 Ciudad de Mexico (Mexico)

2016-09-15

Boiling Water Reactor (BWR) stability studies are generally conducted through nonlinear reduced order models (Rom) employing various techniques such as bifurcation analysis and time domain numerical integration. One of those models used for these studies is the March-Leuba Rom. Such model represents qualitatively the dynamic behavior of a BWR through a one-point reactor kinetics, a one node representation of the heat transfer process in fuel, and a two node representation of the channel Thermal hydraulics to account for the void reactivity feedback. Here, we study the effect of this higher order model on the overall stability of the BWR. The change in the stability boundaries is determined by evaluating the eigenvalues of the Jacobian matrix. The nonlinear model is also integrated numerically to show that in the nonlinear region, the system evolves to stable limit cycles when operating close to the stability boundary. We also applied a new technique based on the Empirical Mode Decomposition (Emd) to estimate a parameter linked with stability in a BWR. This instability parameter is not exactly the classical Decay Ratio (Dr), but it will be linked with it. The proposed method allows decomposing the analyzed signal in different levels or mono-component functions known as intrinsic mode functions (Imf). One or more of these different modes can be associated to the instability problem in BWRs. By tracking the instantaneous frequencies (calculated through Hilbert Huang Transform (HHT) and the autocorrelation function (Acf) of the Imf linked to instability. The estimation of the proposed parameter can be achieved. The current methodology was validated with simulated signals of the studied model. (Author)

Linear and nonlinear stability analysis in BWRs applying a reduced order model

International Nuclear Information System (INIS)

Olvera G, O. A.; Espinosa P, G.; Prieto G, A.

2016-09-01

Boiling Water Reactor (BWR) stability studies are generally conducted through nonlinear reduced order models (Rom) employing various techniques such as bifurcation analysis and time domain numerical integration. One of those models used for these studies is the March-Leuba Rom. Such model represents qualitatively the dynamic behavior of a BWR through a one-point reactor kinetics, a one node representation of the heat transfer process in fuel, and a two node representation of the channel Thermal hydraulics to account for the void reactivity feedback. Here, we study the effect of this higher order model on the overall stability of the BWR. The change in the stability boundaries is determined by evaluating the eigenvalues of the Jacobian matrix. The nonlinear model is also integrated numerically to show that in the nonlinear region, the system evolves to stable limit cycles when operating close to the stability boundary. We also applied a new technique based on the Empirical Mode Decomposition (Emd) to estimate a parameter linked with stability in a BWR. This instability parameter is not exactly the classical Decay Ratio (Dr), but it will be linked with it. The proposed method allows decomposing the analyzed signal in different levels or mono-component functions known as intrinsic mode functions (Imf). One or more of these different modes can be associated to the instability problem in BWRs. By tracking the instantaneous frequencies (calculated through Hilbert Huang Transform (HHT) and the autocorrelation function (Acf) of the Imf linked to instability. The estimation of the proposed parameter can be achieved. The current methodology was validated with simulated signals of the studied model. (Author)

The Lie-Poisson structure of integrable classical non-linear sigma models

International Nuclear Information System (INIS)

Bordemann, M.; Forger, M.; Schaeper, U.; Laartz, J.

1993-01-01

The canonical structure of classical non-linear sigma models on Riemannian symmetric spaces, which constitute the most general class of classical non-linear sigma models known to be integrable, is shown to be governed by a fundamental Poisson bracket relation that fits into the r-s-matrix formalism for non-ultralocal integrable models first discussed by Maillet. The matrices r and s are computed explicitly and, being field dependent, satisfy fundamental Poisson bracket relations of their own, which can be expressed in terms of a new numerical matrix c. It is proposed that all these Poisson brackets taken together are representation conditions for a new kind of algebra which, for this class of models, replaces the classical Yang-Baxter algebra governing the canonical structure of ultralocal models. The Poisson brackets for the transition matrices are also computed, and the notorious regularization problem associated with the definition of the Poisson brackets for the monodromy matrices is discussed. (orig.)

Simplified Qualitative Discrete Numerical Model to Determine Cracking Pattern in Brittle Materials by Means of Finite Element Method

Directory of Open Access Journals (Sweden)

J. Ochoa-Avendaño

2017-01-01

Full Text Available This paper presents the formulation, implementation, and validation of a simplified qualitative model to determine the crack path of solids considering static loads, infinitesimal strain, and plane stress condition. This model is based on finite element method with a special meshing technique, where nonlinear link elements are included between the faces of the linear triangular elements. The stiffness loss of some link elements represents the crack opening. Three experimental tests of bending beams are simulated, where the cracking pattern calculated with the proposed numerical model is similar to experimental result. The advantages of the proposed model compared to discrete crack approaches with interface elements can be the implementation simplicity, the numerical stability, and the very low computational cost. The simulation with greater values of the initial stiffness of the link elements does not affect the discontinuity path and the stability of the numerical solution. The exploded mesh procedure presented in this model avoids a complex nonlinear analysis and regenerative or adaptive meshes.

Numerical simulation of nonlinear dynamical systems driven by commutative noise

International Nuclear Information System (INIS)

Carbonell, F.; Biscay, R.J.; Jimenez, J.C.; Cruz, H. de la

2007-01-01

The local linearization (LL) approach has become an effective technique for the numerical integration of ordinary, random and stochastic differential equations. One of the reasons for this success is that the LL method achieves a convenient trade-off between numerical stability and computational cost. Besides, the LL method reproduces well the dynamics of nonlinear equations for which other classical methods fail. However, in the stochastic case, most of the reported works has been focused in Stochastic Differential Equations (SDE) driven by additive noise. This limits the applicability of the LL method since there is a number of interesting dynamics observed in equations with multiplicative noise. On the other hand, recent results show that commutative noise SDEs can be transformed into a random differential equation (RDE) by means of a random diffeomorfism (conjugacy). This paper takes advantages of such conjugacy property and the LL approach for defining a LL scheme for SDEs driven by commutative noise. The performance of the proposed method is illustrated by means of numerical simulations

On the Validation of a Numerical Model for the Analysis of Soil-Structure Interaction Problems

Directory of Open Access Journals (Sweden)

Jorge Luis Palomino Tamayo

Full Text Available Abstract Modeling and simulation of mechanical response of structures, relies on the use of computational models. Therefore, verification and validation procedures are the primary means of assessing accuracy, confidence and credibility in modeling. This paper is concerned with the validation of a three dimensional numerical model based on the finite element method suitable for the dynamic analysis of soil-structure interaction problems. The soil mass, structure, structure's foundation and the appropriate boundary conditions can be represented altogether in a single model by using a direct approach. The theory of porous media of Biot is used to represent the soil mass as a two-phase material which is considered to be fully saturated with water; meanwhile other parts of the system are treated as one-phase materials. Plasticity of the soil mass is the main source of non-linearity in the problem and therefore an iterative-incremental algorithm based on the Newton-Raphson procedure is used to solve the nonlinear equilibrium equations. For discretization in time, the Generalized Newmark-β method is used. The soil is represented by a plasticity-based, effective-stress constitutive model suitable for liquefaction. Validation of the present numerical model is done by comparing analytical and centrifuge test results of soil and soil-pile systems with those results obtained with the present numerical model. A soil-pile-structure interaction problem is also presented in order to shown the potentiality of the numerical tool.

Numerical Simulations of Self-Focused Pulses Using the Nonlinear Maxwell Equations

Science.gov (United States)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1994-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations. Abstract of a proposed paper for presentation at the meeting NONLINEAR OPTICS: Materials, Fundamentals, and Applications, Hyatt Regency Waikaloa, Waikaloa, Hawaii, July 24-29, 1994, Cosponsored by IEEE/Lasers and Electro-Optics Society and Optical Society of America

Dynamic Behavior for an SIRS Model with Nonlinear Incidence Rate and Treatment

Directory of Open Access Journals (Sweden)

Junhong Li

2013-01-01

Full Text Available This paper considers an SIRS model with nonlinear incidence rate and treatment. It is assumed that susceptible and infectious individuals have constant immigration rates. We investigate the existence of equilibrium and prove the global asymptotical stable results of the endemic equilibrium. We then obtained that the model undergoes a Hopf bifurcation and existences a limit cycle. Some numerical simulations are given to illustrate the analytical results.

Nonlinear models for autoregressive conditional heteroskedasticity

DEFF Research Database (Denmark)

Teräsvirta, Timo

This paper contains a brief survey of nonlinear models of autore- gressive conditional heteroskedasticity. The models in question are parametric nonlinear extensions of the original model by Engle (1982). After presenting the individual models, linearity testing and parameter estimation are discu...

Nonlinear electromechanical modelling and dynamical behavior analysis of a satellite reaction wheel

Science.gov (United States)

Aghalari, Alireza; Shahravi, Morteza

2017-12-01

The present research addresses the satellite reaction wheel (RW) nonlinear electromechanical coupling dynamics including dynamic eccentricity of brushless dc (BLDC) motor and gyroscopic effects, as well as dry friction of shaft-bearing joints (relative small slip) and bearing friction. In contrast to other studies, the rotational velocity of the flywheel is considered to be controllable, so it is possible to study the reaction wheel dynamical behavior in acceleration stages. The RW is modeled as a three-phases BLDC motor as well as flywheel with unbalances on a rigid shaft and flexible bearings. Improved Lagrangian dynamics for electromechanical systems is used to obtain the mathematical model of the system. The developed model can properly describe electromechanical nonlinear coupled dynamical behavior of the satellite RW. Numerical simulations show the effectiveness of the presented approach.

Model tests and numerical analyses on horizontal impedance functions of inclined single piles embedded in cohesionless soil

Science.gov (United States)

Goit, Chandra Shekhar; Saitoh, Masato

2013-03-01

Horizontal impedance functions of inclined single piles are measured experimentally for model soil-pile systems with both the effects of local soil nonlinearity and resonant characteristics. Two practical pile inclinations of 5° and 10° in addition to a vertical pile embedded in cohesionless soil and subjected to lateral harmonic pile head loadings for a wide range of frequencies are considered. Results obtained with low-to-high amplitude of lateral loadings on model soil-pile systems encased in a laminar shear box show that the local nonlinearities have a profound impact on the horizontal impedance functions of piles. Horizontal impedance functions of inclined piles are found to be smaller than the vertical pile and the values decrease as the angle of pile inclination increases. Distinct values of horizontal impedance functions are obtained for the `positive' and `negative' cycles of harmonic loadings, leading to asymmetric force-displacement relationships for the inclined piles. Validation of these experimental results is carried out through three-dimensional nonlinear finite element analyses, and the results from the numerical models are in good agreement with the experimental data. Sensitivity analyses conducted on the numerical models suggest that the consideration of local nonlinearity at the vicinity of the soil-pile interface influence the response of the soil-pile systems.

solveME: fast and reliable solution of nonlinear ME models

DEFF Research Database (Denmark)

Yang, Laurence; Ma, Ding; Ebrahim, Ali

2016-01-01

Background: Genome-scale models of metabolism and macromolecular expression (ME) significantly expand the scope and predictive capabilities of constraint-based modeling. ME models present considerable computational challenges: they are much (>30 times) larger than corresponding metabolic reconstr......Background: Genome-scale models of metabolism and macromolecular expression (ME) significantly expand the scope and predictive capabilities of constraint-based modeling. ME models present considerable computational challenges: they are much (>30 times) larger than corresponding metabolic...... reconstructions (M models), are multiscale, and growth maximization is a nonlinear programming (NLP) problem, mainly due to macromolecule dilution constraints. Results: Here, we address these computational challenges. We develop a fast and numerically reliable solution method for growth maximization in ME models...

2D and 3D numerical modeling of seismic waves from explosion sources

International Nuclear Information System (INIS)

McLaughlin, K.L.; Stevens, J.L.; Barker, T.G.; Shkoller, B.; Day, S.M.

1993-01-01

Over the last decade, nonlinear and linear 2D axisymmetric finite difference codes have been used in conjunction with far-field seismic Green's functions to simulate seismic waves from a variety of sources. In this paper we briefly review some of the results and conclusions that have resulted from numerical simulations and explosion modeling in support of treaty verification research at S-CUBED in the last decade. We then describe in more detail the results from two recent projects. Our goal is to provide a flavor for the kinds of problems that can be examined with numerical methods for modeling excitation of seismic waves from explosions. Two classes of problems have been addressed; nonlinear and linear near-source interactions. In both classes of problems displacements and tractions are saved on a closed surface in the linear region and the representation theorem is used to propagate the seismic waves to the far-field

Taming waveform inversion non-linearity through phase unwrapping of the model and objective functions

KAUST Repository

Alkhalifah, Tariq Ali

2012-09-25

Traveltime inversion focuses on the geometrical features of the waveform (traveltimes), which is generally smooth, and thus, tends to provide averaged (smoothed) information of the model. On other hand, general waveform inversion uses additional elements of the wavefield including amplitudes to extract higher resolution information, but this comes at the cost of introducing non-linearity to the inversion operator, complicating the convergence process. We use unwrapped phase-based objective functions in waveform inversion as a link between the two general types of inversions in a domain in which such contributions to the inversion process can be easily identified and controlled. The instantaneous traveltime is a measure of the average traveltime of the energy in a trace as a function of frequency. It unwraps the phase of wavefields yielding far less non-linearity in the objective function than that experienced with conventional wavefields, yet it still holds most of the critical wavefield information in its frequency dependency. However, it suffers from non-linearity introduced by the model (or reflectivity), as reflections from independent events in our model interact with each other. Unwrapping the phase of such a model can mitigate this non-linearity as well. Specifically, a simple modification to the inverted domain (or model), can reduce the effect of the model-induced non-linearity and, thus, make the inversion more convergent. Simple numerical examples demonstrate these assertions.

Taming waveform inversion non-linearity through phase unwrapping of the model and objective functions

KAUST Repository

Alkhalifah, Tariq Ali; Choi, Yun Seok

2012-01-01

Traveltime inversion focuses on the geometrical features of the waveform (traveltimes), which is generally smooth, and thus, tends to provide averaged (smoothed) information of the model. On other hand, general waveform inversion uses additional elements of the wavefield including amplitudes to extract higher resolution information, but this comes at the cost of introducing non-linearity to the inversion operator, complicating the convergence process. We use unwrapped phase-based objective functions in waveform inversion as a link between the two general types of inversions in a domain in which such contributions to the inversion process can be easily identified and controlled. The instantaneous traveltime is a measure of the average traveltime of the energy in a trace as a function of frequency. It unwraps the phase of wavefields yielding far less non-linearity in the objective function than that experienced with conventional wavefields, yet it still holds most of the critical wavefield information in its frequency dependency. However, it suffers from non-linearity introduced by the model (or reflectivity), as reflections from independent events in our model interact with each other. Unwrapping the phase of such a model can mitigate this non-linearity as well. Specifically, a simple modification to the inverted domain (or model), can reduce the effect of the model-induced non-linearity and, thus, make the inversion more convergent. Simple numerical examples demonstrate these assertions.

A numerical model for simulating electroosmotic micro- and nanochannel flows under non-Boltzmann equilibrium

Energy Technology Data Exchange (ETDEWEB)

Kim, Kyoungjin; Kwak, Ho Sang [School of Mechanical Engineering, Kumoh National Institute of Technology, 1 Yangho, Gumi, Gyeongbuk 730-701 (Korea, Republic of); Song, Tae-Ho, E-mail: kimkj@kumoh.ac.kr, E-mail: hskwak@kumoh.ac.kr, E-mail: thsong@kaist.ac.kr [Department of Mechanical, Aerospace and Systems Engineering, Korea Advanced Institute of Science and Technology, 373-1 Guseong, Yuseong, Daejeon 305-701 (Korea, Republic of)

2011-08-15

This paper describes a numerical model for simulating electroosmotic flows (EOFs) under non-Boltzmann equilibrium in a micro- and nanochannel. The transport of ionic species is represented by employing the Nernst-Planck equation. Modeling issues related to numerical difficulties are discussed, which include the handling of boundary conditions based on surface charge density, the associated treatment of electric potential and the evasion of nonlinearity due to the electric body force. The EOF in the entrance region of a straight channel is examined. The numerical results show that the present model is useful for the prediction of the EOFs requiring a fine resolution of the electric double layer under either the Boltzmann equilibrium or non-equilibrium. Based on the numerical results, the correlation between the surface charge density and the zeta potential is investigated.

Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model

Energy Technology Data Exchange (ETDEWEB)

Litvinenko, Yuri E. [Department of Mathematics, University of Waikato, P. B. 3105, Hamilton 3240 (New Zealand); Fichtner, Horst; Walter, Dominik [Institut für Theoretische Physik IV, Ruhr-Universität Bochum, Universitätsstrasse 150, D-44780 Bochum (Germany)

2017-05-20

We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatment of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.

NUMERICAL MODELLING AS NON-DESTRUCTIVE METHOD FOR THE ANALYSES AND DIAGNOSIS OF STONE STRUCTURES: MODELS AND POSSIBILITIES

Directory of Open Access Journals (Sweden)

Nataša Štambuk-Cvitanović

1999-12-01

Full Text Available Assuming the necessity of analysis, diagnosis and preservation of existing valuable stone masonry structures and ancient monuments in today European urban cores, numerical modelling become an efficient tool for the structural behaviour investigation. It should be supported by experimentally found input data and taken as a part of general combined approach, particularly non-destructive techniques on the structure/model within it. For the structures or their detail which may require more complex analyses three numerical models based upon finite elements technique are suggested: (1 standard linear model; (2 linear model with contact (interface elements; and (3 non-linear elasto-plastic and orthotropic model. The applicability of these models depend upon the accuracy of the approach or type of the problem, and will be presented on some characteristic samples.

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.)

A numerical analysis for non-linear radiation in MHD flow around a cylindrical surface with chemically reactive species

Directory of Open Access Journals (Sweden)

Junaid Ahmad Khan

2018-03-01

Full Text Available Boundary layer flow around a stretchable rough cylinder is modeled by taking into account boundary slip and transverse magnetic field effects. The main concern is to resolve heat/mass transfer problem considering non-linear radiative heat transfer and temperature/concentration jump aspects. Using conventional similarity approach, the equations of motion and heat transfer are converted into a boundary value problem whose solution is computed by shooting method for broad range of slip coefficients. The proposed numerical scheme appears to improve as the strengths of magnetic field and slip coefficients are enhanced. Axial velocity and temperature are considerably influenced by a parameter M which is inversely proportional to the radius of cylinder. A significant change in temperature profile is depicted for growing wall to ambient temperature ratio. Relevant physical quantities such as wall shear stress, local Nusselt number and local Sherwood number are elucidated in detail. Keywords: Stretchable boundary, Thermal radiation, Chemical reaction, Mathematical modeling, Non-linear differential system, Mass transfer

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...

Dynamics modeling for a rigid-flexible coupling system with nonlinear deformation field

International Nuclear Information System (INIS)

Deng Fengyan; He Xingsuo; Li Liang; Zhang Juan

2007-01-01

In this paper, a moving flexible beam, which incorporates the effect of the geometrically nonlinear kinematics of deformation, is investigated. Considering the second-order coupling terms of deformation in the longitudinal and transverse deflections, the exact nonlinear strain-displacement relations for a beam element are described. The shear strains formulated by the present modeling method in this paper are zero, so it is reasonable to use geometrically nonlinear deformation fields to demonstrate and simplify a flexible beam undergoing large overall motions. Then, considering the coupling terms of deformation in two dimensions, finite element shape functions of a beam element and Lagrange's equations are employed for deriving the coupling dynamical formulations. The complete expression of the stiffness matrix and all coupling terms are included in the formulations. A model consisting of a rotating planar flexible beam is presented. Then the frequency and dynamical response are studied, and the differences among the zero-order model, first-order coupling model and the new present model are discussed. Numerical examples demonstrate that a 'stiffening beam' can be obtained, when more coupling terms of deformation are added to the longitudinal and transverse deformation field. It is shown that the traditional zero-order and first-order coupling models may not provide an exact dynamic model in some cases

Numerical approximations of nonlinear fractional differential difference equations by using modified He-Laplace method

Directory of Open Access Journals (Sweden)

J. Prakash

2016-03-01

Full Text Available In this paper, a numerical algorithm based on a modified He-Laplace method (MHLM is proposed to solve space and time nonlinear fractional differential-difference equations (NFDDEs arising in physical phenomena such as wave phenomena in fluids, coupled nonlinear optical waveguides and nanotechnology fields. The modified He-Laplace method is a combined form of the fractional homotopy perturbation method and Laplace transforms method. The nonlinear terms can be easily decomposed by the use of He’s polynomials. This algorithm has been tested against time-fractional differential-difference equations such as the modified Lotka Volterra and discrete (modified KdV equations. The proposed scheme grants the solution in the form of a rapidly convergent series. Three examples have been employed to illustrate the preciseness and effectiveness of the proposed method. The achieved results expose that the MHLM is very accurate, efficient, simple and can be applied to other nonlinear FDDEs.

The inherent complexity in nonlinear business cycle model in resonance

International Nuclear Information System (INIS)

Ma Junhai; Sun Tao; Liu Lixia

2008-01-01

Based on Abraham C.-L. Chian's research, we applied nonlinear dynamic system theory to study the first-order and second-order approximate solutions to one category of the nonlinear business cycle model in resonance condition. We have also analyzed the relation between amplitude and phase of second-order approximate solutions as well as the relation between outer excitements' amplitude, frequency approximate solutions, and system bifurcation parameters. Then we studied the system quasi-periodical solutions, annulus periodical solutions and the path leading to system bifurcation and chaotic state with different parameter combinations. Finally, we conducted some numerical simulations for various complicated circumstances. Therefore this research will lay solid foundation for detecting the complexity of business cycles and systems in the future

An ensemble Kalman filter for statistical estimation of physics constrained nonlinear regression models

International Nuclear Information System (INIS)

Harlim, John; Mahdi, Adam; Majda, Andrew J.

2014-01-01

A central issue in contemporary science is the development of nonlinear data driven statistical–dynamical models for time series of noisy partial observations from nature or a complex model. It has been established recently that ad-hoc quadratic multi-level regression models can have finite-time blow-up of statistical solutions and/or pathological behavior of their invariant measure. Recently, a new class of physics constrained nonlinear regression models were developed to ameliorate this pathological behavior. Here a new finite ensemble Kalman filtering algorithm is developed for estimating the state, the linear and nonlinear model coefficients, the model and the observation noise covariances from available partial noisy observations of the state. Several stringent tests and applications of the method are developed here. In the most complex application, the perfect model has 57 degrees of freedom involving a zonal (east–west) jet, two topographic Rossby waves, and 54 nonlinearly interacting Rossby waves; the perfect model has significant non-Gaussian statistics in the zonal jet with blocked and unblocked regimes and a non-Gaussian skewed distribution due to interaction with the other 56 modes. We only observe the zonal jet contaminated by noise and apply the ensemble filter algorithm for estimation. Numerically, we find that a three dimensional nonlinear stochastic model with one level of memory mimics the statistical effect of the other 56 modes on the zonal jet in an accurate fashion, including the skew non-Gaussian distribution and autocorrelation decay. On the other hand, a similar stochastic model with zero memory levels fails to capture the crucial non-Gaussian behavior of the zonal jet from the perfect 57-mode model

Extinction and persistence of a stochastic nonlinear SIS epidemic model with jumps

Science.gov (United States)

Ge, Qing; Ji, Guilin; Xu, Jiabo; Fan, Xiaolin

2016-11-01

In this paper, Brownian motion and L e ´ vy jumps are introduced to a SIS type epidemic model with nonlinear incidence rate. The dynamical behavior of the considered model is investigated. In order to reveal the extinction and permanence of the disease, two threshold values R˜0 ,R¯0 are showed. We find that if R˜0 1, the disease may be persistent. Finally, the numerical simulations are presented to illustrate our mathematical results.

A numerical approach for assessing effects of shear on equivalent permeability and nonlinear flow characteristics of 2-D fracture networks

Science.gov (United States)

Liu, Richeng; Li, Bo; Jiang, Yujing; Yu, Liyuan

2018-01-01

Hydro-mechanical properties of rock fractures are core issues for many geoscience and geo-engineering practices. Previous experimental and numerical studies have revealed that shear processes could greatly enhance the permeability of single rock fractures, yet the shear effects on hydraulic properties of fractured rock masses have received little attention. In most previous fracture network models, single fractures are typically presumed to be formed by parallel plates and flow is presumed to obey the cubic law. However, related studies have suggested that the parallel plate model cannot realistically represent the surface characters of natural rock fractures, and the relationship between flow rate and pressure drop will no longer be linear at sufficiently large Reynolds numbers. In the present study, a numerical approach was established to assess the effects of shear on the hydraulic properties of 2-D discrete fracture networks (DFNs) in both linear and nonlinear regimes. DFNs considering fracture surface roughness and variation of aperture in space were generated using an originally developed code DFNGEN. Numerical simulations by solving Navier-Stokes equations were performed to simulate the fluid flow through these DFNs. A fracture that cuts through each model was sheared and by varying the shear and normal displacements, effects of shear on equivalent permeability and nonlinear flow characteristics of DFNs were estimated. The results show that the critical condition of quantifying the transition from a linear flow regime to a nonlinear flow regime is: 10-4 〈 J hydraulic gradient. When the fluid flow is in a linear regime (i.e., J reduce the equivalent permeability significantly in the orientation perpendicular to the sheared fracture as much as 53.86% when J = 1, shear displacement Ds = 7 mm, and normal displacement Dn = 1 mm. By fitting the calculated results, the mathematical expression for δ2 is established to help choose proper governing equations when

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay.

Science.gov (United States)

Akimenko, Vitalii; Anguelov, Roumen

2017-12-01

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.

Thin layer model for nonlinear evolution of the Rayleigh-Taylor instability

Science.gov (United States)

Zhao, K. G.; Wang, L. F.; Xue, C.; Ye, W. H.; Wu, J. F.; Ding, Y. K.; Zhang, W. Y.

2018-03-01

On the basis of the thin layer approximation [Ott, Phys. Rev. Lett. 29, 1429 (1972)], a revised thin layer model for incompressible Rayleigh-Taylor instability has been developed to describe the deformation and nonlinear evolution of the perturbed interface. The differential equations for motion are obtained by analyzing the forces (the gravity and pressure difference) of fluid elements (i.e., Newton's second law). The positions of the perturbed interface are obtained from the numerical solution of the motion equations. For the case of vacuum on both sides of the layer, the positions of the upper and lower interfaces obtained from the revised thin layer approximation agree with that from the weakly nonlinear (WN) model of a finite-thickness fluid layer [Wang et al., Phys. Plasmas 21, 122710 (2014)]. For the case considering the fluids on both sides of the layer, the bubble-spike amplitude from the revised thin layer model agrees with that from the WN model [Wang et al., Phys. Plasmas 17, 052305 (2010)] and the expanded Layzer's theory [Goncharov, Phys. Rev. Lett. 88, 134502 (2002)] in the early nonlinear growth regime. Note that the revised thin layer model can be applied to investigate the perturbation growth at arbitrary Atwood numbers. In addition, the large deformation (the large perturbed amplitude and the arbitrary perturbed distributions) in the initial stage can also be described by the present model.

Finite element model for nonlinear shells of revolution

International Nuclear Information System (INIS)

Cook, W.A.

1979-01-01

Nuclear material shipping containers have shells of revolution as basic structural components. Analytically modeling the response of these containers to severe accident impact conditions requires a nonlinear shell-of-revolution model that accounts for both geometric and material nonlinearities. Existing models are limited to large displacements, small rotations, and nonlinear materials. The paper presents a finite element model for a nonlinear shell of revolution that will account for large displacements, large strains, large rotations, and nonlinear materials

Attitude Control of a Single Tilt Tri-Rotor UAV System: Dynamic Modeling and Each Channel's Nonlinear Controllers Design

Directory of Open Access Journals (Sweden)

Juing-Shian Chiou

2013-01-01

Full Text Available This paper has implemented nonlinear control strategy for the single tilt tri-rotor aerial robot. Based on Newton-Euler’s laws, the linear and nonlinear mathematical models of tri-rotor UAVs are obtained. A numerical analysis using Newton-Raphson method is chosen for finding hovering equilibrium point. Back-stepping nonlinear controller design is based on constructing Lyapunov candidate function for closed-loop system. By imitating the linguistic logic of human thought, fuzzy logic controllers (FLCs are designed based on control rules and membership functions, which are much less rigid than the calculations computers generally perform. Effectiveness of the controllers design scheme is shown through nonlinear simulation model on each channel.

Fuzzy Control Model and Simulation for Nonlinear Supply Chain System with Lead Times

Directory of Open Access Journals (Sweden)

Songtao Zhang

2017-01-01

Full Text Available A new fuzzy robust control strategy for the nonlinear supply chain system in the presence of lead times is proposed. Based on Takagi-Sugeno fuzzy control system, the fuzzy control model of the nonlinear supply chain system with lead times is constructed. Additionally, we design a fuzzy robust H∞ control strategy taking the definition of maximal overlapped-rules group into consideration to restrain the impacts such as those caused by lead times, switching actions among submodels, and customers’ stochastic demands. This control strategy can not only guarantee that the nonlinear supply chain system is robustly asymptotically stable but also realize soft switching among subsystems of the nonlinear supply chain to make the less fluctuation of the system variables by introducing the membership function of fuzzy system. The comparisons between the proposed fuzzy robust H∞ control strategy and the robust H∞ control strategy are finally illustrated through numerical simulations on a two-stage nonlinear supply chain with lead times.

Sine sweep and steady-state response of simplified solar array models with nonlinear elements

NARCIS (Netherlands)

Fey, R.H.B.; van Liempt, F.P.H.

2002-01-01

In this paper the nonlinear dynamic behaviour of two simplified solar array systems is investigated experimentally and numerically. A simplified beam model supported by one snubber (a bilinear spring which can only take compressive forces) is used to investigate the dynamics of the extension arm on

A general 3-D nonlinear magnetostrictive constitutive model for soft ferromagnetic materials

International Nuclear Information System (INIS)

Zhou Haomiao; Zhou Youhe; Zheng Xiaojing; Ye Qiang; Wei Jing

2009-01-01

In this paper, a new general nonlinear magnetostrictive constitutive model is proposed for soft ferromagnetic materials, and it can predict magnetostrictive strain and magnetization curves under various pre-stresses. From the viewpoint of magnetic domain, it is based on the important physical fact that a nonlinear part of the elastic strain produced by magnetic domain wall motion under a pre-stress is responsible for the change of the maximum magnetostrictive strain in accordance with the pre-stress. Then the reduction of magnetostrictive strain from the maximum is caused by the domain rotation. Meanwhile, the magnetization under various pre-stresses in this model is introduced by magnetostrictive effect under the same pre-stress. A simplified 3-D model is put forward by means of linearizing the nonlinear function, i.e. the nonlinear part of the elastic strain produced by domain wall motion, and by using the quartic of magnetization to describe domain rotation. Besides, for the convenience of engineering applications, two-dimensional (plate or film) and one-dimensional (rod) models are also given for isotropic materials and their application ranges are discussed too. In comparison with the experimental data of Kuruzar and Jiles, it is found that this model can predict magnetostrictive strain and magnetization curves under various pre-stresses. The numerical simulation further illustrates that the new model can effectively describe the effects of the pre-stress or residual stress on the magnetization and magnetostrictive strain curves. Additionally, this model can be degenerated to the existing magnetostrictive constitutive model for giant magnetostrictive materials (GMM), i.e. a special soft ferromagnetic material

A Nonlinear Diffusion Equation-Based Model for Ultrasound Speckle Noise Removal

Science.gov (United States)

Zhou, Zhenyu; Guo, Zhichang; Zhang, Dazhi; Wu, Boying

2018-04-01

Ultrasound images are contaminated by speckle noise, which brings difficulties in further image analysis and clinical diagnosis. In this paper, we address this problem in the view of nonlinear diffusion equation theories. We develop a nonlinear diffusion equation-based model by taking into account not only the gradient information of the image, but also the information of the gray levels of the image. By utilizing the region indicator as the variable exponent, we can adaptively control the diffusion type which alternates between the Perona-Malik diffusion and the Charbonnier diffusion according to the image gray levels. Furthermore, we analyze the proposed model with respect to the theoretical and numerical properties. Experiments show that the proposed method achieves much better speckle suppression and edge preservation when compared with the traditional despeckling methods, especially in the low gray level and low-contrast regions.

Dynamics in a nonlinear Keynesian good market model

International Nuclear Information System (INIS)

Naimzada, Ahmad; Pireddu, Marina

2014-01-01

In this paper, we show how a rich variety of dynamical behaviors can emerge in the standard Keynesian income-expenditure model when a nonlinearity is introduced, both in the cases with and without endogenous government spending. A specific sigmoidal functional form is used for the adjustment mechanism of income with respect to the excess demand, in order to bound the income variation. With the aid of analytical and numerical tools, we investigate the stability conditions, bifurcations, as well as periodic and chaotic dynamics. Globally, we study multistability phenomena, i.e., the coexistence of different kinds of attractors

Oscillating nonlinear acoustic shock waves

DEFF Research Database (Denmark)

Gaididei, Yuri; Rasmussen, Anders Rønne; Christiansen, Peter Leth

2016-01-01

We investigate oscillating shock waves in a tube using a higher order weakly nonlinear acoustic model. The model includes thermoviscous effects and is non isentropic. The oscillating shock waves are generated at one end of the tube by a sinusoidal driver. Numerical simulations show that at resona......We investigate oscillating shock waves in a tube using a higher order weakly nonlinear acoustic model. The model includes thermoviscous effects and is non isentropic. The oscillating shock waves are generated at one end of the tube by a sinusoidal driver. Numerical simulations show...... polynomial in the space and time variables, we find analytical approximations to the observed single shock waves in an infinitely long tube. Using perturbation theory for the driven acoustic system approximative analytical solutions for the off resonant case are determined....

A general U-block model-based design procedure for nonlinear polynomial control systems

Science.gov (United States)

Zhu, Q. M.; Zhao, D. Y.; Zhang, Jianhua

2016-10-01

The proposition of U-model concept (in terms of 'providing concise and applicable solutions for complex problems') and a corresponding basic U-control design algorithm was originated in the first author's PhD thesis. The term of U-model appeared (not rigorously defined) for the first time in the first author's other journal paper, which established a framework for using linear polynomial control system design approaches to design nonlinear polynomial control systems (in brief, linear polynomial approaches → nonlinear polynomial plants). This paper represents the next milestone work - using linear state-space approaches to design nonlinear polynomial control systems (in brief, linear state-space approaches → nonlinear polynomial plants). The overall aim of the study is to establish a framework, defined as the U-block model, which provides a generic prototype for using linear state-space-based approaches to design the control systems with smooth nonlinear plants/processes described by polynomial models. For analysing the feasibility and effectiveness, sliding mode control design approach is selected as an exemplary case study. Numerical simulation studies provide a user-friendly step-by-step procedure for the readers/users with interest in their ad hoc applications. In formality, this is the first paper to present the U-model-oriented control system design in a formal way and to study the associated properties and theorems. The previous publications, in the main, have been algorithm-based studies and simulation demonstrations. In some sense, this paper can be treated as a landmark for the U-model-based research from intuitive/heuristic stage to rigour/formal/comprehensive studies.

A numerical study of one-dimensional replicating patterns in reaction-diffusion systems with non-linear diffusion coefficients

International Nuclear Information System (INIS)

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

1998-01-01

A numerical study of the dynamics of pattern evolution in reaction-diffusion systems is performed, although limited to one spatial dimension. The diffusion coefficients are nonlinear, based on powers of the scalar variables. The system keeps the dynamics of previous studies in the literature, but the presence of nonlinear diffusion generates a field of strong nonlinear interactions due to the presence of receding travelling waves. This field is limited by the plane of symmetry of the space domain and the last born outgoing travelling wave. These effects are discussed. (author). 10 refs., 7 figs

Superspace formulation of new nonlinear sigma models

International Nuclear Information System (INIS)

Gates, S.J. Jr.

1983-07-01

The superspace formulation of two classes of supersymmetric nonlinear σ-models are presented. Two alternative N=1 superspace formulations are given for the d=2 supersymmetric nonlinear σ-models with Killing vector potentials: (a) formulation uses an active central charge and, (b) formulation uses a spurion superfield without inducing a classical breakdown of supersymmetry. The N=2 vector multiplet is used to construct a new class of d=4 nonlinear σ-models which when reduced to d=2 possess N=4 supersymmetry. Implications of these two classes of nonlinear σ-models for N>=4 superfield supergravity are discussed. (author)

An axisymmetrical non-linear finite element model for induction heating in injection molding tools

DEFF Research Database (Denmark)

Guerrier, Patrick; Nielsen, Kaspar Kirstein; Menotti, Stefano

2016-01-01

To analyze the heating and cooling phase of an induction heated injection molding tool accurately, the temperature dependent magnetic properties, namely the non-linear B-H curves, need to be accounted for in an induction heating simulation. Hence, a finite element model has been developed......, including the non-linear temperature dependent magnetic data described by a three-parameter modified Frohlich equation fitted to the magnetic saturation curve, and solved with an iterative procedure. The numerical calculations are compared with experiments conducted with two types of induction coils, built...... in to the injection molding tool. The model shows very good agreement with the experimental temperature measurements. It is also shown that the non-linearity can be used without the temperature dependency in some cases, and a proposed method is presented of how to estimate an effective linear permeability to use...

3D nonlinear numerical simulation of the current-convective instability in detached diverter plasma

Science.gov (United States)

Stepanenko, Alexander; Krasheninnikov, Sergei

2017-10-01

One of the possible mechanisms responsible for strong radiation fluctuations observed in the recent experiments with detached plasmas at ASDEX Upgrade [Potzel et al., Nuclear Fusion, 2014] can be related to the onset of the current-convective instability (CCI) driven by strong asymmetry of detachment in the inner and outer tokamak divertors [Krasheninnikov and Smolyakov, PoP, 2016]. In this study we present the first results of 3D nonlinear numerical simulations of the CCI in divertor plasma for the conditions relevant to the AUG experiment. The general physical model used to simulate the CCI, qualitative estimates for the instability characteristic growth rate and transverse wavelengths derived for plasma, which is spatially inhomogeneous both across and along the magnetic field lines, are presented. The simulation results, demonstrating nonlinear dynamics of the CCI, provide the frequency spectra of turbulent divertor plasma fluctuations showing good agreement with the available experimental data. This material is based upon the work supported by the U.S. Department of Energy under Award No. DE-FG02-04ER54739 at UCSD and by the Russian Ministry of Education and Science Grant No. 14.Y26.31.0008 at MEPhI.

Hopf bifurcation in love dynamical models with nonlinear couples and time delays

International Nuclear Information System (INIS)

Liao Xiaofeng; Ran Jiouhong

2007-01-01

A love dynamical models with nonlinear couples and two delays is considered. Local stability of this model is studied by analyzing the associated characteristic transcendental equation. We find that the Hopf bifurcation occurs when the sum of the two delays varies and passes a sequence of critical values. The stability and direction of the Hopf bifurcation are determined by applying the normal form theory and the center manifold theorem. Numerical example is given to illustrate our results

Numerical study of propagation properties of surface plasmon polaritons in nonlinear media

KAUST Repository

Sagor, Rakibul Hasan

2016-03-29

We present a time-domain algorithm for simulating nonlinear propagation of surface plasmon polaritons (SPPs) in chalcogenide glass. Due to the high non-linearity property and strong dispersion and confinement chalcogenide glasses are widely known as ultrafast nonlinear materials. We have used the finite difference time domain (FDTD) method to develop the simulation algorithm for the current analysis. We have modeled the frequency dependent dispersion properties and third order nonlinearity property of chalcogenide glass utilizing the general polarization algorithm merged in the auxiliary differential equation (ADE) method. The propagation dynamics of the whole structure with and without third order nonlinearity property of chalcogenide glass have been simulated and the effect of nonlinearity on the propagation properties of SPP has been investigated. © 2016 EDP Sciences, SIF, Springer-Verlag Berlin Heidelberg.

A Two-Step Hybrid Approach for Modeling the Nonlinear Dynamic Response of Piezoelectric Energy Harvesters

Directory of Open Access Journals (Sweden)

Claudio Maruccio

2018-01-01

Full Text Available An effective hybrid computational framework is described here in order to assess the nonlinear dynamic response of piezoelectric energy harvesting devices. The proposed strategy basically consists of two steps. First, fully coupled multiphysics finite element (FE analyses are performed to evaluate the nonlinear static response of the device. An enhanced reduced-order model is then derived, where the global dynamic response is formulated in the state-space using lumped coefficients enriched with the information derived from the FE simulations. The electromechanical response of piezoelectric beams under forced vibrations is studied by means of the proposed approach, which is also validated by comparing numerical predictions with some experimental results. Such numerical and experimental investigations have been carried out with the main aim of studying the influence of material and geometrical parameters on the global nonlinear response. The advantage of the presented approach is that the overall computational and experimental efforts are significantly reduced while preserving a satisfactory accuracy in the assessment of the global behavior.

Rigid-body-spring model numerical analysis of joint performance of engineered cementitious composites and concrete

Science.gov (United States)

Khmurovska, Y.; Štemberk, P.; Křístek, V.

2017-09-01

This paper presents a numerical investigation of effectiveness of using engineered cementitious composites with polyvinyl alcohol fibers for concrete cover layer repair. A numerical model of a monolithic concaved L-shaped concrete structural detail which is strengthened with an engineered cementitious composite layer with polyvinyl alcohol fibers is created and loaded with bending moment. The numerical analysis employs nonlinear 3-D Rigid-Body-Spring Model. The proposed material model shows reliable results and can be used in further studies. The engineered cementitious composite shows extremely good performance in tension due to the strain-hardening effect. Since durability of the bond can be decreased significantly by its degradation due to the thermal loading, this effect should be also taken into account in the future work, as well as the experimental investigation, which should be performed for validation of the proposed numerical model.

Prediction of Mind-Wandering with Electroencephalogram and Non-linear Regression Modeling.

Science.gov (United States)

Kawashima, Issaku; Kumano, Hiroaki

2017-01-01

Mind-wandering (MW), task-unrelated thought, has been examined by researchers in an increasing number of articles using models to predict whether subjects are in MW, using numerous physiological variables. However, these models are not applicable in general situations. Moreover, they output only binary classification. The current study suggests that the combination of electroencephalogram (EEG) variables and non-linear regression modeling can be a good indicator of MW intensity. We recorded EEGs of 50 subjects during the performance of a Sustained Attention to Response Task, including a thought sampling probe that inquired the focus of attention. We calculated the power and coherence value and prepared 35 patterns of variable combinations and applied Support Vector machine Regression (SVR) to them. Finally, we chose four SVR models: two of them non-linear models and the others linear models; two of the four models are composed of a limited number of electrodes to satisfy model usefulness. Examination using the held-out data indicated that all models had robust predictive precision and provided significantly better estimations than a linear regression model using single electrode EEG variables. Furthermore, in limited electrode condition, non-linear SVR model showed significantly better precision than linear SVR model. The method proposed in this study helps investigations into MW in various little-examined situations. Further, by measuring MW with a high temporal resolution EEG, unclear aspects of MW, such as time series variation, are expected to be revealed. Furthermore, our suggestion that a few electrodes can also predict MW contributes to the development of neuro-feedback studies.

Prediction of Mind-Wandering with Electroencephalogram and Non-linear Regression Modeling

Directory of Open Access Journals (Sweden)

Issaku Kawashima

2017-07-01

Full Text Available Mind-wandering (MW, task-unrelated thought, has been examined by researchers in an increasing number of articles using models to predict whether subjects are in MW, using numerous physiological variables. However, these models are not applicable in general situations. Moreover, they output only binary classification. The current study suggests that the combination of electroencephalogram (EEG variables and non-linear regression modeling can be a good indicator of MW intensity. We recorded EEGs of 50 subjects during the performance of a Sustained Attention to Response Task, including a thought sampling probe that inquired the focus of attention. We calculated the power and coherence value and prepared 35 patterns of variable combinations and applied Support Vector machine Regression (SVR to them. Finally, we chose four SVR models: two of them non-linear models and the others linear models; two of the four models are composed of a limited number of electrodes to satisfy model usefulness. Examination using the held-out data indicated that all models had robust predictive precision and provided significantly better estimations than a linear regression model using single electrode EEG variables. Furthermore, in limited electrode condition, non-linear SVR model showed significantly better precision than linear SVR model. The method proposed in this study helps investigations into MW in various little-examined situations. Further, by measuring MW with a high temporal resolution EEG, unclear aspects of MW, such as time series variation, are expected to be revealed. Furthermore, our suggestion that a few electrodes can also predict MW contributes to the development of neuro-feedback studies.

Numerical Investigation on Vortex-Structure Interaction Generating Aerodynamic Noises for Rod-Airfoil Models

Directory of Open Access Journals (Sweden)

FeiFei Liu

2017-01-01

Full Text Available In past several decades, vortex-structure interaction generated aerodynamic noise became one of the main concerns in aircraft design. In order to understand the mechanism, the acoustic analogy method combined with the RANS-based nonlinear acoustics solver (NLAS is investigated. The numerical method is firstly evaluated by the experiment data of the classic rod-airfoil model. Compared with the traditional analogy methods, the RANS/NLAS can capture the nonlinear aerodynamic noise more accurately with lower gird requirements. Then different rod-airfoil configurations were simulated to investigate the aeroacoustic interaction effects. The numerical results are in good agreement with those of the earlier experimental research. It is found that the vortex-shedding crash to the airfoil is the main reason for the noise generation which is dependent on the configurations, distance, and flow conditions.

Numerical Simulation of Entropy Growth for a Nonlinear Evolutionary Model of Random Markets

Directory of Open Access Journals (Sweden)

Mahdi Keshtkar

2016-01-01

Full Text Available In this communication, the generalized continuous economic model for random markets is revisited. In this model for random markets, agents trade by pairs and exchange their money in a random and conservative way. They display the exponential wealth distribution as asymptotic equilibrium, independently of the effectiveness of the transactions and of the limitation of the total wealth. In the current work, entropy of mentioned model is defined and then some theorems on entropy growth of this evolutionary problem are given. Furthermore, the entropy increasing by simulation on some numerical examples is verified.

Numerical Calculation and Measurement of Nonlinear Acoustic Fields in Ultrasound Diagnosis

Science.gov (United States)

Kawagishi, Tetsuya; Saito, Shigemi; Mine, Yoshitaka

2002-05-01

In order to develop a tool for designing on the ultrasonic probe and its peripheral devices for tissue-harmonic-imaging systems, a study is carried out to compare the calculation and observation results of nonlinear acoustic fields for a diagnostic ultrasound system. The pulsed ultrasound with a center frequency of 2.5 MHz is emanated from a weakly focusing sector probe with a 6.5 mm aperture radius and a 50 mm focal length into an agar phantom with an attenuation coefficient of about 0.6 dB/cm/MHz or 1.2 dB/cm/MHz. The nonlinear acoustic field is measured using a needle-type hydrophone. The calculation is based on the Khokhlov-Zabolotskaya-Kuznetsov(KZK) equation which is modified so that the frequency dependence of the attenuation coefficient is the same as that in biological tissue. This equation is numerically solved with the implicit backward method employing the iterative method. The measured and calculated amplitude spectra show good agreement with each other.

Forecasting with nonlinear time series models

DEFF Research Database (Denmark)

Kock, Anders Bredahl; Teräsvirta, Timo

In this paper, nonlinear models are restricted to mean nonlinear parametric models. Several such models popular in time series econo- metrics are presented and some of their properties discussed. This in- cludes two models based on universal approximators: the Kolmogorov- Gabor polynomial model...... applied to economic fore- casting problems, is briefly highlighted. A number of large published studies comparing macroeconomic forecasts obtained using different time series models are discussed, and the paper also contains a small simulation study comparing recursive and direct forecasts in a partic...... and two versions of a simple artificial neural network model. Techniques for generating multi-period forecasts from nonlinear models recursively are considered, and the direct (non-recursive) method for this purpose is mentioned as well. Forecasting with com- plex dynamic systems, albeit less frequently...

Fuzzy model-based servo and model following control for nonlinear systems.

Science.gov (United States)

Ohtake, Hiroshi; Tanaka, Kazuo; Wang, Hua O

2009-12-01

This correspondence presents servo and nonlinear model following controls for a class of nonlinear systems using the Takagi-Sugeno fuzzy model-based control approach. First, the construction method of the augmented fuzzy system for continuous-time nonlinear systems is proposed by differentiating the original nonlinear system. Second, the dynamic fuzzy servo controller and the dynamic fuzzy model following controller, which can make outputs of the nonlinear system converge to target points and to outputs of the reference system, respectively, are introduced. Finally, the servo and model following controller design conditions are given in terms of linear matrix inequalities. Design examples illustrate the utility of this approach.

Nonlinear system identification based on Takagi-Sugeno fuzzy modeling and unscented Kalman filter.

Science.gov (United States)

Vafamand, Navid; Arefi, Mohammad Mehdi; Khayatian, Alireza

2018-03-01

This paper proposes two novel Kalman-based learning algorithms for an online Takagi-Sugeno (TS) fuzzy model identification. The proposed approaches are designed based on the unscented Kalman filter (UKF) and the concept of dual estimation. Contrary to the extended Kalman filter (EKF) which utilizes derivatives of nonlinear functions, the UKF employs the unscented transformation. Consequently, non-differentiable membership functions can be considered in the structure of the TS models. This makes the proposed algorithms to be applicable for the online parameter calculation of wider classes of TS models compared to the recently published papers concerning the same issue. Furthermore, because of the great capability of the UKF in handling severe nonlinear dynamics, the proposed approaches can effectively approximate the nonlinear systems. Finally, numerical and practical examples are provided to show the advantages of the proposed approaches. Simulation results reveal the effectiveness of the proposed methods and performance improvement based on the root mean square (RMS) of the estimation error compared to the existing results. Copyright © 2018 ISA. Published by Elsevier Ltd. All rights reserved.

Identification of nonlinear anelastic models

International Nuclear Information System (INIS)

Draganescu, G E; Bereteu, L; Ercuta, A

2008-01-01

A useful nonlinear identification technique applied to the anelastic and rheologic models is presented in this paper. First introduced by Feldman, the method is based on the Hilbert transform, and is currently used for identification of the nonlinear vibrations

Relation between nonlinear models and gauge ambiguities

International Nuclear Information System (INIS)

Balachandran, A.P.; Ramachandran, R.; Rupertsberger, H.; Skagerstam, B.S.

1980-01-01

We show that the solutions of a class of nonlinear models also generate gauge ambiguities in the vacuum sector of Yang-Mills theories. Our results extend known connections between gauge ambiguities and certain nonlinear sigma-models, and clarify the underlying group theory. For many nonlinear models, we also give a simple, intrinsic parametrization of physical fields (which have values in a homogeneous space of a group). (orig.)

Nonlinear Schrödinger approach to European option pricing

Science.gov (United States)

Wróblewski, Marcin

2017-05-01

This paper deals with numerical option pricing methods based on a Schrödinger model rather than the Black-Scholes model. Nonlinear Schrödinger boundary value problems seem to be alternatives to linear models which better reflect the complexity and behavior of real markets. Therefore, based on the nonlinear Schrödinger option pricing model proposed in the literature, in this paper a model augmented by external atomic potentials is proposed and numerically tested. In terms of statistical physics the developed model describes the option in analogy to a pair of two identical quantum particles occupying the same state. The proposed model is used to price European call options on a stock index. the model is calibrated using the Levenberg-Marquardt algorithm based on market data. A Runge-Kutta method is used to solve the discretized boundary value problem numerically. Numerical results are provided and discussed. It seems that our proposal more accurately models phenomena observed in the real market than do linear models.

Flute-like musical instruments: A toy model investigated through numerical continuation

Science.gov (United States)

Terrien, Soizic; Vergez, Christophe; Fabre, Benoît

2013-07-01

Self-sustained musical instruments (bowed string, woodwind and brass instruments) can be modelled by nonlinear lumped dynamical systems. Among these instruments, flutes and flue organ pipes present the particularity to be modelled as a delay dynamical system. In this paper, such a system, a toy model of flute-like instruments, is studied using numerical continuation. Equilibrium and periodic solutions are explored with respect to the blowing pressure, with focus on amplitude and frequency evolutions along the different solution branches, as well as "jumps" between periodic solution branches. The influence of a second model parameter (namely the inharmonicity) on the behaviour of the system is addressed. It is shown that harmonicity plays a key role in the presence of hysteresis or quasiperiodic regime. Throughout the paper, experimental results on a real instrument are presented to illustrate various phenomena, and allow some qualitative comparisons with numerical results.

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

A numerical algorithm for modelling boson-fermion stars in dilatonic gravity

CERN Document Server

Boyadzhiev, T L; Todorov, M D; Yazadjiev, S S

2002-01-01

We investigate numerically the class of models of the static spherically symmetric boson-fermion stars in the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of such stars is interpreted as a nonlinear two-parametric eigenvalue problem. The first of the parameters is the unknown internal boundary (the radius of the fermionic part of the star) R sub s , and the second one represents the frequency OMEGA of the time oscillations of the boson field. To solve this problem, the whole space [0, infinity) is splitted into two domains: internal [0, R sub s] (inside the star) and external [R sub s , infinity) (outside the star). In each domain the physical model leads to two nonlinear boundary value problems in respect to metric functions, the functions describing the fermionic and bosonic matter, and the dilaton field. These boundary value problems have different dimensions inside and outside the star, respectively. The solutions in these regions are obtained separately and ma...

A Nonlinear Model of Mix Coil Spring – Rubber for Vertical Suspension of Railway Vehicle

Directory of Open Access Journals (Sweden)

Dumitriu Mădălina

2016-03-01

Full Text Available The paper focuses on a nonlinear model to represent the mechanical behaviour of a mix coil spring - rubber used in the secondary suspension of passenger rail vehicles. The principle of the model relies on overlapping of the forces corresponding to three components - the elastic component, the viscous component and the dry friction component. The model has two sources on non-linearity, in the elastic force and the friction force, respectively. The main attributes of the model are made visible by its response to an imposed displacement-type harmonic excitation. The results thus obtained from the applications of numerical simulation show a series of basic properties of the model, namely the dependence on amplitude and the excitation frequency of the model response, as well as of its stiffness and damping.

Estimation methods for nonlinear state-space models in ecology

DEFF Research Database (Denmark)

Pedersen, Martin Wæver; Berg, Casper Willestofte; Thygesen, Uffe Høgsbro

2011-01-01

The use of nonlinear state-space models for analyzing ecological systems is increasing. A wide range of estimation methods for such models are available to ecologists, however it is not always clear, which is the appropriate method to choose. To this end, three approaches to estimation in the theta...... logistic model for population dynamics were benchmarked by Wang (2007). Similarly, we examine and compare the estimation performance of three alternative methods using simulated data. The first approach is to partition the state-space into a finite number of states and formulate the problem as a hidden...... Markov model (HMM). The second method uses the mixed effects modeling and fast numerical integration framework of the AD Model Builder (ADMB) open-source software. The third alternative is to use the popular Bayesian framework of BUGS. The study showed that state and parameter estimation performance...

Physics constrained nonlinear regression models for time series

International Nuclear Information System (INIS)

Majda, Andrew J; Harlim, John

2013-01-01

A central issue in contemporary science is the development of data driven statistical nonlinear dynamical models for time series of partial observations of nature or a complex physical model. It has been established recently that ad hoc quadratic multi-level regression (MLR) models can have finite-time blow up of statistical solutions and/or pathological behaviour of their invariant measure. Here a new class of physics constrained multi-level quadratic regression models are introduced, analysed and applied to build reduced stochastic models from data of nonlinear systems. These models have the advantages of incorporating memory effects in time as well as the nonlinear noise from energy conserving nonlinear interactions. The mathematical guidelines for the performance and behaviour of these physics constrained MLR models as well as filtering algorithms for their implementation are developed here. Data driven applications of these new multi-level nonlinear regression models are developed for test models involving a nonlinear oscillator with memory effects and the difficult test case of the truncated Burgers–Hopf model. These new physics constrained quadratic MLR models are proposed here as process models for Bayesian estimation through Markov chain Monte Carlo algorithms of low frequency behaviour in complex physical data. (paper)

Non-linear finite element modeling

DEFF Research Database (Denmark)

Mikkelsen, Lars Pilgaard

The note is written for courses in "Non-linear finite element method". The note has been used by the author teaching non-linear finite element modeling at Civil Engineering at Aalborg University, Computational Mechanics at Aalborg University Esbjerg, Structural Engineering at the University...

Correlations and Non-Linear Probability Models

DEFF Research Database (Denmark)

Breen, Richard; Holm, Anders; Karlson, Kristian Bernt

2014-01-01

the dependent variable of the latent variable model and its predictor variables. We show how this correlation can be derived from the parameters of non-linear probability models, develop tests for the statistical significance of the derived correlation, and illustrate its usefulness in two applications. Under......Although the parameters of logit and probit and other non-linear probability models are often explained and interpreted in relation to the regression coefficients of an underlying linear latent variable model, we argue that they may also be usefully interpreted in terms of the correlations between...... certain circumstances, which we explain, the derived correlation provides a way of overcoming the problems inherent in cross-sample comparisons of the parameters of non-linear probability models....

Non-linear Loudspeaker Unit Modelling

DEFF Research Database (Denmark)

Pedersen, Bo Rohde; Agerkvist, Finn T.

2008-01-01

Simulations of a 6½-inch loudspeaker unit are performed and compared with a displacement measurement. The non-linear loudspeaker model is based on the major nonlinear functions and expanded with time-varying suspension behaviour and flux modulation. The results are presented with FFT plots of thr...... frequencies and different displacement levels. The model errors are discussed and analysed including a test with loudspeaker unit where the diaphragm is removed....

FDTD modeling of anisotropic nonlinear optical phenomena in silicon waveguides.

Science.gov (United States)

Dissanayake, Chethiya M; Premaratne, Malin; Rukhlenko, Ivan D; Agrawal, Govind P

2010-09-27

A deep insight into the inherent anisotropic optical properties of silicon is required to improve the performance of silicon-waveguide-based photonic devices. It may also lead to novel device concepts and substantially extend the capabilities of silicon photonics in the future. In this paper, for the first time to the best of our knowledge, we present a three-dimensional finite-difference time-domain (FDTD) method for modeling optical phenomena in silicon waveguides, which takes into account fully the anisotropy of the third-order electronic and Raman susceptibilities. We show that, under certain realistic conditions that prevent generation of the longitudinal optical field inside the waveguide, this model is considerably simplified and can be represented by a computationally efficient algorithm, suitable for numerical analysis of complex polarization effects. To demonstrate the versatility of our model, we study polarization dependence for several nonlinear effects, including self-phase modulation, cross-phase modulation, and stimulated Raman scattering. Our FDTD model provides a basis for a full-blown numerical simulator that is restricted neither by the single-mode assumption nor by the slowly varying envelope approximation.

Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in photoelectrochemical (PEC) solar cells

Science.gov (United States)

Harmon, Michael; Gamba, Irene M.; Ren, Kui

2016-12-01

This work concerns the numerical solution of a coupled system of self-consistent reaction-drift-diffusion-Poisson equations that describes the macroscopic dynamics of charge transport in photoelectrochemical (PEC) solar cells with reactive semiconductor and electrolyte interfaces. We present three numerical algorithms, mainly based on a mixed finite element and a local discontinuous Galerkin method for spatial discretization, with carefully chosen numerical fluxes, and implicit-explicit time stepping techniques, for solving the time-dependent nonlinear systems of partial differential equations. We perform computational simulations under various model parameters to demonstrate the performance of the proposed numerical algorithms as well as the impact of these parameters on the solution to the model.

Dynamic modeling of geometrically nonlinear electrostatically actuated microbeams (Corotational Finite Element formulation and analysis)

Energy Technology Data Exchange (ETDEWEB)

Borhan, H; Ahmadian, M T [Sharif University of Technology, Center of Excellence for Design, Robotics and Automation, School of Mechanical Engineering, PO Box 11365-9567, Tehran (Iran, Islamic Republic of)

2006-04-01

In this paper, a complete nonlinear finite element model for coupled-domain MEMS devices with electrostatic actuation and squeeze film effect is developed. For this purpose, a corotational finite element formulation for the dynamic analysis of planer Euler beams is employed. In this method, the internal nodal forces due to deformation and intrinsic residual stresses, the inertial nodal forces, and the damping effect of squeezed air film are systematically derived by consistent linearization of the fully geometrically nonlinear beam theory using d'Alamber and virtual work principles. An incremental-iterative method based on the Newmark direct integration procedure and the Newton-Raphson algorithm is used to solve the nonlinear dynamic equilibrium equations. Numerical examples are presented and compared with experimental findings which indicate properly good agreement.

Technogenic Rock Dumps Physical Properties' Prognosis via Results of the Structure Numerical Modeling

Directory of Open Access Journals (Sweden)

Markov Sergey

2017-01-01

Full Text Available Understanding of internal structure of the technogenic rock dumps (gob dumps is required condition for estimation of using ones as filtration massifs for treatment of mine wastewater. Internal structure of gob piles greatly depends on dumping technology to applying restrictions for use them as filtration massifs. Numerical modelling of gob dumps allows adequately estimate them physical parameters, as a filtration coefficient, density, etc. The gob dumps numerical modelling results given in this article, in particular was examined grain size distribution of determined fractions depend on dump height. Shown, that filtration coefficient is in a nonlinear dependence on amount of several fractions of rock in gob dump. The numerical model adequacy both the gob structure and the dependence of filtration coefficient from gob height acknowledged equality of calculated and real filtration coefficient values. The results of this research can be apply to peripheral dumping technology.

Numerical Analysis of a Class of THM Coupled Model for Porous Materials

Science.gov (United States)

Liu, Tangwei; Zhou, Jingying; Lu, Hongzhi

2018-01-01

We consider the coupled models of the Thermo-hydro-mechanical (THM) problem for porous materials which arises in many engineering applications. Firstly, mathematical models of the THM coupled problem for porous materials were discussed. Secondly, for different cases, some numerical difference schemes of coupled model were constructed, respectively. Finally, aassuming that the original water vapour effect is neglectable and that the volume fraction of liquid phase and the solid phase are constants, the nonlinear equations can be reduced to linear equations. The discrete equations corresponding to the linear equations were solved by the Arnodli method.

The basic approach to age-structured population dynamics models, methods and numerics

CERN Document Server

Iannelli, Mimmo

2017-01-01

This book provides an introduction to age-structured population modeling which emphasises the connection between mathematical theory and underlying biological assumptions. Through the rigorous development of the linear theory and the nonlinear theory alongside numerics, the authors explore classical equations that describe the dynamics of certain ecological systems. Modeling aspects are discussed to show how relevant problems in the fields of demography, ecology, and epidemiology can be formulated and treated within the theory. In particular, the book presents extensions of age-structured modelling to the spread of diseases and epidemics while also addressing the issue of regularity of solutions, the asymptotic behaviour of solutions, and numerical approximation. With sections on transmission models, non-autonomous models and global dynamics, this book fills a gap in the literature on theoretical population dynamics. The Basic Approach to Age-Structured Population Dynamics will appeal to graduate students an...

On the Numerical Modeling of Fluid Instabilities in the Small-Amplitude Limit

Science.gov (United States)

Zalesak, Steven; Schmitt, A. J.; Velikovich, A. L.

2008-11-01

The problem we wish to address is that of accurately modeling the evolution of small-amplitude perturbations to a time- dependent flow, where the unperturbed flow itself exhibits large-amplitude temporal and spatial variations. In particular, we wish to accurately model the evolution of small-amplitude perturbations to an imploding ICF pellet, which is subject to both Richtmyer-Meshkov and Rayleigh-Taylor instabilities. This modeling is difficult despite the expected linear evolution of the perturbations themselves, because these perturbations are embedded in a highly nonlinear, strongly- shocked, and highly complex flow field which in and of itself stresses numerical computation capabilities, and whose simulation often employs numerical techniques which were not designed with the proper treatment of small-amplitude perturbations in mind. We will review some of the techniques that we have found to be of use toward this end, including the imposition of a ``differentiability condition'' on the component numerical algorithms of the codes which implement such modeling, the appropriate representation of interfaces in an Eulerian hydrodynamics context, and the role of exact energy conservation.

A variational numerical method based on finite elements for the nonlinear solution characteristics of the periodically forced Chen system

Science.gov (United States)

Khan, Sabeel M.; Sunny, D. A.; Aqeel, M.

2017-09-01

Nonlinear dynamical systems and their solutions are very sensitive to initial conditions and therefore need to be approximated carefully. In this article, we present and analyze nonlinear solution characteristics of the periodically forced Chen system with the application of a variational method based on the concept of finite time-elements. Our approach is based on the discretization of physical time space into finite elements where each time-element is mapped to a natural time space. The solution of the system is then determined in natural time space using a set of suitable basis functions. The numerical algorithm is presented and implemented to compute and analyze nonlinear behavior at different time-step sizes. The obtained results show an excellent agreement with the classical RK-4 and RK-5 methods. The accuracy and convergence of the method is shown by comparing numerically computed results with the exact solution for a test problem. The presented method has shown a great potential in dealing with the solutions of nonlinear dynamical systems and thus can be utilized in delineating different features and characteristics of their solutions.

Mathematical modeling of zika virus disease with nonlinear incidence and optimal control

Science.gov (United States)

Goswami, Naba Kumar; Srivastav, Akhil Kumar; Ghosh, Mini; Shanmukha, B.

2018-04-01

The Zika virus was first discovered in a rhesus monkey in the Zika Forest of Uganda in 1947, and it was isolated from humans in Nigeria in 1952. Zika virus disease is primarily a mosquito-borne disease, which is transmitted to human primarily through the bite of an infected Aedes species mosquito. However, there is documented evidence of sexual transmission of this disease too. In this paper, a nonlinear mathematical model for Zika virus by considering nonlinear incidence is formulated and analyzed. The equilibria and the basic reproduction number (R0) of the model are found. The stability of the different equilibria of the model is discussed in detail. When the basic reproduction number R0 1, we have endemic equilibrium which is locally stable under some restriction on parameters. Further this model is extended to optimal control model and is analyzed by using Pontryagin’s Maximum Principle. It has been observed that optimal control plays a significant role in reducing the number of zika infectives. Finally, numerical simulation is performed to illustrate the analytical findings.

A Study of Enhanced, Higher Order Boussinesq-Type Equations and Their Numerical Modelling

DEFF Research Database (Denmark)

Banijamali, Babak

model is designated for the solution of higher-order Boussinesq-type equations, formulated in terms of the horizontal velocity at an arbitrary depth vector. Various discretisation techniques and grid definitions have been considered in this endeavour, undertaking a detailed analysis of the selected......This project has encompassed efforts in two separate veins: on the one hand, the acquiring of highly accurate model equations of the Boussinesq-type, and on the other hand, the theoretical and practical work in implementing such equations in the form of conventional numerical models, with obvious...... potential for applications to the realm of numerical modelling in coastal engineering. The derivation and analysis of several forms of higher-order in dispersion and non-linearity Boussinesq-type equations have been undertaken, obtaining and investigating the properties of a new and generalised class...

Numerical study of nonlinear singular fractional differential equations arising in biology by operational matrix of shifted Legendre polynomials

Directory of Open Access Journals (Sweden)

D. Jabari Sabeg

2016-10-01

Full Text Available In this paper, we present a new computational method for solving nonlinear singular boundary value problems of fractional order arising in biology. To this end, we apply the operational matrices of derivatives of shifted Legendre polynomials to reduce such problems to a system of nonlinear algebraic equations. To demonstrate the validity and applicability of the presented method, we present some numerical examples.

Numerical solutions for magnetohydrodynamic flow of nanofluid over a bidirectional non-linear stretching surface with prescribed surface heat flux boundary

International Nuclear Information System (INIS)

Mahanthesh, B.; Gireesha, B.J.; Gorla, R.S. Reddy; Abbasi, F.M.; Shehzad, S.A.

2016-01-01

Numerical solutions of three-dimensional flow over a non-linear stretching surface are developed in this article. An electrically conducting flow of viscous nanoliquid is considered. Heat transfer phenomenon is accounted under thermal radiation, Joule heating and viscous dissipation effects. We considered the variable heat flux condition at the surface of sheet. The governing mathematical equations are reduced to nonlinear ordinary differential systems through suitable dimensionless variables. A well-known shooting technique is implemented to obtain the results of dimensionless velocities and temperature. The obtained results are plotted for multiple values of pertinent parameters to discuss the salient features of these parameters on fluid velocity and temperature. The expressions of skin-friction coefficient and Nusselt number are computed and analyzed comprehensively through numerical values. A comparison of present results with the previous results in absence of nanoparticle volume fraction, mixed convection and magnetic field is computed and an excellent agreement noticed. We also computed the results for both linear and non-linear stretching sheet cases. - Highlights: • Hydromagnetic flow of nanofluid over a bidirectional non-linear stretching surface is examined. • Cu, Al 2 O3 and TiO 2 types nanoparticles are taken into account. • Numerical solutions have been computed and addressed. • The values of skin-friction and Nusselt number are presented.

Numerical models for high beta magnetohydrodynamic flow

International Nuclear Information System (INIS)

Brackbill, J.U.

1987-01-01

The fundamentals of numerical magnetohydrodynamics for highly conducting, high-beta plasmas are outlined. The discussions emphasize the physical properties of the flow, and how elementary concepts in numerical analysis can be applied to the construction of finite difference approximations that capture these features. The linear and nonlinear stability of explicit and implicit differencing in time is examined, the origin and effect of numerical diffusion in the calculation of convective transport is described, and a technique for maintaining solenoidality in the magnetic field is developed. Many of the points are illustrated by numerical examples. The techniques described are applicable to the time-dependent, high-beta flows normally encountered in magnetically confined plasmas, plasma switches, and space and astrophysical plasmas. 40 refs

Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method

Energy Technology Data Exchange (ETDEWEB)

Yang, Xiaofeng, E-mail: xfyang@math.sc.edu [Department of Mathematics, University of South Carolina, Columbia, SC 29208 (United States); Zhao, Jia, E-mail: zhao62@math.sc.edu [Department of Mathematics, University of South Carolina, Columbia, SC 29208 (United States); Department of Mathematics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599 (United States); Wang, Qi, E-mail: qwang@math.sc.edu [Department of Mathematics, University of South Carolina, Columbia, SC 29208 (United States); Beijing Computational Science Research Center, Beijing (China); School of Materials Science and Engineering, Nankai University, Tianjin (China)

2017-03-15

The Molecular Beam Epitaxial model is derived from the variation of a free energy, that consists of either a fourth order Ginzburg–Landau double well potential or a nonlinear logarithmic potential in terms of the gradient of a height function. One challenge in solving the MBE model numerically is how to develop proper temporal discretization for the nonlinear terms in order to preserve energy stability at the time-discrete level. In this paper, we resolve this issue by developing a first and second order time-stepping scheme based on the “Invariant Energy Quadratization” (IEQ) method. The novelty is that all nonlinear terms are treated semi-explicitly, and the resulted semi-discrete equations form a linear system at each time step. Moreover, the linear operator is symmetric positive definite and thus can be solved efficiently. We then prove that all proposed schemes are unconditionally energy stable. The semi-discrete schemes are further discretized in space using finite difference methods and implemented on GPUs for high-performance computing. Various 2D and 3D numerical examples are presented to demonstrate stability and accuracy of the proposed schemes.

Non-linear Model Predictive Control for cooling strings of superconducting magnets using superfluid helium

CERN Document Server

AUTHOR|(SzGeCERN)673023; Blanco Viñuela, Enrique

In each of eight arcs of the 27 km circumference Large Hadron Collider (LHC), 2.5 km long strings of super-conducting magnets are cooled with superfluid Helium II at 1.9 K. The temperature stabilisation is a challenging control problem due to complex non-linear dynamics of the magnets temperature and presence of multiple operational constraints. Strong nonlinearities and variable dead-times of the dynamics originate at strongly heat-flux dependent effective heat conductivity of superfluid that varies three orders of magnitude over the range of possible operational conditions. In order to improve the temperature stabilisation, a proof of concept on-line economic output-feedback Non-linear Model Predictive Controller (NMPC) is presented in this thesis. The controller is based on a novel complex first-principles distributed parameters numerical model of the temperature dynamics over a 214 m long sub-sector of the LHC that is characterized by very low computational cost of simulation needed in real-time optimizat...

A non-Linear transport model for determining shale rock characteristics

Science.gov (United States)

Ali, Iftikhar; Malik, Nadeem

2016-04-01

Unconventional hydrocarbon reservoirs consist of tight porous rocks which are characterised by nano-scale size porous networks with ultra-low permeability [1,2]. Transport of gas through them is not well understood at the present time, and realistic transport models are needed in order to determine rock properties and for estimating future gas pressure distribution in the reservoirs. Here, we consider a recently developed non-linear gas transport equation [3], ∂p-+ U ∂p- = D ∂2p-, t > 0, (1) ∂t ∂x ∂x2 complimented with suitable initial and boundary conditions, in order to determine shale rock properties such as the permeability K, the porosity φ and the tortuosity, τ. In our new model, the apparent convection velocity, U = U(p,px), and the apparent diffusivity D = D(p), are both highly non-linear functions of the pressure. The model incorporate various flow regimes (slip, surface diffusion, transition, continuum) based upon the Knudsen number Kn, and also includes Forchchiemers turbulence correction terms. In application, the model parameters and associated compressibility factors are fully pressure dependent, giving the model more realism than previous models. See [4]. Rock properties are determined by solving an inverse problem, with model parameters adjustment to minimise the error between the model simulation and available data. It is has been found that the proposed model performs better than previous models. Results and details of the model will be presented at the conference. Corresponding author: namalik@kfupm.edu.sa and nadeem_malik@cantab.net References [1] Cui, X., Bustin, A.M. and Bustin, R., "Measurements of gas permeability and diffusivity of tight reservoir rocks: different approaches and their applications", Geofluids 9, 208-223 (2009). [2] Chiba R., Fomin S., Chugunov V., Niibori Y. and Hashida T., "Numerical Simulation of Non Fickian Diffusion and Advection in a Fractured Porous Aquifer", AIP Conference Proceedings 898, 75 (2007

Bäcklund transformation, analytic soliton solutions and numerical simulation for a (2+1)-dimensional complex Ginzburg-Landau equation in a nonlinear fiber

Science.gov (United States)

Yu, Ming-Xiao; Tian, Bo; Chai, Jun; Yin, Hui-Min; Du, Zhong

2017-10-01

In this paper, we investigate a nonlinear fiber described by a (2+1)-dimensional complex Ginzburg-Landau equation with the chromatic dispersion, optical filtering, nonlinear and linear gain. Bäcklund transformation in the bilinear form is constructed. With the modified bilinear method, analytic soliton solutions are obtained. For the soliton, the amplitude can decrease or increase when the absolute value of the nonlinear or linear gain is enlarged, and the width can be compressed or amplified when the absolute value of the chromatic dispersion or optical filtering is enhanced. We study the stability of the numerical solutions numerically by applying the increasing amplitude, embedding the white noise and adding the Gaussian pulse to the initial values based on the analytic solutions, which shows that the numerical solutions are stable, not influenced by the finite initial perturbations.

Exact modelling of the optical bistability in ferroelectics via two-wave mixing: A system with full nonlinearity

Science.gov (United States)

Khushaini, Muhammad Asif A.; Ibrahim, Abdel-Baset M. A.; Choudhury, P. K.

2018-05-01

In this paper, we provide a complete mathematical model of the phenomenon of optical bistability (OB) resulting from the degenerate two-wave mixing (TWM) process of laser beams interacting with a single nonlinear layer of ferroelectric material. Starting with the electromagnetic wave equation for optical wave propagating in nonlinear media, a nonlinear coupled wave (CW) system with both self-phase modulation (SPM) and cross-phase modulation (XPM) sources of nonlinearity are derived. The complete CW system with full nonlinearity is solved numerically and a comparison between both the cases of with and without SPM at various combinations of design parameters is given. Furthermore, to provide a reliable theoretical model for the OB via TWM process, the results obtained theoretically are compared with the available experimental data. We found that the nonlinear system without SPM fails to predict the bistable response at lower combinations of the input parameters. However, at relatively higher values, the solution without SPM shows a reduction in the switching contrast and period in the OB response. A comparison with the experimental results shows better agreement with the system with full nonlinearity.

Continuous limit of a crowd motion and herding model: Analysis and numerical simulations

KAUST Repository

Pietschmann, Jan-Frederik

2011-11-01

In this paper we study the continuum limit of a cellular automaton model used for simulating human crowds with herding behaviour. We derive a system of non-linear partial differential equations resembling the Keller-Segel model for chemotaxis, however with a non-monotone interaction. The latter has interesting consequences on the behaviour of the model\\'s solutions, which we highlight in its analysis. In particular we study the possibility of stationary states, the formation of clusters and explore their connection to congestion. We also introduce an efficient numerical simulation approach based on an appropriate hybrid discontinuous Galerkin method, which in particular allows flexible treatment of complicated geometries. Extensive numerical studies also provide a better understanding of the strengths and shortcomings of the herding model, in particular we examine trapping effects of crowds behind nonconvex obstacles. © American Institute of Mathematical Sciences.

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.

Modeling of Volatility with Non-linear Time Series Model

OpenAIRE

Kim Song Yon; Kim Mun Chol

2013-01-01

In this paper, non-linear time series models are used to describe volatility in financial time series data. To describe volatility, two of the non-linear time series are combined into form TAR (Threshold Auto-Regressive Model) with AARCH (Asymmetric Auto-Regressive Conditional Heteroskedasticity) error term and its parameter estimation is studied.

Monte Carlo simulations with Symanzik's improved actions in the lattice 0(3) non-linear sigma-model

International Nuclear Information System (INIS)

Berg, B.; Montvay, I.; Meyer, S.

1983-10-01

The scaling properties of the lattice 0(3) non-linear delta-model are studied. The mass-gap, energy-momentum dispersion, correlation functions are measured by numerical Monte Carlo methods. Symanzik's tree-level and 1-loop improved actions are compared to the standard (nearest neigbour) action. (orig.)

Numerical analysis targets

International Nuclear Information System (INIS)

Sollogoub, Pierre

2001-01-01

Numerical analyses are needed in different steps of the overall design process. Complex models or non-linear reactor core behaviour are important for qualification and/or comparison of results obtained. Adequate models and test should be defined. Fuel assembly, fuel row, and the complete core should be tested for seismic effects causing LOCA and flow-induced vibrations (FIV)

An extended continuum model considering optimal velocity change with memory and numerical tests

Science.gov (United States)

Qingtao, Zhai; Hongxia, Ge; Rongjun, Cheng

2018-01-01

In this paper, an extended continuum model of traffic flow is proposed with the consideration of optimal velocity changes with memory. The new model's stability condition and KdV-Burgers equation considering the optimal velocities change with memory are deduced through linear stability theory and nonlinear analysis, respectively. Numerical simulation is carried out to study the extended continuum model, which explores how optimal velocity changes with memory affected velocity, density and energy consumption. Numerical results show that when considering the effects of optimal velocity changes with memory, the traffic jams can be suppressed efficiently. Both the memory step and sensitivity parameters of optimal velocity changes with memory will enhance the stability of traffic flow efficiently. Furthermore, numerical results demonstrates that the effect of optimal velocity changes with memory can avoid the disadvantage of historical information, which increases the stability of traffic flow on road, and so it improve the traffic flow stability and minimize cars' energy consumptions.

Numerical modelling of two phase flow with hysteresis in heterogeneous porous media

Energy Technology Data Exchange (ETDEWEB)

Abreu, E. [Instituto Nacional de Matematica Pura e Aplicada (IMPA), Rio de Janeiro, RJ (Brazil); Furtado, F.; Pereira, F. [University of Wyoming, Laramie, WY (United States). Dept. of Mathematicsatics; Souza, G. [Universidade do Estado do Rio de Janeiro (UERJ), RJ (Brazil)

2008-07-01

Numerical simulators are necessary for the understanding of multiphase flow in porous media in order to optimize hydrocarbon recovery. In this work, the immiscible flow of two incompressible phases, a problem very common in waterflooding of petroleum reservoirs, is considered and numerical simulation techniques are presented. The system of equations which describe this type of flow form a coupled, highly nonlinear system of time-dependent partial differential equations (PDEs). The equation for the saturation of the invading fluid is a convection-dominated, degenerate parabolic PDE whose solutions typically exhibit sharp fronts (i.e., internal layers with strong gradients) and is very difficult to approximate numerically. It is well known that accurate modeling of convective and diffusive processes is one of the most daunting tasks in the numerical approximation of PDEs. Particularly difficult is the case where convection dominates diffusion. Specifically, we consider the injection problem for a model of two-phase (water/oil) flow in a core sample of porous rock, taking into account hysteresis effects in the relative permeability of the oil phase. (author)

A time-domain finite element model reduction method for viscoelastic linear and nonlinear systems

Directory of Open Access Journals (Sweden)

Antônio Marcos Gonçalves de Lima

Full Text Available AbstractMany authors have shown that the effective design of viscoelastic systems can be conveniently carried out by using modern mathematical models to represent the frequency- and temperature-dependent behavior of viscoelastic materials. However, in the quest for design procedures of real-word engineering structures, the large number of exact evaluations of the dynamic responses during iterative procedures, combined with the typically high dimensions of large finite element models, makes the numerical analysis very costly, sometimes unfeasible. It is especially true when the viscoelastic materials are used to reduce vibrations of nonlinear systems. As a matter of fact, which the resolution of the resulting nonlinear equations of motion with frequency- and temperature-dependent viscoelastic damping forces is an interesting, but hard-to-solve problem. Those difficulties motivate the present study, in which a time-domain condensation strategy of viscoelastic systems is addressed, where the viscoelastic behavior is modeled by using a four parameter fractional derivative model. After the discussion of various theoretical aspects, the exact and reduced time responses are calculated for a three-layer sandwich plate by considering nonlinear boundary conditions.

Controlled Nonlinear Stochastic Delay Equations: Part II: Approximations and Pipe-Flow Representations

International Nuclear Information System (INIS)

Kushner, Harold J.

2012-01-01

This is the second part of a work dealing with key issues that have not been addressed in the modeling and numerical optimization of nonlinear stochastic delay systems. We consider new classes of models, such as those with nonlinear functions of several controls (such as products), each with is own delay, controlled random Poisson measure driving terms, admissions control with delayed retrials, and others. Part I was concerned with issues concerning the class of admissible controls and their approximations, since the classical definitions are inadequate for our models. This part is concerned with transportation equation representations and their approximations. Such representations of nonlinear stochastic delay models have been crucial in the development of numerical algorithms with much reduced memory and computational requirements. The representations for the new models are not obvious and are developed. They also provide a template for the adaptation of the Markov chain approximation numerical methods.

Comparing coefficients of nested nonlinear probability models

DEFF Research Database (Denmark)

Kohler, Ulrich; Karlson, Kristian Bernt; Holm, Anders

2011-01-01

In a series of recent articles, Karlson, Holm and Breen have developed a method for comparing the estimated coeffcients of two nested nonlinear probability models. This article describes this method and the user-written program khb that implements the method. The KHB-method is a general decomposi......In a series of recent articles, Karlson, Holm and Breen have developed a method for comparing the estimated coeffcients of two nested nonlinear probability models. This article describes this method and the user-written program khb that implements the method. The KHB-method is a general...... decomposition method that is unaffected by the rescaling or attenuation bias that arise in cross-model comparisons in nonlinear models. It recovers the degree to which a control variable, Z, mediates or explains the relationship between X and a latent outcome variable, Y*, underlying the nonlinear probability...

Explicit formulation of second and third order optical nonlinearity in the FDTD framework

Science.gov (United States)

Varin, Charles; Emms, Rhys; Bart, Graeme; Fennel, Thomas; Brabec, Thomas

2018-01-01

The finite-difference time-domain (FDTD) method is a flexible and powerful technique for rigorously solving Maxwell's equations. However, three-dimensional optical nonlinearity in current commercial and research FDTD softwares requires solving iteratively an implicit form of Maxwell's equations over the entire numerical space and at each time step. Reaching numerical convergence demands significant computational resources and practical implementation often requires major modifications to the core FDTD engine. In this paper, we present an explicit method to include second and third order optical nonlinearity in the FDTD framework based on a nonlinear generalization of the Lorentz dispersion model. A formal derivation of the nonlinear Lorentz dispersion equation is equally provided, starting from the quantum mechanical equations describing nonlinear optics in the two-level approximation. With the proposed approach, numerical integration of optical nonlinearity and dispersion in FDTD is intuitive, transparent, and fully explicit. A strong-field formulation is also proposed, which opens an interesting avenue for FDTD-based modelling of the extreme nonlinear optics phenomena involved in laser filamentation and femtosecond micromachining of dielectrics.

Nonlinear dynamics analysis of a modified optically injected semiconductor lasers model

International Nuclear Information System (INIS)

Chu Yandong; Li Xianfeng; Zhang Jiangang; Chang Yingxiang

2009-01-01

In this paper, a new nonlinear autonomous system that was introduced by Chlouverakis and Sprott is studied further, to present very rich and complex nonlinear dynamical behaviors. Some basic dynamical properties are studied either analytically or numerically, such as Poincare mapping, Lyapunov exponents, fractal dimension, continuous power spectrum and so forth. Furthermore, the coexistence of different attractors is discovered on the Poincare maps. Meanwhile, chaotic oscillation of this system is converted into a stable periodic orbit with the method of time-delayed feedback, which demonstrated by numerical simulations and the robustness of this method is proved.

Toward a scalable flexible-order model for 3D nonlinear water waves

DEFF Research Database (Denmark)

Engsig-Karup, Allan Peter; Ducrozet, Guillaume; Bingham, Harry B.

For marine and coastal applications, current work are directed toward the development of a scalable numerical 3D model for fully nonlinear potential water waves over arbitrary depths. The model is high-order accurate, robust and efficient for large-scale problems, and support will be included...... for flexibility in the description of structures by the use of curvilinear boundary-fitted meshes. The mathematical equations for potential waves in the physical domain is transformed through $\\sigma$-mapping(s) to a time-invariant boundary-fitted domain which then becomes a basis for an efficient solution...... strategy on a time-invariant mesh. The 3D numerical model is based on a finite difference method as in the original works \\cite{LiFleming1997,BinghamZhang2007}. Full details and other aspects of an improved 3D solution can be found in \\cite{EBL08}. The new and improved approach for three...

Nonlinear Control of Heartbeat Models

Directory of Open Access Journals (Sweden)

Witt Thanom

2011-02-01

Full Text Available This paper presents a novel application of nonlinear control theory to heartbeat models. Existing heartbeat models are investigated and modified by incorporating the control input as a pacemaker to provide the control channel. A nonlinear feedback linearization technique is applied to force the output of the systems to generate artificial electrocardiogram (ECG signal using discrete data as the reference inputs. The synthetic ECG may serve as a flexible signal source to assess the effectiveness of a diagnostic ECG signal-processing device.

Continuous limit of a crowd motion and herding model: Analysis and numerical simulations

KAUST Repository

Pietschmann, Jan-Frederik; Markowich, Peter Alexander; Burger, Martin

2011-01-01

In this paper we study the continuum limit of a cellular automaton model used for simulating human crowds with herding behaviour. We derive a system of non-linear partial differential equations resembling the Keller-Segel model for chemotaxis, however with a non-monotone interaction. The latter has interesting consequences on the behaviour of the model's solutions, which we highlight in its analysis. In particular we study the possibility of stationary states, the formation of clusters and explore their connection to congestion. We also introduce an efficient numerical simulation approach based on an appropriate hybrid discontinuous Galerkin method, which in particular allows flexible treatment of complicated geometries. Extensive numerical studies also provide a better understanding of the strengths and shortcomings of the herding model, in particular we examine trapping effects of crowds behind nonconvex obstacles. © American Institute of Mathematical Sciences.

Nonlinear modeling and dynamic analysis of hydro-turbine governing system in the process of load rejection transient

International Nuclear Information System (INIS)

Zhang, Hao; Chen, Diyi; Xu, Beibei; Wang, Feifei

2015-01-01

Graphical abstract: Nonlinear dynamic transfer coefficients are introduced to the hydro-turbine governing system. In the process of load reject ion transient, the nonlinear dynamical behaviors of the system are studied in detail. - Highlights: • A novel mathematical model of a hydro-turbine governing system is established. • The process of load rejection transient is considered. • Nonlinear dynamic transfer coefficients are introduced to the system. • The bifurcation diagram with the variable t has better engineering significance. • The nonlinear dynamical behaviors of the system are studied in detail. - Abstract: This article pays attention to the mathematical modeling of a hydro-turbine governing system in the process of load rejection transient. As a pioneer work, the nonlinear dynamic transfer coefficients are introduced in a penstock system. Considering a generator system, a turbine system and a governor system, we present a novel nonlinear dynamical model of a hydro-turbine governing system. Fortunately, for the unchanged of PID parameters, we acquire the stable regions of the governing system in the process of load rejection transient by numerical simulations. Moreover, the nonlinear dynamic behaviors of the governing system are illustrated by bifurcation diagrams, Poincare maps, time waveforms and phase orbits. More importantly, these methods and analytic results will present theoretical groundwork for allowing a hydropower station in the process of load rejection transient

Non-linear numerical studies of the tearing mode

International Nuclear Information System (INIS)

Schnack, D.D. Jr.

1978-01-01

A non-linear, time dependent, hydromagnetic model is developed and applied to the tearing mode, one of a class of instabilities which can occur in a magnetically confined plasma when the constraint of infinite conductivity is relaxed. The model is based on the eight partial differential equations of resistive magnetohydrodynamics (MHD). The equations are expressed as a set of conservation laws which conserves magnetic flux, momentum, mass, and total energy. These equations are then written in general, orthogonal, curvilinear coordinates in two space dimensions, so that the model can readily be applied to a variety of geometries. No assumption about the ordering of terms is made. The resulting equations are then solved by the method of finite differences on an Eulerian mesh. The model is applied to several geometries

Numerical solutions for magnetohydrodynamic flow of nanofluid over a bidirectional non-linear stretching surface with prescribed surface heat flux boundary

Energy Technology Data Exchange (ETDEWEB)

Mahanthesh, B., E-mail: bmanths@gmail.com [Department of Mathematics, AIMS Institutes, Peenya, 560058 Bangalore (India); Department of Studies and Research in Mathematics, Kuvempu University, Shankaraghatta, 577451 Shimoga, Karnataka (India); Gireesha, B.J., E-mail: bjgireesu@rediffmail.com [Department of Studies and Research in Mathematics, Kuvempu University, Shankaraghatta, 577451 Shimoga, Karnataka (India); Department of Mechanical Engineering, Cleveland State University, Cleveland, OH (United States); Gorla, R.S. Reddy, E-mail: r.gorla@csuohio.edu [Department of Mechanical Engineering, Cleveland State University, Cleveland, OH (United States); Abbasi, F.M., E-mail: abbasisarkar@gmail.com [Department of Mathematics, Comsats Institute of Information Technology, Islamabad 44000 (Pakistan); Shehzad, S.A., E-mail: ali_qau70@yahoo.com [Department of Mathematics, Comsats Institute of Information Technology, Sahiwal 57000 (Pakistan)

2016-11-01

Numerical solutions of three-dimensional flow over a non-linear stretching surface are developed in this article. An electrically conducting flow of viscous nanoliquid is considered. Heat transfer phenomenon is accounted under thermal radiation, Joule heating and viscous dissipation effects. We considered the variable heat flux condition at the surface of sheet. The governing mathematical equations are reduced to nonlinear ordinary differential systems through suitable dimensionless variables. A well-known shooting technique is implemented to obtain the results of dimensionless velocities and temperature. The obtained results are plotted for multiple values of pertinent parameters to discuss the salient features of these parameters on fluid velocity and temperature. The expressions of skin-friction coefficient and Nusselt number are computed and analyzed comprehensively through numerical values. A comparison of present results with the previous results in absence of nanoparticle volume fraction, mixed convection and magnetic field is computed and an excellent agreement noticed. We also computed the results for both linear and non-linear stretching sheet cases. - Highlights: • Hydromagnetic flow of nanofluid over a bidirectional non-linear stretching surface is examined. • Cu, Al{sub 2}O3 and TiO{sub 2} types nanoparticles are taken into account. • Numerical solutions have been computed and addressed. • The values of skin-friction and Nusselt number are presented.

Nonlinear finite element modeling of concrete deep beams with openings strengthened with externally-bonded composites

International Nuclear Information System (INIS)

Hawileh, Rami A.; El-Maaddawy, Tamer A.; Naser, Mohannad Z.

2012-01-01

Highlights: ► A 3D nonlinear FE model is developed of RC deep beams with web openings. ► We used cohesion elements to simulate bond. ► The developed FE model is suitable for analysis of such complex structures. -- Abstract: This paper aims to develop 3D nonlinear finite element (FE) models for reinforced concrete (RC) deep beams containing web openings and strengthened in shear with carbon fiber reinforced polymer (CFRP) composite sheets. The web openings interrupted the natural load path either fully or partially. The FE models adopted realistic materials constitutive laws that account for the nonlinear behavior of materials. In the FE models, solid elements for concrete, multi-layer shell elements for CFRP and link elements for steel reinforcement were used to simulate the physical models. Special interface elements were implemented in the FE models to simulate the interfacial bond behavior between the concrete and CFRP composites. A comparison between the FE results and experimental data published in the literature demonstrated the validity of the computational models in capturing the structural response for both unstrengthened and CFRP-strengthened deep beams with openings. The developed FE models can serve as a numerical platform for performance prediction of RC deep beams with openings strengthened in shear with CFRP composites.

Analytical and numerical modelling of thermoviscous shocks in their interactions in nonlinear fluids including dissipation

DEFF Research Database (Denmark)

Rasmussen, Anders Rønne; Sørensen, Mads Peter; Gaididei, Yuri Borisovich

2010-01-01

A wave equation, that governs finite amplitude acoustic disturbances in a thermoviscous Newtonian fluid, and includes nonlinear terms up to second order, is proposed. The equation preserves the Hamiltonian structure of the fundamental fluid dynamical equations in the non dissipative limit. An exact...... thermoviscous shock solution is derived. This solution is, in an overall sense, equivalent to the Taylor shock solution of the Burgers equation. However, in contrast to the Burgers equation, the model equation considered here is capable to describe waves propagating in opposite directions. Studies of head...

Influence of Non-linear Radiation Heat Flux on Rotating Maxwell Fluid over a Deformable Surface: A Numerical Study

Science.gov (United States)

Mustafa, M.; Mushtaq, A.; Hayat, T.; Alsaedi, A.

2018-04-01

Mathematical model for Maxwell fluid flow in rotating frame induced by an isothermal stretching wall is explored numerically. Scale analysis based boundary layer approximations are applied to simplify the conservation relations which are later converted to similar forms via appropriate substitutions. A numerical approach is utilized to derive similarity solutions for broad range of Deborah number. The results predict that velocity distributions are inversely proportional to the stress relaxation time. This outcome is different from that observed for the elastic parameter of second grade fluid. Unlike non-rotating frame, the solution curves are oscillatory decaying functions of similarity variable. As angular velocity enlarges, temperature rises and significant drop in the heat transfer coefficient occurs. We note that the wall slope of temperature has an asymptotically decaying profile against the wall to ambient ratio parameter. From the qualitative view point, temperature ratio parameter and radiation parameter have similar effect on the thermal boundary layer. Furthermore, radiation parameter has a definite role in improving the cooling process of the stretching boundary. A comparative study of current numerical computations and those from the existing studies is also presented in a limiting case. To our knowledge, the phenomenon of non-linear radiation in rotating viscoelastic flow due to linearly stretched plate is just modeled here.

Nonlinear flight dynamics and stability of hovering model insects

Science.gov (United States)

Liang, Bin; Sun, Mao

2013-01-01

Current analyses on insect dynamic flight stability are based on linear theory and limited to small disturbance motions. However, insects' aerial environment is filled with swirling eddies and wind gusts, and large disturbances are common. Here, we numerically solve the equations of motion coupled with the Navier–Stokes equations to simulate the large disturbance motions and analyse the nonlinear flight dynamics of hovering model insects. We consider two representative model insects, a model hawkmoth (large size, low wingbeat frequency) and a model dronefly (small size, high wingbeat frequency). For small and large initial disturbances, the disturbance motion grows with time, and the insects tumble and never return to the equilibrium state; the hovering flight is inherently (passively) unstable. The instability is caused by a pitch moment produced by forward/backward motion and/or a roll moment produced by side motion of the insect. PMID:23697714

Nonlinear structural mechanics theory, dynamical phenomena and modeling

CERN Document Server

Lacarbonara, Walter

2013-01-01

Nonlinear Structural Mechanics: Theory, Dynamical Phenomena and Modeling offers a concise, coherent presentation of the theoretical framework of nonlinear structural mechanics, computational methods, applications, parametric investigations of nonlinear phenomena and their mechanical interpretation towards design. The theoretical and computational tools that enable the formulation, solution, and interpretation of nonlinear structures are presented in a systematic fashion so as to gradually attain an increasing level of complexity of structural behaviors, under the prevailing assumptions on the geometry of deformation, the constitutive aspects and the loading scenarios. Readers will find a treatment of the foundations of nonlinear structural mechanics towards advanced reduced models, unified with modern computational tools in the framework of the prominent nonlinear structural dynamic phenomena while tackling both the mathematical and applied sciences. Nonlinear Structural Mechanics: Theory, Dynamical Phenomena...

Threshold Dynamics in Stochastic SIRS Epidemic Models with Nonlinear Incidence and Vaccination

Directory of Open Access Journals (Sweden)

Lei Wang

2017-01-01

Full Text Available In this paper, the dynamical behaviors for a stochastic SIRS epidemic model with nonlinear incidence and vaccination are investigated. In the models, the disease transmission coefficient and the removal rates are all affected by noise. Some new basic properties of the models are found. Applying these properties, we establish a series of new threshold conditions on the stochastically exponential extinction, stochastic persistence, and permanence in the mean of the disease with probability one for the models. Furthermore, we obtain a sufficient condition on the existence of unique stationary distribution for the model. Finally, a series of numerical examples are introduced to illustrate our main theoretical results and some conjectures are further proposed.

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

Analysis of nonlinear systems using ARMA [autoregressive moving average] models

International Nuclear Information System (INIS)

Hunter, N.F. Jr.

1990-01-01

While many vibration systems exhibit primarily linear behavior, a significant percentage of the systems encountered in vibration and model testing are mildly to severely nonlinear. Analysis methods for such nonlinear systems are not yet well developed and the response of such systems is not accurately predicted by linear models. Nonlinear ARMA (autoregressive moving average) models are one method for the analysis and response prediction of nonlinear vibratory systems. In this paper we review the background of linear and nonlinear ARMA models, and illustrate the application of these models to nonlinear vibration systems. We conclude by summarizing the advantages and disadvantages of ARMA models and emphasizing prospects for future development. 14 refs., 11 figs

On the kinematic criterion for the inception of breaking in surface gravity waves: Fully nonlinear numerical simulations and experimental verification

Science.gov (United States)

Khait, A.; Shemer, L.

2018-05-01

The evolution of unidirectional wave trains containing a wave that gradually becomes steep is evaluated experimentally and numerically using the Boundary Element Method (BEM). The boundary conditions for the nonlinear numerical simulations corresponded to the actual movements of the wavemaker paddle as recorded in the physical experiments, allowing direct comparison between the measured in experiments' characteristics of the wave train and the numerical predictions. The high level of qualitative and quantitative agreement between the measurements and simulations validated the kinematic criterion for the inception of breaking and the location of the spilling breaker, on the basis of the BEM computations and associated experiments. The breaking inception is associated with the fluid particle at the crest of the steep wave that has been accelerated to match and surpass the crest velocity. The previously observed significant slow-down of the crest while approaching breaking is verified numerically; both narrow-/broad-banded wave trains are considered. Finally, the relative importance of linear and nonlinear contributions is analyzed.

Nonlinear bioheat transfer models and multi-objective numerical optimization of the cryosurgery operations

Energy Technology Data Exchange (ETDEWEB)

Kudryashov, Nikolay A.; Shilnikov, Kirill E. [National Research Nuclear University MEPhI, Department of Applied Mathematics, Moscow (Russian Federation)

2016-06-08

Numerical computation of the three dimensional problem of the freezing interface propagation during the cryosurgery coupled with the multi-objective optimization methods is used in order to improve the efficiency and safety of the cryosurgery operations performing. Prostate cancer treatment and cutaneous cryosurgery are considered. The heat transfer in soft tissue during the thermal exposure to low temperature is described by the Pennes bioheat model and is coupled with an enthalpy method for blurred phase change computations. The finite volume method combined with the control volume approximation of the heat fluxes is applied for the cryosurgery numerical modeling on the tumor tissue of a quite arbitrary shape. The flux relaxation approach is used for the stability improvement of the explicit finite difference schemes. The method of the additional heating elements mounting is studied as an approach to control the cellular necrosis front propagation. Whereas the undestucted tumor tissue and destucted healthy tissue volumes are considered as objective functions, the locations of additional heating elements in cutaneous cryosurgery and cryotips in prostate cancer cryotreatment are considered as objective variables in multi-objective problem. The quasi-gradient method is proposed for the searching of the Pareto front segments as the multi-objective optimization problem solutions.

Nonlinear plasma wave models in 3D fluid simulations of laser-plasma interaction

Science.gov (United States)

Chapman, Thomas; Berger, Richard; Arrighi, Bill; Langer, Steve; Banks, Jeffrey; Brunner, Stephan

2017-10-01

Simulations of laser-plasma interaction (LPI) in inertial confinement fusion (ICF) conditions require multi-mm spatial scales due to the typical laser beam size and durations of order 100 ps in order for numerical laser reflectivities to converge. To be computationally achievable, these scales necessitate a fluid-like treatment of light and plasma waves with a spatial grid size on the order of the light wave length. Plasma waves experience many nonlinear phenomena not naturally described by a fluid treatment, such as frequency shifts induced by trapping, a nonlinear (typically suppressed) Landau damping, and mode couplings leading to instabilities that can cause the plasma wave to decay rapidly. These processes affect the onset and saturation of stimulated Raman and Brillouin scattering, and are of direct interest to the modeling and prediction of deleterious LPI in ICF. It is not currently computationally feasible to simulate these Debye length-scale phenomena in 3D across experimental scales. Analytically-derived and/or numerically benchmarked models of processes occurring at scales finer than the fluid simulation grid offer a path forward. We demonstrate the impact of a range of kinetic processes on plasma reflectivity via models included in the LPI simulation code pF3D. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

A discontinuous Galerkin approach for conservative modeling of fully nonlinear and weakly dispersive wave transformations

Science.gov (United States)

Sharifian, Mohammad Kazem; Kesserwani, Georges; Hassanzadeh, Yousef

2018-05-01

This work extends a robust second-order Runge-Kutta Discontinuous Galerkin (RKDG2) method to solve the fully nonlinear and weakly dispersive flows, within a scope to simultaneously address accuracy, conservativeness, cost-efficiency and practical needs. The mathematical model governing such flows is based on a variant form of the Green-Naghdi (GN) equations decomposed as a hyperbolic shallow water system with an elliptic source term. Practical features of relevance (i.e. conservative modeling over irregular terrain with wetting and drying and local slope limiting) have been restored from an RKDG2 solver to the Nonlinear Shallow Water (NSW) equations, alongside new considerations to integrate elliptic source terms (i.e. via a fourth-order local discretization of the topography) and to enable local capturing of breaking waves (i.e. via adding a detector for switching off the dispersive terms). Numerical results are presented, demonstrating the overall capability of the proposed approach in achieving realistic prediction of nearshore wave processes involving both nonlinearity and dispersion effects within a single model.

Nonlinear State Space Modeling and System Identification for Electrohydraulic Control

Directory of Open Access Journals (Sweden)

Jun Yan

2013-01-01

Full Text Available The paper deals with nonlinear modeling and identification of an electrohydraulic control system for improving its tracking performance. We build the nonlinear state space model for analyzing the highly nonlinear system and then develop a Hammerstein-Wiener (H-W model which consists of a static input nonlinear block with two-segment polynomial nonlinearities, a linear time-invariant dynamic block, and a static output nonlinear block with single polynomial nonlinearity to describe it. We simplify the H-W model into a linear-in-parameters structure by using the key term separation principle and then use a modified recursive least square method with iterative estimation of internal variables to identify all the unknown parameters simultaneously. It is found that the proposed H-W model approximates the actual system better than the independent Hammerstein, Wiener, and ARX models. The prediction error of the H-W model is about 13%, 54%, and 58% less than the Hammerstein, Wiener, and ARX models, respectively.

Seismic analysis of piping with nonlinear supports

International Nuclear Information System (INIS)

Barta, D.A.; Huang, S.N.; Severud, L.K.

1980-01-01

The modeling and results of nonlinear time-history seismic analyses for three sizes of pipelines restrained by mechanical snubbes are presented. Numerous parametric analyses were conducted to obtain sensitivity information which identifies relative importance of the model and analysis ingredients. Special considerations for modeling the pipe clamps and the mechanical snubbers based on experimental characterization data are discussed. Comparisions are also given of seismic responses, loads and pipe stresses predicted by standard response spectra methods and the nonlinear time-history methods

Nonlinear Poisson equation for heterogeneous media.

Science.gov (United States)

Hu, Langhua; Wei, Guo-Wei

2012-08-22

The Poisson equation is a widely accepted model for electrostatic analysis. However, the Poisson equation is derived based on electric polarizations in a linear, isotropic, and homogeneous dielectric medium. This article introduces a nonlinear Poisson equation to take into consideration of hyperpolarization effects due to intensive charges and possible nonlinear, anisotropic, and heterogeneous media. Variational principle is utilized to derive the nonlinear Poisson model from an electrostatic energy functional. To apply the proposed nonlinear Poisson equation for the solvation analysis, we also construct a nonpolar solvation energy functional based on the nonlinear Poisson equation by using the geometric measure theory. At a fixed temperature, the proposed nonlinear Poisson theory is extensively validated by the electrostatic analysis of the Kirkwood model and a set of 20 proteins, and the solvation analysis of a set of 17 small molecules whose experimental measurements are also available for a comparison. Moreover, the nonlinear Poisson equation is further applied to the solvation analysis of 21 compounds at different temperatures. Numerical results are compared to theoretical prediction, experimental measurements, and those obtained from other theoretical methods in the literature. A good agreement between our results and experimental data as well as theoretical results suggests that the proposed nonlinear Poisson model is a potentially useful model for electrostatic analysis involving hyperpolarization effects. Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.

Numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow

KAUST Repository

Zhu, Guangpu

2018-04-17

In this paper, we consider the numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow. The nonlinearly coupled model consists of two Cahn-Hilliard type equations and incompressible Navier-Stokes equations. Using the Invariant Energy Quadratization (IEQ) approach, the governing system is transformed into an equivalent form, which allows the nonlinear potentials to be treated efficiently and semi-explicitly. we construct a first and a second-order time marching schemes, which are extremely efficient and easy-to-implement, for the transformed governing system. At each time step, the schemes involve solving a sequence of linear elliptic equations, and computations of phase variables, velocity and pressure are totally decoupled. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient and unconditionally energy stable. Using our schemes, we investigate the effect of surfactants on droplet deformation and collision under a shear flow. The increase of surfactant concentration can enhance droplet deformation and inhibit droplet coalescence.

Numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow

KAUST Repository

Zhu, Guangpu; Kou, Jisheng; Sun, Shuyu; Yao, Jun; Li, Aifen

2018-01-01

In this paper, we consider the numerical approximation of a binary fluid-surfactant phase field model of two-phase incompressible flow. The nonlinearly coupled model consists of two Cahn-Hilliard type equations and incompressible Navier-Stokes equations. Using the Invariant Energy Quadratization (IEQ) approach, the governing system is transformed into an equivalent form, which allows the nonlinear potentials to be treated efficiently and semi-explicitly. we construct a first and a second-order time marching schemes, which are extremely efficient and easy-to-implement, for the transformed governing system. At each time step, the schemes involve solving a sequence of linear elliptic equations, and computations of phase variables, velocity and pressure are totally decoupled. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient and unconditionally energy stable. Using our schemes, we investigate the effect of surfactants on droplet deformation and collision under a shear flow. The increase of surfactant concentration can enhance droplet deformation and inhibit droplet coalescence.

The cultural implications of growth: Modeling nonlinear interaction of trait selection and population dynamics.

Science.gov (United States)

Antoci, Angelo; Galeotti, Marcello; Russu, Paolo; Luigi Sacco, Pier

2018-05-01

In this paper, we study a nonlinear model of the interaction between trait selection and population dynamics, building on previous work of Ghirlanda et al. [Theor. Popul. Biol. 77, 181-188 (2010)] and Antoci et al. [Commun. Nonlinear Sci. Numer. Simul. 58, 92-106 (2018)]. We establish some basic properties of the model dynamics and present some simulations of the fine-grained structure of alternative dynamic regimes for chosen combinations of parameters. The role of the parameters that govern the reinforcement/corruption of maladaptive vs. adaptive traits is of special importance in determining the model's dynamic evolution. The main implication of this result is the need to pay special attention to the structural forces that may favor the emergence and consolidation of maladaptive traits in contemporary socio-economies, as it is the case, for example, for the stimulation of dysfunctional consumption habits and lifestyles in the pursuit of short-term profits.

The cultural implications of growth: Modeling nonlinear interaction of trait selection and population dynamics

Science.gov (United States)

Antoci, Angelo; Galeotti, Marcello; Russu, Paolo; Luigi Sacco, Pier

2018-05-01

In this paper, we study a nonlinear model of the interaction between trait selection and population dynamics, building on previous work of Ghirlanda et al. [Theor. Popul. Biol. 77, 181-188 (2010)] and Antoci et al. [Commun. Nonlinear Sci. Numer. Simul. 58, 92-106 (2018)]. We establish some basic properties of the model dynamics and present some simulations of the fine-grained structure of alternative dynamic regimes for chosen combinations of parameters. The role of the parameters that govern the reinforcement/corruption of maladaptive vs. adaptive traits is of special importance in determining the model's dynamic evolution. The main implication of this result is the need to pay special attention to the structural forces that may favor the emergence and consolidation of maladaptive traits in contemporary socio-economies, as it is the case, for example, for the stimulation of dysfunctional consumption habits and lifestyles in the pursuit of short-term profits.

Numerical modeling of the 2017 active seismic infrasound balloon experiment

Science.gov (United States)

Brissaud, Q.; Komjathy, A.; Garcia, R.; Cutts, J. A.; Pauken, M.; Krishnamoorthy, S.; Mimoun, D.; Jackson, J. M.; Lai, V. H.; Kedar, S.; Levillain, E.

2017-12-01

We have developed a numerical tool to propagate acoustic and gravity waves in a coupled solid-fluid medium with topography. It is a hybrid method between a continuous Galerkin and a discontinuous Galerkin method that accounts for non-linear atmospheric waves, visco-elastic waves and topography. We apply this method to a recent experiment that took place in the Nevada desert to study acoustic waves from seismic events. This experiment, developed by JPL and its partners, wants to demonstrate the viability of a new approach to probe seismic-induced acoustic waves from a balloon platform. To the best of our knowledge, this could be the only way, for planetary missions, to perform tomography when one faces challenging surface conditions, with high pressure and temperature (e.g. Venus), and thus when it is impossible to use conventional electronics routinely employed on Earth. To fully demonstrate the effectiveness of such a technique one should also be able to reconstruct the observed signals from numerical modeling. To model the seismic hammer experiment and the subsequent acoustic wave propagation, we rely on a subsurface seismic model constructed from the seismometers measurements during the 2017 Nevada experiment and an atmospheric model built from meteorological data. The source is considered as a Gaussian point source located at the surface. Comparison between the numerical modeling and the experimental data could help future mission designs and provide great insights into the planet's interior structure.

Nonlinear Convective Models of RR Lyrae Stars

Science.gov (United States)

Feuchtinger, M.; Dorfi, E. A.

The nonlinear behavior of RR Lyrae pulsations is investigated using a state-of-the-art numerical technique solving the full time-dependent system of radiation hydrodynamics. Grey radiative transfer is included by a variable Eddington-factor method and we use the time-dependent turbulent convection model according to Kuhfuss (1986, A&A 160, 116) in the version of Wuchterl (1995, Comp. Phys. Comm. 89, 19). OPAL opacities extended by the Alexander molecule opacities at temperatures below 6000 K and an equation of state according to Wuchterl (1990, A&A 238, 83) close the system. The resulting nonlinear system is discretized on an adaptive mesh developed by Dorfi & Drury (1987, J. Comp. Phys. 69, 175), which is important to provide the necessary spatial resolution in critical regions like ionization zones and shock waves. Additionally, we employ a second order advection scheme, a time centered temporal discretizaton and an artificial tensor viscosity in order to treat discontinuities. We compute fundamental as well first overtone models of RR Lyrae stars for a grid of stellar parameters both with and without convective energy transport in order to give a detailed picture of the pulsation-convection interaction. In order to investigate the influence of the different features of the convection model calculations with and without overshooting, turbulent pressure and turbulent viscosity are performed and compared with each other. A standard Fourier decomposition is used to confront the resulting light and radial velocity variations with recent observations and we show that the well known RR Lyrae phase discrepancy problem (Simon 1985, ApJ 299, 723) can be resolved with these stellar pulsation computations.

Exact travelling wave solutions for some important nonlinear

Indian Academy of Sciences (India)

The two-dimensional nonlinear physical models and coupled nonlinear systems such as Maccari equations, Higgs equations and Schrödinger–KdV equations have been widely applied in many branches of physics. So, finding exact travelling wave solutions of such equations are very helpful in the theories and numerical ...

Simulation of nonlinear wave run-up with a high-order Boussinesq model

DEFF Research Database (Denmark)

Fuhrman, David R.; Madsen, Per A.

2008-01-01

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions....... As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions...... cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability...

NONLINEAR PLANT PIECEWISE-CONTINUOUS MODEL MATRIX PARAMETERS ESTIMATION

Directory of Open Access Journals (Sweden)

Roman L. Leibov

2017-09-01

Full Text Available This paper presents a nonlinear plant piecewise-continuous model matrix parameters estimation technique using nonlinear model time responses and random search method. One of piecewise-continuous model application areas is defined. The results of proposed approach application for aircraft turbofan engine piecewisecontinuous model formation are presented

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.

Simulation of nonlinear random vibrations using artificial neural networks

Energy Technology Data Exchange (ETDEWEB)

Paez, T.L.; Tucker, S.; O`Gorman, C.

1997-02-01

The simulation of mechanical system random vibrations is important in structural dynamics, but it is particularly difficult when the system under consideration is nonlinear. Artificial neural networks provide a useful tool for the modeling of nonlinear systems, however, such modeling may be inefficient or insufficiently accurate when the system under consideration is complex. This paper shows that there are several transformations that can be used to uncouple and simplify the components of motion of a complex nonlinear system, thereby making its modeling and random vibration simulation, via component modeling with artificial neural networks, a much simpler problem. A numerical example is presented.

Mathematical and numerical analysis of a few hydrodynamic and kinetic models of plasma physics

International Nuclear Information System (INIS)

Buet, C.

2005-01-01

My research work deals mainly with the mathematical modelling and the numerical simulation of plasma physics. This document is divided into 3 parts. The first one is a summary of the works done for the numerical solving of collision operators. The common thread of this part is obtaining numerical schemes preserving operators' properties namely physical invariants like mass, momentum and energy, equilibrium states and entropy decrease. These properties are generally checked formally for continuous operators, may give rise to some difficulties for discrete operators. In the second part I present a summary of the works regarding moments methods applied to radiative transfer and the numerical issues dealing with their discretization. The common thread of this part is how to get numerical schemes preserving asymptotic scattering and invariant domains for Lorentz models and also for non-linear telegraph-type equations involved in radiative transfer or electronic plasma. In the third part I present 2 themes linked to collision operators: multi-fluid ionization and the non-existence of linear monotone schemes for some linear parabolic equations

Fractional-Order Nonlinear Systems Modeling, Analysis and Simulation

CERN Document Server

Petráš, Ivo

2011-01-01

"Fractional-Order Nonlinear Systems: Modeling, Analysis and Simulation" presents a study of fractional-order chaotic systems accompanied by Matlab programs for simulating their state space trajectories, which are shown in the illustrations in the book. Description of the chaotic systems is clearly presented and their analysis and numerical solution are done in an easy-to-follow manner. Simulink models for the selected fractional-order systems are also presented. The readers will understand the fundamentals of the fractional calculus, how real dynamical systems can be described using fractional derivatives and fractional differential equations, how such equations can be solved, and how to simulate and explore chaotic systems of fractional order. The book addresses to mathematicians, physicists, engineers, and other scientists interested in chaos phenomena or in fractional-order systems. It can be used in courses on dynamical systems, control theory, and applied mathematics at graduate or postgraduate level. ...

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

Science.gov (United States)

Goorjian, Peter M.; Silberberg, Yaron; Kwak, Dochan (Technical Monitor)

1995-01-01

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Study of the nonlinear imperfect software debugging model

International Nuclear Information System (INIS)

Wang, Jinyong; Wu, Zhibo

2016-01-01

In recent years there has been a dramatic proliferation of research on imperfect software debugging phenomena. Software debugging is a complex process and is affected by a variety of factors, including the environment, resources, personnel skills, and personnel psychologies. Therefore, the simple assumption that debugging is perfect is inconsistent with the actual software debugging process, wherein a new fault can be introduced when removing a fault. Furthermore, the fault introduction process is nonlinear, and the cumulative number of nonlinearly introduced faults increases over time. Thus, this paper proposes a nonlinear, NHPP imperfect software debugging model in consideration of the fact that fault introduction is a nonlinear process. The fitting and predictive power of the NHPP-based proposed model are validated through related experiments. Experimental results show that this model displays better fitting and predicting performance than the traditional NHPP-based perfect and imperfect software debugging models. S-confidence bounds are set to analyze the performance of the proposed model. This study also examines and discusses optimal software release-time policy comprehensively. In addition, this research on the nonlinear process of fault introduction is significant given the recent surge of studies on software-intensive products, such as cloud computing and big data. - Highlights: • Fault introduction is a nonlinear changing process during the debugging phase. • The assumption that the process of fault introduction is nonlinear is credible. • Our proposed model can better fit and accurately predict software failure behavior. • Research on fault introduction case is significant to software-intensive products.

A Nonlinear differential equation model of Asthma effect of environmental pollution using LHAM

Science.gov (United States)

Joseph, G. Arul; Balamuralitharan, S.

2018-04-01

In this paper, we investigated a nonlinear differential equation mathematical model to study the spread of asthma in the environmental pollutants from industry and mainly from tobacco smoke from smokers in different type of population. Smoking is the main cause to spread Asthma in the environment. Numerical simulation is also discussed. Finally by using Liao’s Homotopy analysis Method (LHAM), we found that the approximate analytical solution of Asthmatic disease in the environmental.

Heterotic sigma models and non-linear strings

International Nuclear Information System (INIS)

Hull, C.M.

1986-01-01

The two-dimensional supersymmetric non-linear sigma models are examined with respect to the heterotic string. The paper was presented at the workshop on :Supersymmetry and its applications', Cambridge, United Kingdom, 1985. The non-linear sigma model with Wess-Zumino-type term, the coupling of the fermionic superfields to the sigma model, super-conformal invariance, and the supersymmetric string, are all discussed. (U.K.)

Nonlinear Multigrid for Reservoir Simulation

DEFF Research Database (Denmark)

Christensen, Max la Cour; Eskildsen, Klaus Langgren; Engsig-Karup, Allan Peter

2016-01-01

efficiency for a black-oil model. Furthermore, the use of the FAS method enables a significant reduction in memory usage compared with conventional techniques, which suggests new possibilities for improved large-scale reservoir simulation and numerical efficiency. Last, nonlinear multilevel preconditioning...

A NURBS-based finite element model applied to geometrically nonlinear elastodynamics using a corotational approach

KAUST Repository

Espath, L. F R

2015-02-03

A numerical model to deal with nonlinear elastodynamics involving large rotations within the framework of the finite element based on NURBS (Non-Uniform Rational B-Spline) basis is presented. A comprehensive kinematical description using a corotational approach and an orthogonal tensor given by the exact polar decomposition is adopted. The state equation is written in terms of corotational variables according to the hypoelastic theory, relating the Jaumann derivative of the Cauchy stress to the Eulerian strain rate.The generalized-α method (Gα) method and Generalized Energy-Momentum Method with an additional parameter (GEMM+ξ) are employed in order to obtain a stable and controllable dissipative time-stepping scheme with algorithmic conservative properties for nonlinear dynamic analyses.The main contribution is to show that the energy-momentum conservation properties and numerical stability may be improved once a NURBS-based FEM in the spatial discretization is used. Also it is shown that high continuity can postpone the numerical instability when GEMM+ξ with consistent mass is employed; likewise, increasing the continuity class yields a decrease in the numerical dissipation. A parametric study is carried out in order to show the stability and energy budget in terms of several properties such as continuity class, spectral radius and lumped as well as consistent mass matrices.

A NURBS-based finite element model applied to geometrically nonlinear elastodynamics using a corotational approach

KAUST Repository

Espath, L. F R; Braun, Alexandre Luis; Awruch, Armando Miguel; Dalcin, Lisandro

2015-01-01

A numerical model to deal with nonlinear elastodynamics involving large rotations within the framework of the finite element based on NURBS (Non-Uniform Rational B-Spline) basis is presented. A comprehensive kinematical description using a corotational approach and an orthogonal tensor given by the exact polar decomposition is adopted. The state equation is written in terms of corotational variables according to the hypoelastic theory, relating the Jaumann derivative of the Cauchy stress to the Eulerian strain rate.The generalized-α method (Gα) method and Generalized Energy-Momentum Method with an additional parameter (GEMM+ξ) are employed in order to obtain a stable and controllable dissipative time-stepping scheme with algorithmic conservative properties for nonlinear dynamic analyses.The main contribution is to show that the energy-momentum conservation properties and numerical stability may be improved once a NURBS-based FEM in the spatial discretization is used. Also it is shown that high continuity can postpone the numerical instability when GEMM+ξ with consistent mass is employed; likewise, increasing the continuity class yields a decrease in the numerical dissipation. A parametric study is carried out in order to show the stability and energy budget in terms of several properties such as continuity class, spectral radius and lumped as well as consistent mass matrices.

Numerical Modelling of Electrical Discharges

International Nuclear Information System (INIS)

Durán-Olivencia, F J; Pontiga, F; Castellanos, A

2014-01-01

The problem of the propagation of an electrical discharge between a spherical electrode and a plane has been solved by means of finite element methods (FEM) using a fluid approximation and assuming weak ionization and local equilibrium with the electric field. The numerical simulation of this type of problems presents the usual difficulties of convection-diffusion-reaction problems, in addition to those associated with the nonlinearities of the charged species velocities, the formation of steep gradients of the electric field and particle densities, and the coexistence of very different temporal scales. The effect of using different temporal discretizations for the numerical integration of the corresponding system of partial differential equations will be here investigated. In particular, the so-called θ-methods will be used, which allows to implement implicit, semi-explicit and fully explicit schemes in a simple way

Geometrically Nonlinear Shell Analysis of Wrinkled Thin-Film Membranes with Stress Concentrations

Science.gov (United States)

Tessler, Alexander; Sleight, David W.

2006-01-01

Geometrically nonlinear shell finite element analysis has recently been applied to solar-sail membrane problems in order to model the out-of-plane deformations due to structural wrinkling. Whereas certain problems lend themselves to achieving converged nonlinear solutions that compare favorably with experimental observations, solutions to tensioned membranes exhibiting high stress concentrations have been difficult to obtain even with the best nonlinear finite element codes and advanced shell element technology. In this paper, two numerical studies are presented that pave the way to improving the modeling of this class of nonlinear problems. The studies address the issues of mesh refinement and stress-concentration alleviation, and the effects of these modeling strategies on the ability to attain converged nonlinear deformations due to wrinkling. The numerical studies demonstrate that excessive mesh refinement in the regions of stress concentration may be disadvantageous to achieving wrinkled equilibrium states, causing the nonlinear solution to lock in the membrane response mode, while totally discarding the very low-energy bending response that is necessary to cause wrinkling deformation patterns.

The effect of nonlinear propagation on heating of tissue: A numerical model of diagnostic ultrasound beams

Science.gov (United States)

Cahill, Mark D.; Humphrey, Victor F.; Doody, Claire

2000-07-01

Thermal safety indices for diagnostic ultrasound beams are calculated under the assumption that the sound propagates under linear conditions. A non-axisymmetric finite difference model is used to solve the KZK equation, and so to model the beam of a diagnostic scanner in pulsed Doppler mode. Beams from both a uniform focused rectangular source and a linear array are considered. Calculations are performed in water, and in attenuating media with tissue-like characteristics. Attenuating media are found to exhibit significant nonlinear effects for finite-amplitude beams. The resulting loss of intensity by the beam is then used as the source term in a model of tissue heating to estimate the maximum temperature rises. These are compared with the thermal indices, derived from the properties of the water-propagated beams.

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.

Validation of a RANS transition model using a high-order weighted compact nonlinear scheme

Science.gov (United States)

Tu, GuoHua; Deng, XiaoGang; Mao, MeiLiang

2013-04-01

A modified transition model is given based on the shear stress transport (SST) turbulence model and an intermittency transport equation. The energy gradient term in the original model is replaced by flow strain rate to saving computational costs. The model employs local variables only, and then it can be conveniently implemented in modern computational fluid dynamics codes. The fifth-order weighted compact nonlinear scheme and the fourth-order staggered scheme are applied to discrete the governing equations for the purpose of minimizing discretization errors, so as to mitigate the confusion between numerical errors and transition model errors. The high-order package is compared with a second-order TVD method on simulating the transitional flow of a flat plate. Numerical results indicate that the high-order package give better grid convergence property than that of the second-order method. Validation of the transition model is performed for transitional flows ranging from low speed to hypersonic speed.

A novel numerical model for estimating the collapse pressure of flexible pipes

Energy Technology Data Exchange (ETDEWEB)

Nogueira, Victor P.P.; Antoun Netto, Theodoro [Universidade Federal do Rio de Janeiro (COPPE/UFRJ), RJ (Brazil). Coordenacao dos Programas de Pos-graduacao em Engenharia], e-mail: victor@lts.coppe.ufrj.br

2009-07-01

As the worldwide oil and gas industry operational environments move to ultra-deep waters, failure mechanisms in flexible pipes such as instability of the armor layers under compression and hydrostatic collapse are more likely to occur. Therefore, it is important to develop reliable numerical tools to reproduce the failure mechanisms that may occur in flexible pipes. This work presents a representative finite element model of flexible pipe capable to reproduce its pre and post-collapse behavior under hydrostatic pressure. The model, developed in the scope of this work, uses beam elements and includes nonlinear kinematics and material behavior influences. The dependability of the numerical results is assessed in light of experimental tests on flexible pipes with 4 inches and 8 inches nominal diameter available in the literature (Souza, 2002). The applied methodology provided coherent values regarding the estimation of the collapse pressures and results have shown that the proposed model is capable to reproduce experimental results. (author)

Numerical studies of non-linear evolution of kink and tearing modes in tokamaks

International Nuclear Information System (INIS)

White, R.; Monticello, D.; Rosenbluth, M.N.; Strauss, H.; Kadomtsev, B.B.

1975-01-01

A set of numerical techniques for investigating the full nonlinear unstable behavior of low β kink modes of given helical symmetry in Tokamaks is presented. Uniform current density plasmas display complicated deformations including the formation of large vacuum bubbles provided that the safety factor q is sufficiently close to integral. Fairly large m = 1 deformations, but not bubble formation, persist for a plasma with a parabolic current density profile (and hence shear). Deformations for m greater than or equal to 2 are, however, greatly suppressed. (auth)

Nonlinear distortion in wireless systems modeling and simulation with Matlab

CERN Document Server

Gharaibeh, Khaled M

2011-01-01

This book covers the principles of modeling and simulation of nonlinear distortion in wireless communication systems with MATLAB simulations and techniques In this book, the author describes the principles of modeling and simulation of nonlinear distortion in single and multichannel wireless communication systems using both deterministic and stochastic signals. Models and simulation methods of nonlinear amplifiers explain in detail how to analyze and evaluate the performance of data communication links under nonlinear amplification. The book addresses the analysis of nonlinear systems

Numerical Investigation of Corrugated Wire Mesh Laminate

Directory of Open Access Journals (Sweden)

Jeongho Choi

2013-01-01

Full Text Available The aim of this work is to develop a numerical model of Corrugated Wire Mesh Laminate (CWML capturing all its complexities such as nonlinear material properties, nonlinear geometry and large deformation behaviour, and frictional behaviour. Development of such a model will facilitate numerical simulation of the mechanical behaviour of the wire mesh structure under various types of loading as well as the variation of the CWML configuration parameters to tailor its mechanical properties to suit the intended application. Starting with a single strand truss model consisting of four waves with a bilinear stress-strain model to represent the plastic behaviour of stainless steel, the finite element model is gradually built up to study single-layer structures with 18 strands of corrugated wire meshes consistency and double- and quadruple-layered laminates with alternating crossply orientations. The compressive behaviour of the CWML model is simulated using contact elements to model friction and is compared to the load-deflection behaviour determined experimentally in uniaxial compression tests. The numerical model of the CWML is then employed to conduct the aim of establishing the upper and lower bounds of stiffness and load capacity achievable by such structures.

Nonlinear electron transport in magnetized laser plasmas

International Nuclear Information System (INIS)

Kho, T.H.; Haines, M.G.

1986-01-01

Electron transport in a magnetized plasma heated by inverse bremsstrahlung is studied numerically using a nonlinear Fokker--Planck model with self-consistent E and B fields. The numerical scheme is described. Nonlocal transport is found to alter many of the transport coefficients derived from linear transport theory, in particular, the Nernst and Righi--Leduc effects, in addition to the perpendicular heat flux q/sub perpendicular/, are substantially reduced near critical surface. The magnetic field, however, remains strongly coupled to the nonlinear q/sub perpendicular/ and, as has been found in hydrosimulations, convective amplification of the magnetic field occurs in the overdense plasma

Mathematical and numerical methods for the nonlinear hyperbolic propagation problem: d2GAMMA/dt2 = d/dz [dGAMMA/dt dGAMMA/dz

International Nuclear Information System (INIS)

Stanco, L.; Vaccaro, V.G.; Funk, U.; Krueger, U.; Mika, K.; Wuestefeld, G.

1982-03-01

In the first part of this report a physical model is presented, which describes the deforming of a bunch in a storage ring influenced only by its own space charge field. A system of two differential equations for the density and the momentum of the particles is set up, which is independent of any special machine parameter. Due to the sign of the inductance of the chamber walls and the sign of the dispersion of the revolution frequency, we distinguish between a de-bunching and a self-bunching situation. The de-bunching corresponds to a nonlinear hyperbolic propagation problem well-known in gas dynamics, and the self-bunching to a nonlinear elliptic initial value problem. The second part deals with a mathematical and numerical treatment of an approximate equation for the hyperbolic case. For this nonlinear second order partial differential equation we first present three particular integrals: the solution by separating the variables, the similarity solution, and the solution for a parabolic initial distribution of the density. For a more realistic initial condition, we must resort to other methods: Results are obtained in three different ways, first from a highly accurate Taylor series expansion, second from a common finite difference method, and thirdly from the numerical method of characteristics. The appearance of a shock discontinuity is furthermore established in each of these cases. (orig.)

Modelling nonlinear viscoelastic behaviours of loudspeaker suspensions-like structures

Science.gov (United States)

Maillou, Balbine; Lotton, Pierrick; Novak, Antonin; Simon, Laurent

2018-03-01

Mechanical properties of an electrodynamic loudspeaker are mainly determined by its suspensions (surround and spider) that behave nonlinearly and typically exhibit frequency dependent viscoelastic properties such as creep effect. The paper aims at characterizing the mechanical behaviour of electrodynamic loudspeaker suspensions at low frequencies using nonlinear identification techniques developed in recent years. A Generalized Hammerstein based model can take into account both frequency dependency and nonlinear properties. As shown in the paper, the model generalizes existing nonlinear or viscoelastic models commonly used for loudspeaker modelling. It is further experimentally shown that a possible input-dependent law may play a key role in suspension characterization.

Collapse arresting in an inhomogeneous quintic nonlinear Schrodinger model

DEFF Research Database (Denmark)

Gaididei, Yuri Borisovich; Schjødt-Eriksen, Jens; Christiansen, Peter Leth

1999-01-01

Collapse of (1 + 1)-dimensional beams in the inhomogeneous one-dimensional quintic nonlinear Schrodinger equation is analyzed both numerically and analytically. It is shown that in the vicinity of a narrow attractive inhomogeneity, the collapse of beams in which the homogeneous medium would blow up...

Estimation of biological parameters of marine organisms using linear and nonlinear acoustic scattering model-based inversion methods.

Science.gov (United States)

Chu, Dezhang; Lawson, Gareth L; Wiebe, Peter H

2016-05-01

The linear inversion commonly used in fisheries and zooplankton acoustics assumes a constant inversion kernel and ignores the uncertainties associated with the shape and behavior of the scattering targets, as well as other relevant animal parameters. Here, errors of the linear inversion due to uncertainty associated with the inversion kernel are quantified. A scattering model-based nonlinear inversion method is presented that takes into account the nonlinearity of the inverse problem and is able to estimate simultaneously animal abundance and the parameters associated with the scattering model inherent to the kernel. It uses sophisticated scattering models to estimate first, the abundance, and second, the relevant shape and behavioral parameters of the target organisms. Numerical simulations demonstrate that the abundance, size, and behavior (tilt angle) parameters of marine animals (fish or zooplankton) can be accurately inferred from the inversion by using multi-frequency acoustic data. The influence of the singularity and uncertainty in the inversion kernel on the inversion results can be mitigated by examining the singular values for linear inverse problems and employing a non-linear inversion involving a scattering model-based kernel.

A Novel Machine Learning Strategy Based on Two-Dimensional Numerical Models in Financial Engineering

Directory of Open Access Journals (Sweden)

Qingzhen Xu

2013-01-01

Full Text Available Machine learning is the most commonly used technique to address larger and more complex tasks by analyzing the most relevant information already present in databases. In order to better predict the future trend of the index, this paper proposes a two-dimensional numerical model for machine learning to simulate major U.S. stock market index and uses a nonlinear implicit finite-difference method to find numerical solutions of the two-dimensional simulation model. The proposed machine learning method uses partial differential equations to predict the stock market and can be extensively used to accelerate large-scale data processing on the history database. The experimental results show that the proposed algorithm reduces the prediction error and improves forecasting precision.

New non-linear model of groundwater recharge: Inclusion of memory, heterogeneity and visco-elasticity

Directory of Open Access Journals (Sweden)

Spannenberg Jescica

2017-09-01

Full Text Available Fractional differentiation has adequate use for investigating real world scenarios related to geological formations associated with elasticity, heterogeneity, viscoelasticity, and the memory effect. Since groundwater systems exist in these geological formations, modelling groundwater recharge as a real world scenario is a challenging task to do because existing recharge estimation methods are governed by linear equations which make use of constant field parameters. This is inadequate because in reality these parameters are a function of both space and time. This study therefore concentrates on modifying the recharge equation governing the EARTH model, by application of the Eton approach. Accordingly, this paper presents a modified equation which is non-linear, and accounts for parameters in a way that it is a function of both space and time. To be more specific, herein, recharge and drainage resistance which are parameters within the equation, became a function of both space and time. Additionally, the study entailed solving the non-linear equation using an iterative method as well as numerical solutions by means of the Crank-Nicolson scheme. The numerical solutions were used alongside the Riemann-Liouville, Caputo-Fabrizio, and Atangana-Baleanu derivatives, so that account was taken for elasticity, heterogeneity, viscoelasticity, and the memory effect. In essence, this paper presents a more adequate model for recharge estimation.

Analysis of nonlinear parabolic equations modeling plasma diffusion across a magnetic field

International Nuclear Information System (INIS)

Hyman, J.M.; Rosenau, P.

1984-01-01

We analyse the evolutionary behavior of the solution of a pair of coupled quasilinear parabolic equations modeling the diffusion of heat and mass of a magnetically confined plasma. The solutions's behavior, due to the nonlinear diffusion coefficients, exhibits many new phenomena. In short time, the solution converges into a highly organized symmetric pattern that is almost completely independent of initial data. The asymptotic dynamics then become very simple and take place in a finite dimensional space. These conclusions are backed by extensive numerical experimentation

Non-Linear Slosh Damping Model Development and Validation

Science.gov (United States)

Yang, H. Q.; West, Jeff

2015-01-01

Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can

Efficient uncertainty quantification of a fully nonlinear and dispersive water wave model with random inputs

DEFF Research Database (Denmark)

Bigoni, Daniele; Engsig-Karup, Allan Peter; Eskilsson, Claes

2016-01-01

A major challenge in next-generation industrial applications is to improve numerical analysis by quantifying uncertainties in predictions. In this work we present a formulation of a fully nonlinear and dispersive potential flow water wave model with random inputs for the probabilistic description...... at different points in the parameter space, allowing for the reuse of existing simulation software. The choice of the applied methods is driven by the number of uncertain input parameters and by the fact that finding the solution of the considered model is computationally intensive. We revisit experimental...... benchmarks often used for validation of deterministic water wave models. Based on numerical experiments and assumed uncertainties in boundary data, our analysis reveals that some of the known discrepancies from deterministic simulation in comparison with experimental measurements could be partially explained...

Model reduction of nonlinear systems subject to input disturbances

KAUST Repository

Ndoye, Ibrahima

2017-07-10

The method of convex optimization is used as a tool for model reduction of a class of nonlinear systems in the presence of disturbances. It is shown that under some conditions the nonlinear disturbed system can be approximated by a reduced order nonlinear system with similar disturbance-output properties to the original plant. The proposed model reduction strategy preserves the nonlinearity and the input disturbance nature of the model. It guarantees a sufficiently small error between the outputs of the original and the reduced-order systems, and also maintains the properties of input-to-state stability. The matrices of the reduced order system are given in terms of a set of linear matrix inequalities (LMIs). The paper concludes with a demonstration of the proposed approach on model reduction of a nonlinear electronic circuit with additive disturbances.

A NONLINEAR MATHEMATICAL MODEL FOR ASTHMA: EFFECT OF ENVIRONMENTAL POLLUTION

Directory of Open Access Journals (Sweden)

NARESHA RAM

2009-04-01

Full Text Available In this paper, we explore a nonlinear mathematical model to study the spread of asthma due to inhaled pollutants from industry as well as tobacco smoke from smokers in a variable size population. The model is analyzed using stability theory of differential equations and computer simulation. It is shown that with an increase in the level of air pollutants concentration, the asthmatic (diseased population increases. It is also shown that along with pollutants present in the environment, smoking (active or passive also helps in the spread of asthma. Moreover, with the increase in the rate of interaction between susceptibles and smokers, the persistence of the spread of asthma is higher. A numerical study of the model is also performed to see the role of certain key parameters on the spread of asthma and to support the analytical results.

Non-linear Growth Models in Mplus and SAS

Science.gov (United States)

Grimm, Kevin J.; Ram, Nilam

2013-01-01

Non-linear growth curves or growth curves that follow a specified non-linear function in time enable researchers to model complex developmental patterns with parameters that are easily interpretable. In this paper we describe how a variety of sigmoid curves can be fit using the Mplus structural modeling program and the non-linear mixed-effects modeling procedure NLMIXED in SAS. Using longitudinal achievement data collected as part of a study examining the effects of preschool instruction on academic gain we illustrate the procedures for fitting growth models of logistic, Gompertz, and Richards functions. Brief notes regarding the practical benefits, limitations, and choices faced in the fitting and estimation of such models are included. PMID:23882134

Finite difference modelling of the temperature rise in non-linear medical ultrasound fields.

Science.gov (United States)

Divall, S A; Humphrey, V F

2000-03-01

Non-linear propagation of ultrasound can lead to increased heat generation in medical diagnostic imaging due to the preferential absorption of harmonics of the original frequency. A numerical model has been developed and tested that is capable of predicting the temperature rise due to a high amplitude ultrasound field. The acoustic field is modelled using a numerical solution to the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, known as the Bergen Code, which is implemented in cylindrical symmetric form. A finite difference representation of the thermal equations is used to calculate the resulting temperature rises. The model allows for the inclusion of a number of layers of tissue with different acoustic and thermal properties and accounts for the effects of non-linear propagation, direct heating by the transducer, thermal diffusion and perfusion in different tissues. The effect of temperature-dependent skin perfusion and variation in background temperature between the skin and deeper layers of the body are included. The model has been tested against analytic solutions for simple configurations and then used to estimate temperature rises in realistic obstetric situations. A pulsed 3 MHz transducer operating with an average acoustic power of 200 mW leads to a maximum steady state temperature rise inside the foetus of 1.25 degrees C compared with a 0.6 degree C rise for the same transmitted power under linear propagation conditions. The largest temperature rise occurs at the skin surface, with the temperature rise at the foetus limited to less than 2 degrees C for the range of conditions considered.

Nonlinear mirror mode dynamics: Simulations and modeling

Czech Academy of Sciences Publication Activity Database

Califano, F.; Hellinger, Petr; Kuznetsov, E.; Passot, T.; Sulem, P. L.; Trávníček, Pavel

2008-01-01

Roč. 113, - (2008), A08219/1-A08219/20 ISSN 0148-0227 R&D Projects: GA AV ČR IAA300420702; GA AV ČR IAA300420602 Grant - others:PECS(CZ) 98024 Institutional research plan: CEZ:AV0Z30420517 Keywords : mirror instability * nonlinear evolution * numerical simulations * magnetic holes * mirror structures * kinetic plasma instabilities Subject RIV: BL - Plasma and Gas Discharge Physics Impact factor: 3.147, year: 2008

Model of anisotropic nonlinearity in self-defocusing photorefractive media.

Science.gov (United States)

Barsi, C; Fleischer, J W

2015-09-21

We develop a phenomenological model of anisotropy in self-defocusing photorefractive crystals. In addition to an independent term due to nonlinear susceptibility, we introduce a nonlinear, non-separable correction to the spectral diffraction operator. The model successfully describes the crossover between photovoltaic and photorefractive responses and the spatially dispersive shock wave behavior of a nonlinearly spreading Gaussian input beam. It should prove useful for characterizing internal charge dynamics in complex materials and for accurate image reconstruction through nonlinear media.

Analytical model for nonlinear piezoelectric energy harvesting devices

International Nuclear Information System (INIS)

Neiss, S; Goldschmidtboeing, F; M Kroener; Woias, P

2014-01-01

In this work we propose analytical expressions for the jump-up and jump-down point of a nonlinear piezoelectric energy harvester. In addition, analytical expressions for the maximum power output at optimal resistive load and the 3 dB-bandwidth are derived. So far, only numerical models have been used to describe the physics of a piezoelectric energy harvester. However, this approach is not suitable to quickly evaluate different geometrical designs or piezoelectric materials in the harvester design process. In addition, the analytical expressions could be used to predict the jump-frequencies of a harvester during operation. In combination with a tuning mechanism, this would allow the design of an efficient control algorithm to ensure that the harvester is always working on the oscillator's high energy attractor. (paper)

Numerical study of surface plasmon enhanced nonlinear absorption and refraction.

Science.gov (United States)

Kohlgraf-Owens, Dana C; Kik, Pieter G

2008-07-07

Maxwell Garnett effective medium theory is used to study the influence of silver nanoparticle induced field enhancement on the nonlinear response of a Kerr-type nonlinear host. We show that the composite nonlinear absorption coefficient, beta(c), can be enhanced relative to the host nonlinear absorption coefficient near the surface plasmon resonance of silver nanoparticles. This enhancement is not due to a resonant enhancement of the host nonlinear absorption, but rather due to a phase shifted enhancement of the host nonlinear refractive response. The enhancement occurs at the expense of introducing linear absorption, alpha(c), which leads to an overall reduced figure of merit beta(c)/alpha(c) for nonlinear absorption. For thin (< 1 microm) composites, the use of surface plasmons is found to result in an increased nonlinear absorption response compared to that of the host material.

Nonlinear Control of Hydraulic Manipulator for Decommissioning Nuclear Reactor

Energy Technology Data Exchange (ETDEWEB)

Kim, Myoung-Ho; Lee, Sung-Uk; Kim, Chang-Hoi; Choi, Byung-Seon; Moon, Jei-Kwon [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2016-10-15

Robot technique is need to decommission nuclear reactor because of high radiation environment. Especially, Manipulator systems are useful for dismantling complex structure in a nuclear facility. In addition, Hydraulic system is applied to handle heavy duty object. Since hydraulic system can demonstrate high power. The manipulator with hydraulic power is already developed. To solve this problem, various nonlinear control method includes acceleration control. But, it is difficult because acceleration value is highly noisy. In this paper, the nonlinear control algorithm without acceleration control is studied. To verify, the hydraulic manipulator model had been developed. Furthermore, the numerical simulation is carried out. The nonlinear control without acceleration parameter method is developed for hydraulic manipulator. To verify control algorithm, the manipulator is modeled by MBD and the hydraulic servo system is also derived. In addition, the numerical simulation is also carried out. Especially, PID gain is determined though TDC algorithm. In the result of numerical simulation, tracking performance is good without acceleration control. Thus, the PID though TDC with SMC is good for hydraulic manipulator control.

Nonlinear Control of Hydraulic Manipulator for Decommissioning Nuclear Reactor

International Nuclear Information System (INIS)

Kim, Myoung-Ho; Lee, Sung-Uk; Kim, Chang-Hoi; Choi, Byung-Seon; Moon, Jei-Kwon

2016-01-01

Robot technique is need to decommission nuclear reactor because of high radiation environment. Especially, Manipulator systems are useful for dismantling complex structure in a nuclear facility. In addition, Hydraulic system is applied to handle heavy duty object. Since hydraulic system can demonstrate high power. The manipulator with hydraulic power is already developed. To solve this problem, various nonlinear control method includes acceleration control. But, it is difficult because acceleration value is highly noisy. In this paper, the nonlinear control algorithm without acceleration control is studied. To verify, the hydraulic manipulator model had been developed. Furthermore, the numerical simulation is carried out. The nonlinear control without acceleration parameter method is developed for hydraulic manipulator. To verify control algorithm, the manipulator is modeled by MBD and the hydraulic servo system is also derived. In addition, the numerical simulation is also carried out. Especially, PID gain is determined though TDC algorithm. In the result of numerical simulation, tracking performance is good without acceleration control. Thus, the PID though TDC with SMC is good for hydraulic manipulator control

Non-Linear Numerical Modeling and Experimental Testing of a Point Absorber Wave Energy Converter

DEFF Research Database (Denmark)

Zurkinden, Andrew Stephen; Ferri, Francesco; Beatty, S.

2014-01-01

the calculation of the non-linear hydrostatic restoring moment by a cubic polynomial function fit to laboratory test results. Moreover, moments due to viscous drag are evaluated on the oscillating hemisphere considering the horizontal and vertical drag force components. The influence on the motions of this non.......e. H/λ≤0.02. For steep waves, H/λ≥0.04 however, the relative velocities between the body and the waves increase thus requiring inclusion of the non-linear hydrostatic restoring moment to effectively predict the dynamics of the wave energy converter. For operation of the device with a passively damping...

Nonlinear model of epidemic spreading in a complex social network.

Science.gov (United States)

Kosiński, Robert A; Grabowski, A

2007-10-01

The epidemic spreading in a human society is a complex process, which can be described on the basis of a nonlinear mathematical model. In such an approach the complex and hierarchical structure of social network (which has implications for the spreading of pathogens and can be treated as a complex network), can be taken into account. In our model each individual has one of the four permitted states: susceptible, infected, infective, unsusceptible or dead. This refers to the SEIR model used in epidemiology. The state of an individual changes in time, depending on the previous state and the interactions with other individuals. The description of the interpersonal contacts is based on the experimental observations of the social relations in the community. It includes spatial localization of the individuals and hierarchical structure of interpersonal interactions. Numerical simulations were performed for different types of epidemics, giving the progress of a spreading process and typical relationships (e.g. range of epidemic in time, the epidemic curve). The spreading process has a complex and spatially chaotic character. The time dependence of the number of infective individuals shows the nonlinear character of the spreading process. We investigate the influence of the preventive vaccinations on the spreading process. In particular, for a critical value of preventively vaccinated individuals the percolation threshold is observed and the epidemic is suppressed.

Nonlinear model predictive control theory and algorithms

CERN Document Server

Grüne, Lars

2017-01-01

This book offers readers a thorough and rigorous introduction to nonlinear model predictive control (NMPC) for discrete-time and sampled-data systems. NMPC schemes with and without stabilizing terminal constraints are detailed, and intuitive examples illustrate the performance of different NMPC variants. NMPC is interpreted as an approximation of infinite-horizon optimal control so that important properties like closed-loop stability, inverse optimality and suboptimality can be derived in a uniform manner. These results are complemented by discussions of feasibility and robustness. An introduction to nonlinear optimal control algorithms yields essential insights into how the nonlinear optimization routine—the core of any nonlinear model predictive controller—works. Accompanying software in MATLAB® and C++ (downloadable from extras.springer.com/), together with an explanatory appendix in the book itself, enables readers to perform computer experiments exploring the possibilities and limitations of NMPC. T...

Analytical modeling of soliton interactions in a nonlocal nonlinear medium analogous to gravitational force

Science.gov (United States)

Zeng, Shihao; Chen, Manna; Zhang, Ting; Hu, Wei; Guo, Qi; Lu, Daquan

2018-01-01

We illuminate an analytical model of soliton interactions in lead glass by analogizing to a gravitational force system. The orbits of spiraling solitons under a long-range interaction are given explicitly and demonstrated to follow Newton's second law of motion and the Binet equation by numerical simulations. The condition for circular orbits is obtained and the oscillating orbits are proved not to be closed. We prove the analogy between the nonlocal nonlinear optical system and gravitational system and specify the quantitative relation of the quantity between the two models.

Numerical Transducer Modeling

DEFF Research Database (Denmark)

Henriquez, Vicente Cutanda

This thesis describes the development of a numerical model of the propagation of sound waves in fluids with viscous and thermal losses, with application to the simulation of acoustic transducers, in particular condenser microphones for measurement. The theoretical basis is presented, numerical...... manipulations are developed to satisfy the more complicated boundary conditions, and a model of a condenser microphone with a coupled membrane is developed. The model is tested against measurements of ¼ inch condenser microphones and analytical calculations. A detailed discussion of the results is given....

Nonlinear analysis of flexible plates lying on elastic foundation

Directory of Open Access Journals (Sweden)

Trushin Sergey

2017-01-01

Full Text Available This article describes numerical procedures for analysis of flexible rectangular plates lying on elastic foundation. Computing models are based on the theory of plates with account of transverse shear deformations. The finite difference energy method of discretization is used for reducing the initial continuum problem to finite dimensional problem. Solution procedures for nonlinear problem are based on Newton-Raphson method. This theory of plates and numerical methods have been used for investigation of nonlinear behavior of flexible plates on elastic foundation with different properties.

Numerical assessment of factors affecting nonlinear internal waves in the South China Sea

Science.gov (United States)

Li, Qiang

2014-02-01

Nonlinear internal waves in the South China Sea exhibit diverse characteristics, which are associated with the complex conditions in Luzon Strait, such as the double ridge topography, the Earth’s rotation, variations in stratification and the background current induced by the Kuroshio. These effects are individually assessed using the MITgcm. The performance of the model is first validated through comparison with field observations. Because of in-phased ray interaction, the western ridge in Luzon Strait intensifies the semidiurnal internal tides generated from the eastern ridge, thus reinforcing the formation of nonlinear internal waves. However, the ray interaction for K1 forcing becomes anti-phased so that the K1 internal tide generation is reduced by the western ridge. Not only does the rotational dispersion suppress internal tide generation, it also inhibits nonlinear steepening and consequent internal solitary wave formation. As a joint effect, the double ridges and the rotational dispersion result in a paradoxical phenomenon: diurnal barotropic tidal forcing is dominant in Luzon Strait, but semidiurnal internal tides prevail in the deep basin of the South China Sea. The seasonal variation of the Kuroshio is consistent with the seasonal appearance of nonlinear internal waves in the South China Sea. The model results show that the westward inflow due to the Kuroshio intrusion reduces the amplitude of internal tides in the South China Sea, causing the weakening or absence of internal solitary waves. Winter stratification cannot account for the significant reduction of nonlinear internal waves, because the amplitude growth of internal tides due to increased thermocline tilting counteracts the reduced nonlinearity caused by thermocline deepening.

Parameter Estimation of Nonlinear Models in Forestry.

OpenAIRE

Fekedulegn, Desta; Mac Siúrtáin, Máirtín Pádraig; Colbert, Jim J.

1999-01-01

Partial derivatives of the negative exponential, monomolecular, Mitcherlich, Gompertz, logistic, Chapman-Richards, von Bertalanffy, Weibull and the Richard’s nonlinear growth models are presented. The application of these partial derivatives in estimating the model parameters is illustrated. The parameters are estimated using the Marquardt iterative method of nonlinear regression relating top height to age of Norway spruce (Picea abies L.) from the Bowmont Norway Spruce Thinnin...

Nonlinearities in modified gravity cosmology: Signatures of modified gravity in the nonlinear matter power spectrum

International Nuclear Information System (INIS)

Cui Weiguang; Zhang Pengjie; Yang Xiaohu

2010-01-01

A large fraction of cosmological information on dark energy and gravity is encoded in the nonlinear regime. Precision cosmology thus requires precision modeling of nonlinearities in general dark energy and modified gravity models. We modify the Gadget-2 code and run a series of N-body simulations on modified gravity cosmology to study the nonlinearities. The modified gravity model that we investigate in the present paper is characterized by a single parameter ζ, which determines the enhancement of particle acceleration with respect to general relativity (GR), given the identical mass distribution (ζ=1 in GR). The first nonlinear statistics we investigate is the nonlinear matter power spectrum at k < or approx. 3h/Mpc, which is the relevant range for robust weak lensing power spectrum modeling at l < or approx. 2000. In this study, we focus on the relative difference in the nonlinear power spectra at corresponding redshifts where different gravity models have the same linear power spectra. This particular statistics highlights the imprint of modified gravity in the nonlinear regime and the importance of including the nonlinear regime in testing GR. By design, it is less susceptible to the sample variance and numerical artifacts. We adopt a mass assignment method based on wavelet to improve the power spectrum measurement. We run a series of tests to determine the suitable simulation specifications (particle number, box size, and initial redshift). We find that, the nonlinear power spectra can differ by ∼30% for 10% deviation from GR (|ζ-1|=0.1) where the rms density fluctuations reach 10. This large difference, on one hand, shows the richness of information on gravity in the corresponding scales, and on the other hand, invalidates simple extrapolations of some existing fitting formulae to modified gravity cosmology.

Multiplicative noise removal through fractional order tv-based model and fast numerical schemes for its approximation

Science.gov (United States)

Ullah, Asmat; Chen, Wen; Khan, Mushtaq Ahmad

2017-07-01

This paper introduces a fractional order total variation (FOTV) based model with three different weights in the fractional order derivative definition for multiplicative noise removal purpose. The fractional-order Euler Lagrange equation which is a highly non-linear partial differential equation (PDE) is obtained by the minimization of the energy functional for image restoration. Two numerical schemes namely an iterative scheme based on the dual theory and majorization- minimization algorithm (MMA) are used. To improve the restoration results, we opt for an adaptive parameter selection procedure for the proposed model by applying the trial and error method. We report numerical simulations which show the validity and state of the art performance of the fractional-order model in visual improvement as well as an increase in the peak signal to noise ratio comparing to corresponding methods. Numerical experiments also demonstrate that MMAbased methodology is slightly better than that of an iterative scheme.

Nonlinear integral equations for the sausage model

Science.gov (United States)

Ahn, Changrim; Balog, Janos; Ravanini, Francesco

2017-08-01

The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the O(3) sigma model preserving integrability. The target space is deformed from the sphere to ‘sausage’ shape by a deformation parameter ν. This model is defined by a factorizable S-matrix which is obtained by deforming that of the O(3) sigma model by a parameter λ. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of 1/λ integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonlinear integral equations (NLIEs), which are applicable to generic value of λ. Our derivation is based on T-Q relations extracted from the truncated TBA equations. For a consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between ν and λ. Dedicated to the memory of Petr Petrovich Kulish.

Towards high fidelity numerical wave tanks for modelling coastal and ocean engineering processes

Science.gov (United States)

Cozzuto, G.; Dimakopoulos, A.; de Lataillade, T.; Kees, C. E.

2017-12-01

With the increasing availability of computational resources, the engineering and research community is gradually moving towards using high fidelity Comutational Fluid Mechanics (CFD) models to perform numerical tests for improving the understanding of physical processes pertaining to wave propapagation and interaction with the coastal environment and morphology, either physical or man-made. It is therefore important to be able to reproduce in these models the conditions that drive these processes. So far, in CFD models the norm is to use regular (linear or nonlinear) waves for performing numerical tests, however, only random waves exist in nature. In this work, we will initially present the verification and validation of numerical wave tanks based on Proteus, an open-soruce computational toolkit based on finite element analysis, with respect to the generation, propagation and absorption of random sea states comprising of long non-repeating wave sequences. Statistical and spectral processing of results demonstrate that the methodologies employed (including relaxation zone methods and moving wave paddles) are capable of producing results of similar quality to the wave tanks used in laboratories (Figure 1). Subsequently cases studies of modelling complex process relevant to coastal defences and floating structures such as sliding and overturning of composite breakwaters, heave and roll response of floating caissons are presented. Figure 1: Wave spectra in the numerical wave tank (coloured symbols), compared against the JONSWAP distribution

A Simple Model for Nonlinear Confocal Ultrasonic Beams

Science.gov (United States)

Zhang, Dong; Zhou, Lin; Si, Li-Sheng; Gong, Xiu-Fen

2007-01-01

A confocally and coaxially arranged pair of focused transmitter and receiver represents one of the best geometries for medical ultrasonic imaging and non-invasive detection. We develop a simple theoretical model for describing the nonlinear propagation of a confocal ultrasonic beam in biological tissues. On the basis of the parabolic approximation and quasi-linear approximation, the nonlinear Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation is solved by using the angular spectrum approach. Gaussian superposition technique is applied to simplify the solution, and an analytical solution for the second harmonics in the confocal ultrasonic beam is presented. Measurements are performed to examine the validity of the theoretical model. This model provides a preliminary model for acoustic nonlinear microscopy.

An Effective Numerical Method and Its Utilization to Solution of Fractional Models Used in Bioengineering Applications

Directory of Open Access Journals (Sweden)

Petráš Ivo

2011-01-01

Full Text Available This paper deals with the fractional-order linear and nonlinear models used in bioengineering applications and an effective method for their numerical solution. The proposed method is based on the power series expansion of a generating function. Numerical solution is in the form of the difference equation, which can be simply applied in the Matlab/Simulink to simulate the dynamics of system. Several illustrative examples are presented, which can be widely used in bioengineering as well as in the other disciplines, where the fractional calculus is often used.

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.

Bayesian nonlinear structural FE model and seismic input identification for damage assessment of civil structures

Science.gov (United States)

Astroza, Rodrigo; Ebrahimian, Hamed; Li, Yong; Conte, Joel P.

2017-09-01

A methodology is proposed to update mechanics-based nonlinear finite element (FE) models of civil structures subjected to unknown input excitation. The approach allows to jointly estimate unknown time-invariant model parameters of a nonlinear FE model of the structure and the unknown time histories of input excitations using spatially-sparse output response measurements recorded during an earthquake event. The unscented Kalman filter, which circumvents the computation of FE response sensitivities with respect to the unknown model parameters and unknown input excitations by using a deterministic sampling approach, is employed as the estimation tool. The use of measurement data obtained from arrays of heterogeneous sensors, including accelerometers, displacement sensors, and strain gauges is investigated. Based on the estimated FE model parameters and input excitations, the updated nonlinear FE model can be interrogated to detect, localize, classify, and assess damage in the structure. Numerically simulated response data of a three-dimensional 4-story 2-by-1 bay steel frame structure with six unknown model parameters subjected to unknown bi-directional horizontal seismic excitation, and a three-dimensional 5-story 2-by-1 bay reinforced concrete frame structure with nine unknown model parameters subjected to unknown bi-directional horizontal seismic excitation are used to illustrate and validate the proposed methodology. The results of the validation studies show the excellent performance and robustness of the proposed algorithm to jointly estimate unknown FE model parameters and unknown input excitations.

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

Numerical reconstruction of wave field spatial distributions at the output and input planes of nonlinear medium with use of digital holography

International Nuclear Information System (INIS)

Nalegaev, S S; Petrov, N V; Bespalov, V G

2014-01-01

A numerical reconstruction of spatial distributions of optical radiation propagating through a volume of nonlinear medium at input and output planes of the medium was demonstrated using a scheme of digital holography. A nonlinear Schrodinger equation with Fourier Split-Step method was used as a tool to propagate wavefront in the volume of the medium. Time dependence of the refractive index change was not taken into account.

Numerical modeling for the retrofit of the hydraulic cooling subsystems in operating power plant

Science.gov (United States)

AlSaqoor, S.; Alahmer, A.; Al Quran, F.; Andruszkiewicz, A.; Kubas, K.; Regucki, P.; Wędrychowicz, W.

2017-08-01

This paper presents the possibility of using the numerical methods to analyze the work of hydraulic systems on the example of a cooling system of a power boiler auxiliary devices. The variety of conditions at which hydraulic system that operated in specific engineering subsystems requires an individualized approach to the model solutions that have been developed for these systems modernizing. A mathematical model of a series-parallel propagation for the cooling water was derived and iterative methods were used to solve the system of nonlinear equations. The results of numerical calculations made it possible to analyze different variants of a modernization of the studied system and to indicate its critical elements. An economic analysis of different options allows an investor to choose an optimal variant of a reconstruction of the installation.

Nonlinear Kalman Filtering in Affine Term Structure Models

DEFF Research Database (Denmark)

Christoffersen, Peter; Dorion, Christian; Jacobs, Kris

When the relationship between security prices and state variables in dynamic term structure models is nonlinear, existing studies usually linearize this relationship because nonlinear fi…ltering is computationally demanding. We conduct an extensive investigation of this linearization and analyze...... the potential of the unscented Kalman …filter to properly capture nonlinearities. To illustrate the advantages of the unscented Kalman …filter, we analyze the cross section of swap rates, which are relatively simple non-linear instruments, and cap prices, which are highly nonlinear in the states. An extensive...

A nonlinear beam model to describe the postbuckling of wide neo-Hookean beams

Science.gov (United States)

Lubbers, Luuk A.; van Hecke, Martin; Coulais, Corentin

2017-09-01

Wide beams can exhibit subcritical buckling, i.e. the slope of the force-displacement curve can become negative in the postbuckling regime. In this paper, we capture this intriguing behaviour by constructing a 1D nonlinear beam model, where the central ingredient is the nonlinearity in the stress-strain relation of the beams constitutive material. First, we present experimental and numerical evidence of a transition to subcritical buckling for wide neo-Hookean hyperelastic beams, when their width-to-length ratio exceeds a critical value of 12%. Second, we construct an effective 1D energy density by combining the Mindlin-Reissner kinematics with a nonlinearity in the stress-strain relation. Finally, we establish and solve the governing beam equations to analytically determine the slope of the force-displacement curve in the postbuckling regime. We find, without any adjustable parameters, excellent agreement between the 1D theory, experiments and simulations. Our work extends the understanding of the postbuckling of structures made of wide elastic beams and opens up avenues for the reverse-engineering of instabilities in soft and metamaterials.

On the instability of a 3-dimensional attachment line boundary layer: Weakly nonlinear theory and a numerical approach

Science.gov (United States)

Hall, P.; Malik, M. R.

1984-01-01

The instability of a three dimensional attachment line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time dependent Navier-Stokes equations for the attachment line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment line boundary layer is also investigated.

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

Science.gov (United States)

Hall, P.; Malik, M. R.

1986-01-01

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Nonlinear Growth Models in M"plus" and SAS

Science.gov (United States)

Grimm, Kevin J.; Ram, Nilam

2009-01-01

Nonlinear growth curves or growth curves that follow a specified nonlinear function in time enable researchers to model complex developmental patterns with parameters that are easily interpretable. In this article we describe how a variety of sigmoid curves can be fit using the M"plus" structural modeling program and the nonlinear…

Nonlinear response and bistability of driven ion acoustic waves

Science.gov (United States)

Akbari-Moghanjoughi, M.

2017-08-01

The hydrodynamic model is used to obtain a generalized pseudoforce equation through which the nonlinear response of periodically driven ion acoustic waves is studied in an electron-ion plasma with isothermal and adiabatic ion fluids. The pseudotime series, corresponding to different driving frequencies, indicates that nonlinearity effects appear more strongly for smaller frequency values. The existence of extra harmonic resonances in the nonlinear amplitude spectrum is a clear indication of the interaction of an external force with harmonic components of the nonlinear ion acoustic waves. It is shown that many plasma parameters significantly and differently affect the nonlinear resonance spectrum of ion acoustic excitations. A heuristic but accurate model for the foldover effect is used which quite satisfactorily predicts the bistability of driven plasma oscillations. It is remarked that the characteristic resonance peak of isothermal ion plasma oscillations appears at lower frequencies but is stronger compared to that of adiabatic ions. Comparison of the exact numerical results for fully nonlinear and approximate (weakly nonlinear) models indicates that a weakly nonlinear model exaggerates the hysteresis and jump phenomenon for higher values of the external force amplitude.

Numerical Simulation of Coupled Nonlinear Schrödinger Equations Using the Generalized Differential Quadrature Method

International Nuclear Information System (INIS)

Mokhtari, R.; Toodar, A. Samadi; Chegini, N. G.

2011-01-01

We the extend application of the generalized differential quadrature method (GDQM) to solve some coupled nonlinear Schrödinger equations. The cosine-based GDQM is employed and the obtained system of ordinary differential equations is solved via the fourth order Runge—Kutta method. The numerical solutions coincide with the exact solutions in desired machine precision and invariant quantities are conserved sensibly. Some comparisons with the methods applied in the literature are carried out. (general)

Nonlinear Dynamic Models in Advanced Life Support

Science.gov (United States)

Jones, Harry

2002-01-01

To facilitate analysis, ALS systems are often assumed to be linear and time invariant, but they usually have important nonlinear and dynamic aspects. Nonlinear dynamic behavior can be caused by time varying inputs, changes in system parameters, nonlinear system functions, closed loop feedback delays, and limits on buffer storage or processing rates. Dynamic models are usually cataloged according to the number of state variables. The simplest dynamic models are linear, using only integration, multiplication, addition, and subtraction of the state variables. A general linear model with only two state variables can produce all the possible dynamic behavior of linear systems with many state variables, including stability, oscillation, or exponential growth and decay. Linear systems can be described using mathematical analysis. Nonlinear dynamics can be fully explored only by computer simulations of models. Unexpected behavior is produced by simple models having only two or three state variables with simple mathematical relations between them. Closed loop feedback delays are a major source of system instability. Exceeding limits on buffer storage or processing rates forces systems to change operating mode. Different equilibrium points may be reached from different initial conditions. Instead of one stable equilibrium point, the system may have several equilibrium points, oscillate at different frequencies, or even behave chaotically, depending on the system inputs and initial conditions. The frequency spectrum of an output oscillation may contain harmonics and the sums and differences of input frequencies, but it may also contain a stable limit cycle oscillation not related to input frequencies. We must investigate the nonlinear dynamic aspects of advanced life support systems to understand and counter undesirable behavior.

Nonlinear Eddy Viscosity Models applied to Wind Turbine Wakes

DEFF Research Database (Denmark)

Laan, van der, Paul Maarten; Sørensen, Niels N.; Réthoré, Pierre-Elouan

2013-01-01

The linear k−ε eddy viscosity model and modified versions of two existing nonlinear eddy viscosity models are applied to single wind turbine wake simulations using a Reynolds Averaged Navier-Stokes code. Results are compared with field wake measurements. The nonlinear models give better results...

Nonlinear dynamic behavior of an assembly of tubes under transverse fluid flow

International Nuclear Information System (INIS)

Beaufils, B.; Axisa, F.; Antunes, J.

1989-01-01

The mechanical vibrations induced by a transverse fluid flow passing through an assembly of cylindrical tubes is investigated. Studies on the numerical modeling of such phenomena are presented. The purpose of the work is to allow the evaluation of the risks induced by the vibrations in industrial heat exchangers. The methods for the analysis of nonlinear problems and numerical calculations of the nonlinear dynamic behavior are performed [fr

Topological approximation of the nonlinear Anderson model

Science.gov (United States)

Milovanov, Alexander V.; Iomin, Alexander

2014-06-01

We study the phenomena of Anderson localization in the presence of nonlinear interaction on a lattice. A class of nonlinear Schrödinger models with arbitrary power nonlinearity is analyzed. We conceive the various regimes of behavior, depending on the topology of resonance overlap in phase space, ranging from a fully developed chaos involving Lévy flights to pseudochaotic dynamics at the onset of delocalization. It is demonstrated that the quadratic nonlinearity plays a dynamically very distinguished role in that it is the only type of power nonlinearity permitting an abrupt localization-delocalization transition with unlimited spreading already at the delocalization border. We describe this localization-delocalization transition as a percolation transition on the infinite Cayley tree (Bethe lattice). It is found in the vicinity of the criticality that the spreading of the wave field is subdiffusive in the limit t →+∞. The second moment of the associated probability distribution grows with time as a power law ∝ tα, with the exponent α =1/3 exactly. Also we find for superquadratic nonlinearity that the analog pseudochaotic regime at the edge of chaos is self-controlling in that it has feedback on the topology of the structure on which the transport processes concentrate. Then the system automatically (without tuning of parameters) develops its percolation point. We classify this type of behavior in terms of self-organized criticality dynamics in Hilbert space. For subquadratic nonlinearities, the behavior is shown to be sensitive to the details of definition of the nonlinear term. A transport model is proposed based on modified nonlinearity, using the idea of "stripes" propagating the wave process to large distances. Theoretical investigations, presented here, are the basis for consistency analysis of the different localization-delocalization patterns in systems with many coupled degrees of freedom in association with the asymptotic properties of the

A non-linear state space approach to model groundwater fluctuations

NARCIS (Netherlands)

Berendrecht, W.L.; Heemink, A.W.; Geer, F.C. van; Gehrels, J.C.

2006-01-01

A non-linear state space model is developed for describing groundwater fluctuations. Non-linearity is introduced by modeling the (unobserved) degree of water saturation of the root zone. The non-linear relations are based on physical concepts describing the dependence of both the actual

Augmented twin-nonlinear two-box behavioral models for multicarrier LTE power amplifiers.

Science.gov (United States)

Hammi, Oualid

2014-01-01

A novel class of behavioral models is proposed for LTE-driven Doherty power amplifiers with strong memory effects. The proposed models, labeled augmented twin-nonlinear two-box models, are built by cascading a highly nonlinear memoryless function with a mildly nonlinear memory polynomial with cross terms. Experimental validation on gallium nitride based Doherty power amplifiers illustrates the accuracy enhancement and complexity reduction achieved by the proposed models. When strong memory effects are observed, the augmented twin-nonlinear two-box models can improve the normalized mean square error by up to 3 dB for the same number of coefficients when compared to state-of-the-art twin-nonlinear two-box models. Furthermore, the augmented twin-nonlinear two-box models lead to the same performance as previously reported twin-nonlinear two-box models while requiring up to 80% less coefficients.

Multi-physics modeling and numerical simulation of weld pool in GTA welding

International Nuclear Information System (INIS)

Nguyen, Minh-Chien

2015-01-01

In this work, we develop a 3D physical and numerical model of the GTA (Gas Tungsten Arc) welding process in order to predict, for given welding parameters, useful quantities for the designer of welded assembly: weld bead shape, fluid flow in the weld pool as well as thermal distribution in the work piece. The model is developed in the Cast3M (http://www-cast3m.cea.fr/) finite element software and takes into account the main physical phenomena acting in the work piece and particularly in the weld pool, subject to source terms modeling the arc part of the welding process. A steady solution of this model is thought for and involves the coupling of the nonlinear thermohydraulics and electromagnetic equations together with the displacement of the deformable free surface of the weld pool. A first step in the development consisted in modeling the electromagnetic phenomena with two different numerical methods, in comparing the numerical results obtained with those of the literature and in quantifying the importance of the Lorentz force and the Joule effect compared to the other mechanical and thermal sources by computing power balances. Then, in order to assess the predictive capability of the model, simulations of various welding configurations are performed: variation in the chemical composition of the material, of the welding speed, of the prescribed arc pressure and of the welding positions, which is a focus of this work, are studied. A good agreement is obtained between the results of our model and other experimental and numerical results of the literature. Eventually, a model accounting for metal filling is proposed and its results are discussed. Thus, our complete model can be seen as a solid foundation towards future totally-coupled 3D welding models including the arc and it will be included in WPROCESS the in-house CEA software dedicated to the numerical simulation of welding. (author) [fr

Numerical Methods for a Multicomponent Two-Phase Interface Model with Geometric Mean Influence Parameters

KAUST Repository

Kou, Jisheng

2015-07-16

In this paper, we consider an interface model for multicomponent two-phase fluids with geometric mean influence parameters, which is popularly used to model and predict surface tension in practical applications. For this model, there are two major challenges in theoretical analysis and numerical simulation: the first one is that the influence parameter matrix is not positive definite; the second one is the complicated structure of the energy function, which requires us to find out a physically consistent treatment. To overcome these two challenging problems, we reduce the formulation of the energy function by employing a linear transformation and a weighted molar density, and furthermore, we propose a local minimum grand potential energy condition to establish the relation between the weighted molar density and mixture compositions. From this, we prove the existence of the solution under proper conditions and prove the maximum principle of the weighted molar density. For numerical simulation, we propose a modified Newton\\'s method for solving this nonlinear model and analyze its properties; we also analyze a finite element method with a physical-based adaptive mesh-refinement technique. Numerical examples are tested to verify the theoretical results and the efficiency of the proposed methods.

Numerical Methods for a Multicomponent Two-Phase Interface Model with Geometric Mean Influence Parameters

KAUST Repository

Kou, Jisheng; Sun, Shuyu

2015-01-01

In this paper, we consider an interface model for multicomponent two-phase fluids with geometric mean influence parameters, which is popularly used to model and predict surface tension in practical applications. For this model, there are two major challenges in theoretical analysis and numerical simulation: the first one is that the influence parameter matrix is not positive definite; the second one is the complicated structure of the energy function, which requires us to find out a physically consistent treatment. To overcome these two challenging problems, we reduce the formulation of the energy function by employing a linear transformation and a weighted molar density, and furthermore, we propose a local minimum grand potential energy condition to establish the relation between the weighted molar density and mixture compositions. From this, we prove the existence of the solution under proper conditions and prove the maximum principle of the weighted molar density. For numerical simulation, we propose a modified Newton's method for solving this nonlinear model and analyze its properties; we also analyze a finite element method with a physical-based adaptive mesh-refinement technique. Numerical examples are tested to verify the theoretical results and the efficiency of the proposed methods.

Dynamic nonlinear interaction of elastic plates on discrete supports

International Nuclear Information System (INIS)

Coutinho, A.L.G.A.; Landau, L.; Lima, E.C.P. de; Ebecken, N.F.F.

1984-01-01

A study on the dynamic nonlinear interaction of elastic plates using the finite element method is presented. The elastic plate is discretized by 4-node isoparametric Mindlin elements. The constitutive relation of the discrete supports can be any nonlinear curve given by pairs of force-displacement points. The nonlinear behaviour is represented by the overlay approach. This model also allows the simulation of a progressive decrease on the supports stiffnesses during load cycles. The dynamic nonlinear incremental movement equations are integrated by the Newmark implicit operator. Two alternatives for the incremental-iterative formulation are compared. The paper ends with a discussion of the advantages and limitations of the presented numerical models. (Author) [pt

A nonlinear model for AC induced corrosion

Directory of Open Access Journals (Sweden)

N. Ida

2012-09-01

Full Text Available The modeling of corrosion poses particular difficulties. The understanding of corrosion as an electrochemical process has led to simple capacitive-resistive models that take into account the resistance of the electrolytic cell and the capacitive effect of the surface potential at the interface between conductors and the electrolyte. In some models nonlinear conduction effects have been added to account for more complex observed behavior. While these models are sufficient to describe the behavior in systems with cathodic protection, the behavior in the presence of induced AC currents from power lines and from RF sources cannot be accounted for and are insufficient to describe the effects observed in the field. Field observations have shown that a rectifying effect exists that affects the cathodic protection potential and this effect is responsible for corrosion in the presence of AC currents. The rectifying effects of the metal-corrosion interface are totally missing from current models. This work proposes a nonlinear model based on finite element analysis that takes into account the nonlinear behavior of the metal-oxide interface and promises to improve modeling by including the rectification effects at the interface.

A Lie-admissible method of integration of Fokker-Planck equations with non-linear coefficients (exact and numerical solutions)

International Nuclear Information System (INIS)

Fronteau, J.; Combis, P.

1984-08-01

A Lagrangian method is introduced for the integration of non-linear Fokker-Planck equations. Examples of exact solutions obtained in this way are given, and also the explicit scheme used for the computation of numerical solutions. The method is, in addition, shown to be of a Lie-admissible type

Numerical Analysis of Ginzburg-Landau Models for Superconductivity.