E. Luckins, James M. Oliver, C. Please, Benjamin M. Sloman, A. Valderhaug, R. V. Van Gorder
Modelling the production of silicon in a submerged arc furnace (SAF) requires accounting for the wide range of timescales of the different physical and chemical processes: the electric current which is used to heat the furnace varies over a timescale of around $10^{-2},$ s, whereas the flow and chemical consumption of the raw materials in the furnace occurs over several hours. Models for the silicon furnace generally either include only the fast-timescale, or only the slow-timescale processes. In a prior work, we developed a model incorporating effects on both the fast and slow timescales, and used a multiple-timescales analysis to homogenise the fast variations, deriving an averaged model for the slow evolution of the raw materials. For simplicity, in the previous work we focussed on the electrical behaviour around the base of a single electrode, and prescribed the current in this electrode to be sinusoidal, with given amplitude. In this paper, we extend our previous analysis to include the full electrical system, modelled using an equivalent circuit system. In this way, we demonstrate how the two furnace-modelling approaches (on the fast and slow timescales) may be combined in a computationally efficient way. Our previously derived model for the arc resistance is based on the assumption that the dominant heat loss from the arc is by radiation (we will refer to this as the radiation model). Alternative arc models include the empirical Cassie and Mayr models, which are commonly used in the SAF literature. We compare these various arc models, explore the dependence of the solution of our model on the model parameters and compare our solutions with measurements from an operational silicon furnace. In particular, we show that only the radiation arc model has a rising current-voltage characteristic at high currents. Simulations of the model show that there is an upper limit on the length of the furnace arc, above which all the current bypasses the arc and flows through the surrounding material.
{"title":"Modelling alternating current effects in a submerged arc furnace","authors":"E. Luckins, James M. Oliver, C. Please, Benjamin M. Sloman, A. Valderhaug, R. V. Van Gorder","doi":"10.1093/imamat/hxac012","DOIUrl":"https://doi.org/10.1093/imamat/hxac012","url":null,"abstract":"\u0000 Modelling the production of silicon in a submerged arc furnace (SAF) requires accounting for the wide range of timescales of the different physical and chemical processes: the electric current which is used to heat the furnace varies over a timescale of around $10^{-2},$ s, whereas the flow and chemical consumption of the raw materials in the furnace occurs over several hours. Models for the silicon furnace generally either include only the fast-timescale, or only the slow-timescale processes. In a prior work, we developed a model incorporating effects on both the fast and slow timescales, and used a multiple-timescales analysis to homogenise the fast variations, deriving an averaged model for the slow evolution of the raw materials. For simplicity, in the previous work we focussed on the electrical behaviour around the base of a single electrode, and prescribed the current in this electrode to be sinusoidal, with given amplitude. In this paper, we extend our previous analysis to include the full electrical system, modelled using an equivalent circuit system. In this way, we demonstrate how the two furnace-modelling approaches (on the fast and slow timescales) may be combined in a computationally efficient way. Our previously derived model for the arc resistance is based on the assumption that the dominant heat loss from the arc is by radiation (we will refer to this as the radiation model). Alternative arc models include the empirical Cassie and Mayr models, which are commonly used in the SAF literature. We compare these various arc models, explore the dependence of the solution of our model on the model parameters and compare our solutions with measurements from an operational silicon furnace. In particular, we show that only the radiation arc model has a rising current-voltage characteristic at high currents. Simulations of the model show that there is an upper limit on the length of the furnace arc, above which all the current bypasses the arc and flows through the surrounding material.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":"1 1","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44163717","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We examine the analytic extension of solutions of linear, constant-coefficient initial-boundary value problems outside their spatial domain of definition. We use the Unified Transform Method or Method of Fokas, which gives a representation for solutions to half-line and finite-interval initial-boundary value problems as integrals of kernels with explicit spatial and temporal dependence. These solution representations are defined within the spatial domain of the problem. We obtain the extension of these representation formulae via Taylor series outside these spatial domains and find the extension of the initial condition that gives rise to a whole-line initial-value problem solved by the extended solution. In general, the extended initial condition is not differentiable or continuous unless the boundary and initial conditions satisfy compatibility conditions. We analyze dissipative and dispersive problems, and problems with continuous and discrete spatial variables.
{"title":"The analytic extension of solutions to initial-boundary value problems outside their domain of definition","authors":"Matthew Farkas, J. Cisneros, B. Deconinck","doi":"10.1093/imamat/hxad007","DOIUrl":"https://doi.org/10.1093/imamat/hxad007","url":null,"abstract":"\u0000 We examine the analytic extension of solutions of linear, constant-coefficient initial-boundary value problems outside their spatial domain of definition. We use the Unified Transform Method or Method of Fokas, which gives a representation for solutions to half-line and finite-interval initial-boundary value problems as integrals of kernels with explicit spatial and temporal dependence. These solution representations are defined within the spatial domain of the problem. We obtain the extension of these representation formulae via Taylor series outside these spatial domains and find the extension of the initial condition that gives rise to a whole-line initial-value problem solved by the extended solution. In general, the extended initial condition is not differentiable or continuous unless the boundary and initial conditions satisfy compatibility conditions. We analyze dissipative and dispersive problems, and problems with continuous and discrete spatial variables.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49047344","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We look at the periodic behaviour of the Earth’s glacial cycles and the transitions between different periodic states when either external parameters (such as $omega $) or internal parameters (such as $d$) are varied. We model this using the PP04 model of climate change. This is a forced discontinuous Filippov (non-smooth) dynamical system. When periodically forced this has coexisting periodic orbits. We find that the transitions in this system are mainly due to grazing events, leading to grazing bifurcations. An analysis of the grazing bifurcations is given and the impact of these on the domains of attraction and regions of existence of the periodic orbits is determined under various changes in the parameters of the system. Grazing transitions arise for general variations in the parameters (both internal and external) of the PP04 model. We find that the grazing transitions between the period orbits resemble those of the Mid-Pleistocene-Transition.
{"title":"Grazing bifurcations and transitions between periodic states of the PP04 model for the glacial cycle","authors":"Chris J Budd Kgomotso S. Morupisi","doi":"10.1093/imamat/hxac013","DOIUrl":"https://doi.org/10.1093/imamat/hxac013","url":null,"abstract":"\u0000 We look at the periodic behaviour of the Earth’s glacial cycles and the transitions between different periodic states when either external parameters (such as $omega $) or internal parameters (such as $d$) are varied. We model this using the PP04 model of climate change. This is a forced discontinuous Filippov (non-smooth) dynamical system. When periodically forced this has coexisting periodic orbits. We find that the transitions in this system are mainly due to grazing events, leading to grazing bifurcations. An analysis of the grazing bifurcations is given and the impact of these on the domains of attraction and regions of existence of the periodic orbits is determined under various changes in the parameters of the system. Grazing transitions arise for general variations in the parameters (both internal and external) of the PP04 model. We find that the grazing transitions between the period orbits resemble those of the Mid-Pleistocene-Transition.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-05-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44181526","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nastaran Naghshineh, W. Reinberger, N. Barlow, M. Samaha, S. J. Weinstein
We examine the power series solutions of two classical nonlinear ordinary differential equations of fluid mechanics that are mathematically related by their large-distance asymptotic behaviors in semi-infinite domains. The first problem is that of the “Sakiadis” boundary layer over a moving flat wall, for which no exact analytic solution has been put forward. The second problem is that of a static air–liquid meniscus with surface tension that intersects a flat wall at a given contact angle and limits to a flat pool away from the wall. For the latter problem, the exact analytic solution—given as distance from the wall as function of meniscus height—has long been known (Batchelor, 1967). Here, we provide an explicit solution as meniscus height vs. distance from the wall to elucidate structural similarities to the Sakiadis boundary layer. Although power series solutions are readily obtainable to the governing nonlinear ordinary differential equations, we show that—in both problems—the series diverge due to non-physical complex or negative real-valued singularities. In both cases, these singularities can be moved by expanding in exponential gauge functions motivated by their respective large distance asymptotic behaviors to enable series convergence over their full semi-infinite domains. For the Sakiadis problem, this not only provides a convergent Taylor series (and conjectured exact) solution to the ODE, but also a means to evaluate the wall shear parameter (and other properties) to within any desired precision. Although the nature of nonlinear ODEs precludes general conclusions, our results indicate that asymptotic behaviors can be useful when proposing variable transformations to overcome power series divergence. Sakiadis boundary layer; meniscus; asymptotic expansion; summation of series
{"title":"On the use of asymptotically motivated gauge functions to obtain convergent series solutions to nonlinear ODEs","authors":"Nastaran Naghshineh, W. Reinberger, N. Barlow, M. Samaha, S. J. Weinstein","doi":"10.1093/imamat/hxad006","DOIUrl":"https://doi.org/10.1093/imamat/hxad006","url":null,"abstract":"\u0000 We examine the power series solutions of two classical nonlinear ordinary differential equations of fluid mechanics that are mathematically related by their large-distance asymptotic behaviors in semi-infinite domains. The first problem is that of the “Sakiadis” boundary layer over a moving flat wall, for which no exact analytic solution has been put forward. The second problem is that of a static air–liquid meniscus with surface tension that intersects a flat wall at a given contact angle and limits to a flat pool away from the wall. For the latter problem, the exact analytic solution—given as distance from the wall as function of meniscus height—has long been known (Batchelor, 1967). Here, we provide an explicit solution as meniscus height vs. distance from the wall to elucidate structural similarities to the Sakiadis boundary layer. Although power series solutions are readily obtainable to the governing nonlinear ordinary differential equations, we show that—in both problems—the series diverge due to non-physical complex or negative real-valued singularities. In both cases, these singularities can be moved by expanding in exponential gauge functions motivated by their respective large distance asymptotic behaviors to enable series convergence over their full semi-infinite domains. For the Sakiadis problem, this not only provides a convergent Taylor series (and conjectured exact) solution to the ODE, but also a means to evaluate the wall shear parameter (and other properties) to within any desired precision. Although the nature of nonlinear ODEs precludes general conclusions, our results indicate that asymptotic behaviors can be useful when proposing variable transformations to overcome power series divergence. Sakiadis boundary layer; meniscus; asymptotic expansion; summation of series","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46709109","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper we consider the inverse scattering problem of determining the shape of a two-layered cavity with conductive boundary condition from sources and measurements placed on a curve inside the cavity. First, we show the well-posedness of the direct scattering problem by using the boundary integral equation method. Then, we prove that the factorization method can be applied to reconstruct the interface of the two-layered cavity from near-field data. Some numerical experiments are also presented to demonstrate the feasibility and effectiveness of the factorization method.
{"title":"The factorization method for inverse scattering by a two-layered cavity with conductive boundary condition","authors":"Jianguo Ye, G. Yan","doi":"10.1093/imamat/hxac005","DOIUrl":"https://doi.org/10.1093/imamat/hxac005","url":null,"abstract":"\u0000 In this paper we consider the inverse scattering problem of determining the shape of a two-layered cavity with conductive boundary condition from sources and measurements placed on a curve inside the cavity. First, we show the well-posedness of the direct scattering problem by using the boundary integral equation method. Then, we prove that the factorization method can be applied to reconstruct the interface of the two-layered cavity from near-field data. Some numerical experiments are also presented to demonstrate the feasibility and effectiveness of the factorization method.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44618289","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Topologically protected wave motion has attracted considerable research interest due to its chirality and potential applications in many applied fields. We construct quasi-traveling wave solutions to the two-dimensional Dirac equation with a domain wall mass in this work. It is known that the system admits exact and explicit traveling wave solutions, which are termed edge states if the interface is a straight line. By modifying such explicit solutions, we construct quasi-traveling-wave solutions if the interface is nearly straight. The approximate solutions in two scenarios are given. One is the circular edge with a large radius, and the second is a straight line edge with the slowly varying along the perpendicular direction. We show the quasi-traveling wave solutions are valid in a long lifespan by energy estimates. Numerical simulations are provided to support our analysis both qualitatively and quantitatively.
{"title":"Traveling edge states in massive Dirac equations along slowly varying edges","authors":"Pipi Hu, Peng Xie, Yi Zhu","doi":"10.1093/imamat/hxad015","DOIUrl":"https://doi.org/10.1093/imamat/hxad015","url":null,"abstract":"\u0000 Topologically protected wave motion has attracted considerable research interest due to its chirality and potential applications in many applied fields. We construct quasi-traveling wave solutions to the two-dimensional Dirac equation with a domain wall mass in this work. It is known that the system admits exact and explicit traveling wave solutions, which are termed edge states if the interface is a straight line. By modifying such explicit solutions, we construct quasi-traveling-wave solutions if the interface is nearly straight. The approximate solutions in two scenarios are given. One is the circular edge with a large radius, and the second is a straight line edge with the slowly varying along the perpendicular direction. We show the quasi-traveling wave solutions are valid in a long lifespan by energy estimates. Numerical simulations are provided to support our analysis both qualitatively and quantitatively.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2022-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49077028","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We study the homogenization of the Poisson equation with a reaction term and of the eigenvalue problem associated to the generator of multiscale Langevin dynamics. Our analysis extends the theory of two-scale convergence to the case of weighted Sobolev spaces in unbounded domains. We provide convergence results for the solution of the multiscale problems above to their homogenized surrogate. A series of numerical examples corroborate our analysis.
{"title":"Homogenization results for the generator of multiscale Langevin dynamics in weighted Sobolev spaces","authors":"Andrea Zanoni","doi":"10.1093/imamat/hxad003","DOIUrl":"https://doi.org/10.1093/imamat/hxad003","url":null,"abstract":"\u0000 We study the homogenization of the Poisson equation with a reaction term and of the eigenvalue problem associated to the generator of multiscale Langevin dynamics. Our analysis extends the theory of two-scale convergence to the case of weighted Sobolev spaces in unbounded domains. We provide convergence results for the solution of the multiscale problems above to their homogenized surrogate. A series of numerical examples corroborate our analysis.","PeriodicalId":56297,"journal":{"name":"IMA Journal of Applied Mathematics","volume":" ","pages":""},"PeriodicalIF":1.2,"publicationDate":"2021-12-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41566465","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We consider an initial-value problem based on a class of scalar nonlinear hyperbolic reaction–diffusion equations of the general form $$begin{align*} & u_{tautau}+u_{tau}=u_{{xx}}+varepsilon (F(u)+F(u)_{tau} ), end{align*}$$