Roy Baty - Academia.edu (original) (raw)
Papers by Roy Baty
Journal of Sound and Vibration, May 1, 1997
Bulletin of the American Physical Society, 2015
considers the group invariance properties of the inviscid compressible flow equations (Euler equa... more considers the group invariance properties of the inviscid compressible flow equations (Euler equations) under the assumptions of one-dimensional symmetry and a modified Tait equation of state (EOS) closure model. When written in terms of an adiabatic bulk modulus, a transformed version of these equations is found to be identical to that for an ideal gas EOS. As a result, the Lie group invariance structure of these equations and their subsequent reduction to a lower-order system is identical to the published results for the ideal gas case. Following the reduction of the Euler equations to a system of ordinary differential equations, a variety of elementary closed-form solutions are derived. These solutions are then used in conjunction with the Rankine-Hugoniot conditions to construct discontinuous shock wave and free surface solutions that are analogous to the classical Noh, Sedov, Guderley, and Hunter similarity solutions of the Euler equations for an ideal gas EOS. The versions of these problems for the modified Tait EOS are found to be semi-analytic in that a transcendental root extraction (and in some cases numerical integration of ordinary differential equations) enables solution of the relevant equations.
Quarterly Journal of Mechanics and Applied Mathematics, 2017
We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindri... more We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindrical or spherical shock wave from the ideal gas case to the case of flows with an arbitrary equation of state closure model, giving necessary conditions for the existence of a solution. The necessary condition is a thermodynamic one, namely that the adiabatic bulk modulus, KS, of the fluid be of the form pf (ρ) where p is pressure, ρ is mass density, and f is an arbitrary function. Although this condition has appeared in the literature before, here we give a more rigorous and extensive treatment. Of particular interest is our novel analysis of the governing ordinary differential equations (ODEs), which shows that, in general, the Guderley problem is always an eigenvalue problem. The need for an eigenvalue arises from basic shock stability principles-an interesting connection to the existing literature on the relationship between self-similarity of the second kind and stability. We also investigate a special case, usually neglected by previous authors, where assuming constant shock velocity yields a reduction to ODEs for every material, but those ODEs never have a bounded, differentiable solution. This theoretical work is motivated by the need for more realistic test problems in the verification of inviscid compressible flow codes that simulate flows in a variety of non-ideal gas materials.
Bulletin of the American Physical Society, 2016
The Quarterly Journal of Mechanics and Applied Mathematics, 2019
As modern hydrodynamic codes increase in sophistication, the availability of realistic test probl... more As modern hydrodynamic codes increase in sophistication, the availability of realistic test problems becomes increasingly important. In gas dynamics, one common unrealistic aspect of most test problems is the ideal gas assumption, which is unsuited to many real applications, especially those involving high pressure and speed metal deformation. Our work considers the collapsing cavity and converging shock test problems, showing to what extent the ideal gas assumption can be removed from their specification. It is found that while most materials simply do not admit simple (that is scaling) solutions in this context, there are infinite-dimensional families of materials which do admit such solutions. We characterize such materials, derive the appropriate ordinary differential equations and analyze the associated nonlinear eigenvalue problem. It is shown that there is an inherent tension between boundedness of the solution, boundedness of its derivatives and the entropy condition. The sp...
This article applies nonstandard analysis to derive jump conditions for one-dimensional, divergin... more This article applies nonstandard analysis to derive jump conditions for one-dimensional, diverging, magnetogasdynamic shock waves emerging on the surface of a star. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions for the flow parameters occur smoothly across this interval. Predistributions of the Heaviside function and the Dirac delta measure are used to model the flow variables across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the flow parameters. It is shown here that the equations modeling a family of magnetogasdynamic shock waves yield products of generalized functions that may be analyzed consistently using nonstandard predistributions.
We investigate the inviscid compressible flow (Euler) equations constrained by an "isentropic" eq... more We investigate the inviscid compressible flow (Euler) equations constrained by an "isentropic" equation of state (EOS), whose functional form in pressure is an arbitrary function of density alone. Under the aforementioned condition, we interrogate using symmetry methods the scale-invariance of the homentropic inviscid Euler equations. We find that under general conditions, we can reduce the inviscid Euler equations into a system of two coupled ordinary differential equations. To exemplify the utility of these results, we formulate two example scale-invariant, self-similar solutions. The first example includes a shock-free expanding bubble scenario, featuring a modified Tait EOS. The second example features the classical Noh problem, coupled to an arbitrary isentropic EOS. In this case, in order to satisfy the conditions set forth in the classical Noh problem, we find that the solution for the flow is given by a transcendental algebraic equation in the shocked density.
International Journal of Aeroacoustics
This article applies nonstandard analysis to study the generalized solutions of entropy and energ... more This article applies nonstandard analysis to study the generalized solutions of entropy and energy across one-dimensional shock waves in a compressible, inviscid, perfect gas. Nonstandard analysis is an area of modern mathematics that studies number systems that contain both infinitely large and infinitely small numbers. For an inviscid shock wave, it is assumed that the shock thickness occurs on an infinitesimal interval and that the jump functions for the field variables are smoothly defined on this interval. A weak converse to the existence of the entropy peak is derived and discussed. Generalized solutions of the Euler equations for entropy and energy are then derived for both theoretical and realistic normalized velocity profiles.
Physics of Fluids
The problem of a one-dimensional (1D) cylindrically or spherically symmetric shock wave convergin... more The problem of a one-dimensional (1D) cylindrically or spherically symmetric shock wave converging into an inviscid, ideal gas was first investigated by Guderley[Starke kugelige und zylinrische verdichtungsstosse in der nahe des kugelmitterpunktes bzw. Der zylinderachse," Luftfahrtforschung 19, 302 (1942)]. In the time since, many authors have discussed the practical notion of how Guderley-like flows might be generated. One candidate is a constant velocity, converging "cylindrical or spherical piston," giving rise to a converging shock wave in the spirit of its classical, planar counterpart. A limitation of pre-existing analyses along these lines is the restriction to flows in materials described by an ideal gas equation of state (EOS) constitutive law. This choice is of course necessary for the direct comparison with the classical Guderley solution, which also features an ideal gas EOS. However, the ideal gas EOS is limited in its utility in describing a wide variety of physical phenomena and, in particular, the shock compression of solid materials. This work is thus intended to provide an extension of previous work to a nonideal EOS. The stiff gas EOS is chosen as a logical starting point due to not only its close resemblance to the ideal gas law but also its relevance to the shock compression of various liquid and solid materials. Using this choice of EOS, the solution of a 1D planar piston problem is constructed and subsequently used as the lowest order term in a quasi-self-similar series expansion intended to capture both curvilinear and nonideal EOS effects. The solution associated with this procedure provides correction terms to the 1D planar solution so that the expected accelerating shock trajectory and nontrivially varying state variable profiles can be obtained. This solution is further examined in the limit as the converging shock wave approaches the 1D curvilinear origin. Given the stiff gas EOS is not otherwise expected to admit a Guderley-like solution when coupled to the inviscid Euler equations, this work thus provides the semianalytical limiting behavior of a flow that cannot be otherwise captured using self-similar analysis.
Journal of Mathematical Physics
A brief review of the theory of exterior differential systems and isovector symmetry analysis met... more A brief review of the theory of exterior differential systems and isovector symmetry analysis methods is presented in the context of the one-dimensional inviscid compressible flow equations. These equations are formulated as an exterior differential system with equation of state (EOS) closure provided in terms of an adiabatic bulk modulus. The scaling symmetry generators-and corresponding EOS constraints-otherwise appearing in the existing literature are recovered through the application of and invariance under Lie derivative dragging operations.
International Journal of Aeroacoustics
Journal of Geometric Mechanics
Physics of Fluids
Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equation... more Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Grüneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Grüneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie-Grüneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Grüneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Grüneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical, Guderley-like imploding shock solution in a metal sphere, and to determine if and when the solution may be valid for the underlying Mie-Grüneisen EOS.
The Quarterly Journal of Mechanics and Applied Mathematics
Cavity-expansion approximations are widely-used in the study of penetration mechanics and indenta... more Cavity-expansion approximations are widely-used in the study of penetration mechanics and indentation phenomena. We apply the isovector method to a well-established model in the literature for elastic-plastic cavity-expansion to systematically demonstrate the existence of Lie symmetries corresponding to scale-invariant solutions. We use the symmetries obtained from the equations of motion to determine compatible auxiliary conditions describing the cavity wall trajectory and the elastic-plastic material interface. The admissible conditions are then compared with specific similarity solutions in the literature.
Journal of Sound and Vibration, May 1, 1997
Bulletin of the American Physical Society, 2015
considers the group invariance properties of the inviscid compressible flow equations (Euler equa... more considers the group invariance properties of the inviscid compressible flow equations (Euler equations) under the assumptions of one-dimensional symmetry and a modified Tait equation of state (EOS) closure model. When written in terms of an adiabatic bulk modulus, a transformed version of these equations is found to be identical to that for an ideal gas EOS. As a result, the Lie group invariance structure of these equations and their subsequent reduction to a lower-order system is identical to the published results for the ideal gas case. Following the reduction of the Euler equations to a system of ordinary differential equations, a variety of elementary closed-form solutions are derived. These solutions are then used in conjunction with the Rankine-Hugoniot conditions to construct discontinuous shock wave and free surface solutions that are analogous to the classical Noh, Sedov, Guderley, and Hunter similarity solutions of the Euler equations for an ideal gas EOS. The versions of these problems for the modified Tait EOS are found to be semi-analytic in that a transcendental root extraction (and in some cases numerical integration of ordinary differential equations) enables solution of the relevant equations.
Quarterly Journal of Mechanics and Applied Mathematics, 2017
We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindri... more We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindrical or spherical shock wave from the ideal gas case to the case of flows with an arbitrary equation of state closure model, giving necessary conditions for the existence of a solution. The necessary condition is a thermodynamic one, namely that the adiabatic bulk modulus, KS, of the fluid be of the form pf (ρ) where p is pressure, ρ is mass density, and f is an arbitrary function. Although this condition has appeared in the literature before, here we give a more rigorous and extensive treatment. Of particular interest is our novel analysis of the governing ordinary differential equations (ODEs), which shows that, in general, the Guderley problem is always an eigenvalue problem. The need for an eigenvalue arises from basic shock stability principles-an interesting connection to the existing literature on the relationship between self-similarity of the second kind and stability. We also investigate a special case, usually neglected by previous authors, where assuming constant shock velocity yields a reduction to ODEs for every material, but those ODEs never have a bounded, differentiable solution. This theoretical work is motivated by the need for more realistic test problems in the verification of inviscid compressible flow codes that simulate flows in a variety of non-ideal gas materials.
Bulletin of the American Physical Society, 2016
The Quarterly Journal of Mechanics and Applied Mathematics, 2019
As modern hydrodynamic codes increase in sophistication, the availability of realistic test probl... more As modern hydrodynamic codes increase in sophistication, the availability of realistic test problems becomes increasingly important. In gas dynamics, one common unrealistic aspect of most test problems is the ideal gas assumption, which is unsuited to many real applications, especially those involving high pressure and speed metal deformation. Our work considers the collapsing cavity and converging shock test problems, showing to what extent the ideal gas assumption can be removed from their specification. It is found that while most materials simply do not admit simple (that is scaling) solutions in this context, there are infinite-dimensional families of materials which do admit such solutions. We characterize such materials, derive the appropriate ordinary differential equations and analyze the associated nonlinear eigenvalue problem. It is shown that there is an inherent tension between boundedness of the solution, boundedness of its derivatives and the entropy condition. The sp...
This article applies nonstandard analysis to derive jump conditions for one-dimensional, divergin... more This article applies nonstandard analysis to derive jump conditions for one-dimensional, diverging, magnetogasdynamic shock waves emerging on the surface of a star. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions for the flow parameters occur smoothly across this interval. Predistributions of the Heaviside function and the Dirac delta measure are used to model the flow variables across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the flow parameters. It is shown here that the equations modeling a family of magnetogasdynamic shock waves yield products of generalized functions that may be analyzed consistently using nonstandard predistributions.
We investigate the inviscid compressible flow (Euler) equations constrained by an "isentropic" eq... more We investigate the inviscid compressible flow (Euler) equations constrained by an "isentropic" equation of state (EOS), whose functional form in pressure is an arbitrary function of density alone. Under the aforementioned condition, we interrogate using symmetry methods the scale-invariance of the homentropic inviscid Euler equations. We find that under general conditions, we can reduce the inviscid Euler equations into a system of two coupled ordinary differential equations. To exemplify the utility of these results, we formulate two example scale-invariant, self-similar solutions. The first example includes a shock-free expanding bubble scenario, featuring a modified Tait EOS. The second example features the classical Noh problem, coupled to an arbitrary isentropic EOS. In this case, in order to satisfy the conditions set forth in the classical Noh problem, we find that the solution for the flow is given by a transcendental algebraic equation in the shocked density.
International Journal of Aeroacoustics
This article applies nonstandard analysis to study the generalized solutions of entropy and energ... more This article applies nonstandard analysis to study the generalized solutions of entropy and energy across one-dimensional shock waves in a compressible, inviscid, perfect gas. Nonstandard analysis is an area of modern mathematics that studies number systems that contain both infinitely large and infinitely small numbers. For an inviscid shock wave, it is assumed that the shock thickness occurs on an infinitesimal interval and that the jump functions for the field variables are smoothly defined on this interval. A weak converse to the existence of the entropy peak is derived and discussed. Generalized solutions of the Euler equations for entropy and energy are then derived for both theoretical and realistic normalized velocity profiles.
Physics of Fluids
The problem of a one-dimensional (1D) cylindrically or spherically symmetric shock wave convergin... more The problem of a one-dimensional (1D) cylindrically or spherically symmetric shock wave converging into an inviscid, ideal gas was first investigated by Guderley[Starke kugelige und zylinrische verdichtungsstosse in der nahe des kugelmitterpunktes bzw. Der zylinderachse," Luftfahrtforschung 19, 302 (1942)]. In the time since, many authors have discussed the practical notion of how Guderley-like flows might be generated. One candidate is a constant velocity, converging "cylindrical or spherical piston," giving rise to a converging shock wave in the spirit of its classical, planar counterpart. A limitation of pre-existing analyses along these lines is the restriction to flows in materials described by an ideal gas equation of state (EOS) constitutive law. This choice is of course necessary for the direct comparison with the classical Guderley solution, which also features an ideal gas EOS. However, the ideal gas EOS is limited in its utility in describing a wide variety of physical phenomena and, in particular, the shock compression of solid materials. This work is thus intended to provide an extension of previous work to a nonideal EOS. The stiff gas EOS is chosen as a logical starting point due to not only its close resemblance to the ideal gas law but also its relevance to the shock compression of various liquid and solid materials. Using this choice of EOS, the solution of a 1D planar piston problem is constructed and subsequently used as the lowest order term in a quasi-self-similar series expansion intended to capture both curvilinear and nonideal EOS effects. The solution associated with this procedure provides correction terms to the 1D planar solution so that the expected accelerating shock trajectory and nontrivially varying state variable profiles can be obtained. This solution is further examined in the limit as the converging shock wave approaches the 1D curvilinear origin. Given the stiff gas EOS is not otherwise expected to admit a Guderley-like solution when coupled to the inviscid Euler equations, this work thus provides the semianalytical limiting behavior of a flow that cannot be otherwise captured using self-similar analysis.
Journal of Mathematical Physics
A brief review of the theory of exterior differential systems and isovector symmetry analysis met... more A brief review of the theory of exterior differential systems and isovector symmetry analysis methods is presented in the context of the one-dimensional inviscid compressible flow equations. These equations are formulated as an exterior differential system with equation of state (EOS) closure provided in terms of an adiabatic bulk modulus. The scaling symmetry generators-and corresponding EOS constraints-otherwise appearing in the existing literature are recovered through the application of and invariance under Lie derivative dragging operations.
International Journal of Aeroacoustics
Journal of Geometric Mechanics
Physics of Fluids
Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equation... more Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Grüneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Grüneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie-Grüneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Grüneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Grüneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical, Guderley-like imploding shock solution in a metal sphere, and to determine if and when the solution may be valid for the underlying Mie-Grüneisen EOS.
The Quarterly Journal of Mechanics and Applied Mathematics
Cavity-expansion approximations are widely-used in the study of penetration mechanics and indenta... more Cavity-expansion approximations are widely-used in the study of penetration mechanics and indentation phenomena. We apply the isovector method to a well-established model in the literature for elastic-plastic cavity-expansion to systematically demonstrate the existence of Lie symmetries corresponding to scale-invariant solutions. We use the symmetries obtained from the equations of motion to determine compatible auxiliary conditions describing the cavity wall trajectory and the elastic-plastic material interface. The admissible conditions are then compared with specific similarity solutions in the literature.